rtkcmn.c
lsq中的矩阵运算
matmul
- 函数调用
matmul("NN",n,1,m,1.0,A,y,0.0,Ay); /* Ay=A*y */
- 函数定义
extern void matmul(const char *tr, int n, int k, int m, double alpha,
const double *A, const double *B, double beta, double *C)
- 函数声明
/* multiply matrix (wrapper of blas dgemm) ----------------*
multiply matrix by matrix (C=alpha*A*B+beta*C)
* args :
* char *tr I transpose flags("N":normal,"T":transpose)
* int n,k,m I size of (transposed) matrix A,B
* double alpha I alpha
* double *A,*B I (transposed) matrix A (n x m), B (m x k)
* double beta I beta
* double *C IO matrix C (n x k)* return : none
* -------------------------------------------------*/