source:
CIVL/mods/dev.civl.com/examples/compare/CholeskyDecomposition/cholsl.c@
cb4d4f4
| Last change on this file since cb4d4f4 was aad342c, checked in by , 3 years ago | |
|---|---|
|
|
| File size: 304 bytes | |
| Rev | Line | |
|---|---|---|
| [9765af0] | 1 | |
| 2 | void cholsl(float **a, int n, float p[], float b[], float x[]) { | |
| 3 | int i, k; | |
| 4 | float sum; | |
| 5 | ||
| 6 | for (i=0; i<=n; i++) { | |
| 7 | for (sum=b[i], k=i-1; k>=1; k--) sum -= a[i][k]*x[k]; | |
| 8 | x[i]=sum/p[i]; | |
| 9 | } | |
| 10 | for (i=n;i>=1;i--) { | |
| 11 | for (sum=x[i], k=i+1; k<=n; k++) sum -= a[k][i]*x[k]; | |
| 12 | x[i]=sum/p[i]; | |
| 13 | } | |
| 14 | } |
Note:
See TracBrowser
for help on using the repository browser.
