atimes.c

适合大型数值计算代码 现在网络上已经找不到了 购买需要20$

extern unsigned long ija[];
extern double sa[];

void atimes(unsigned long n, double x[], double r[], int itrnsp)
{
	void dsprsax(double sa[], unsigned long ija[], double x[], double b[],
		unsigned long n);
	void dsprstx(double sa[], unsigned long ija[], double x[], double b[],
		unsigned long n);

	if (itrnsp) dsprstx(sa,ija,x,r,n);
	else dsprsax(sa,ija,x,r,n);
}