相关代码 |
|
double spdmtrace(M1,M2,n)/* trace of M1*M2, M1 double, M2 short */ register double *M1; register short *M2; register int n; { register int i; register double trace; double spdot(); trace=0; for(i=0;i trace+=spdot(M1,1,M2,n,n); return(trace); }