cholsl.c

来自「统计模式识别算法包」· C语言 代码 · 共 28 行

/* This was taken from the new version of the book */

#include <math.h>

void cholsl(a,n,p,b,x)
float **a,*p,*b,*x;
int n;

{

int i,k;
float sum;


for(i=1;i<=n;i++)
	{
	for(sum=b[i],k=i-1;k>=1;k--)
		sum -= a[i][k]*x[k];
	x[i]=sum/p[i];
	}
for(i=n;i>=1;i--)
	{
	for(sum=x[i],k=i+1;k<=n;k++)
		sum -= a[k][i]*x[k];
	x[i]=sum/p[i];
	}
}