| 72 | |
| 73 | |
| 74 | void dMultiply1 (dReal *A, const dReal *B, const dReal *C, int p, int q, int r) |
| 75 | { |
| 76 | int i,j,k,pskip,rskip; |
| 77 | dReal sum; |
| 78 | dAASSERT (A && B && C && p>0 && q>0 && r>0); |
| 79 | pskip = dPAD(p); |
| 80 | rskip = dPAD(r); |
| 81 | for (i=0; i<p; i++) { |
| 82 | for (j=0; j<r; j++) { |
| 83 | sum = 0; |
| 84 | for (k=0; k<q; k++) sum += B[i+k*pskip] * C[j+k*rskip]; |
| 85 | A[i*rskip+j] = sum; |
| 86 | } |
| 87 | } |
| 88 | } |
| 89 | |
| 90 | |
| 91 | void dMultiply2 (dReal *A, const dReal *B, const dReal *C, int p, int q, int r) |
no outgoing calls
no test coverage detected