Metropoli BBS
VIEWER: invertc.h MODE: TEXT (ASCII)
/*Copyright (C) 1993, 1996 by Thomas Glen Smith.  All Rights Reserved.*/
/* invertc.h APL2 V1.0.0 **********************************************
* Included in invert to make elements 0 in column i below row i.      *
**********************************************************************/
		if (i<(rows-1))
			for (j=i+1;j<rows;j++)
				if ((w=*(t=(m+j*rows+i)))!=0.0) {
					*t=0.0;
					t=m+j*rows;
					v=m+n;
					for (k=i+1;k<rows;k++) {
						s=t+k;
						*s=*s/w-*(v+k);
					}
					t=z+j*rows;
					v=z+n;
					for (k=0;k<rows;k++) {
						s=t+k;
						*s=*s/w-*(v+k);
					}
				}
[ RETURN TO DIRECTORY ]