1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
// This is mul/mbl/mbl_matrix_products.h
#ifndef mbl_matrix_product_h
#define mbl_matrix_product_h
//:
// \file
// \author Tim Cootes
// \date 25-Apr-2001
// \brief Various specialised versions of matrix product operations
#include <vnl/vnl_fwd.h>
//: Compute product AB = A * B
void mbl_matrix_product(vnl_matrix<double>& AB, const vnl_matrix<double>& A,
const vnl_matrix<double>& B);
//: Compute product ABt = A * B.transpose()
void mbl_matrix_product_a_bt(vnl_matrix<double>& ABt, const vnl_matrix<double>& A,
const vnl_matrix<double>& B);
//: Compute ABt = A * B.transpose() using only first n_cols columns of A and B
void mbl_matrix_product_a_bt(vnl_matrix<double>& ABt,
const vnl_matrix<double>& A,
const vnl_matrix<double>& B,
int n_cols);
//: Compute AAt = A * A.transpose(), using only first nr x nc partition of A
// Uses symmetry of result to improve speed
void mbl_matrix_product_a_at(vnl_matrix<double>& AAt,
const vnl_matrix<double>& A,
unsigned nr, unsigned nc);
//: Compute AAt = A * A.transpose()
// Uses symmetry of result to improve speed
void mbl_matrix_product_a_at(vnl_matrix<double>& AAt,
const vnl_matrix<double>& A);
//: Compute product AtB = A.transpose() * B
// Result is of size A.cols() x B.cols()
void mbl_matrix_product_at_b(vnl_matrix<double>& AtB, const vnl_matrix<double>& A,
const vnl_matrix<double>& B);
//: Compute AtB = A.transpose() * B , using first ncols_a cols of A
// Result is of size ncols_a x B.cols()
void mbl_matrix_product_at_b(vnl_matrix<double>& AtB, const vnl_matrix<double>& A,
const vnl_matrix<double>& B, int ncols_a);
//: Compute product AtA = A.transpose() * A using nc cols of A
void mbl_matrix_product_at_a(vnl_matrix<double>& AtA,
const vnl_matrix<double>& A,
unsigned nc);
//: Compute product AtA = A.transpose() * A
void mbl_matrix_product_at_a(vnl_matrix<double>& AtA,
const vnl_matrix<double>& A);
//: ADB = A * D * B where D is diagonal with elements d
// A, d, and B should be compatible sizes, ADB will
// be resized to fit.
void mbl_matrix_product_adb(vnl_matrix<double>& ADB,
const vnl_matrix<double>& A,
const vnl_vector<double>& d,
const vnl_matrix<double>& B);
#endif // mbl_matrix_product_h
|
