sprsax.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 16 行

void sprsax(sa,ija,x,b,n)
float b[],sa[],x[];
unsigned long ija[],n;
{
	void nrerror();
	unsigned long i,k;

	if (ija[1] != n+2) nrerror("sprsax: mismatched vector and matrix");
	for (i=1;i<=n;i++) {
		b[i]=sa[i]*x[i];
		for (k=ija[i];k<=ija[i+1]-1;k++)
			b[i] += sa[k]*x[ija[k]];

		}
}