banmul.c

来自「Numerical Recipes Software 提供的算法子程序集」· C语言 代码 · 共 18 行

#include "nrutil.h"

void banmul(a,n,m1,m2,x,b)
float **a,b[],x[];
int m1,m2;
unsigned long n;
{
	unsigned long i,j,k,tmploop;

	for (i=1;i<=n;i++) {
		k=i-m1-1;
		tmploop=LMIN(m1+m2+1,n-k);
		b[i]=0.0;
		for (j=LMAX(1,1-k);j<=tmploop;j++) b[i] += a[i][j]*x[j+k];
	}
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */