| 89 | |
| 90 | |
| 91 | void dMultiply2 (dReal *A, const dReal *B, const dReal *C, int p, int q, int r) |
| 92 | { |
| 93 | int i,j,k,z,rpad,qskip; |
| 94 | dReal sum; |
| 95 | const dReal *bb,*cc; |
| 96 | dAASSERT (A && B && C && p>0 && q>0 && r>0); |
| 97 | rpad = dPAD(r) - r; |
| 98 | qskip = dPAD(q); |
| 99 | bb = B; |
| 100 | for (i=p; i; i--) { |
| 101 | cc = C; |
| 102 | for (j=r; j; j--) { |
| 103 | z = 0; |
| 104 | sum = 0; |
| 105 | for (k=q; k; k--,z++) sum += bb[z] * cc[z]; |
| 106 | *(A++) = sum; |
| 107 | cc += qskip; |
| 108 | } |
| 109 | A += rpad; |
| 110 | bb += qskip; |
| 111 | } |
| 112 | } |
| 113 | |
| 114 | |
| 115 | int dFactorCholesky (dReal *A, int n) |
no outgoing calls
no test coverage detected