svdvar.c

来自「数值算法」· C语言 代码 · 共 22 行

void svdvar(v,ma,w,cvm)
float **v,w[],**cvm;
int ma;
{
	int k,j,i;
	float sum,*wti,*vector();
	void free_vector();

	wti=vector(1,ma);
	for (i=1;i<=ma;i++) {
		wti[i]=0.0;
		if (w[i]) wti[i]=1.0/(w[i]*w[i]);
	}
	for (i=1;i<=ma;i++) {
		for (j=1;j<=i;j++) {
			for (sum=0.0,k=1;k<=ma;k++) sum += v[i][k]*v[j][k]*wti[k];
			cvm[j][i]=cvm[i][j]=sum;
		}
	}
	free_vector(wti,1,ma);
}