Metropoli BBS
VIEWER: invertd.h MODE: TEXT (ASCII)
/*Copyright (C) 1993, 1996 by Thomas Glen Smith.  All Rights Reserved.*/
/* invertd.h APL2 V1.0.0 **********************************************
* Included in invert to make elements 0 in column i above row i.      *
**********************************************************************/
		/* now make all elements in this column above the diagonal 0 */
		if (i>0)
			for (j=i-1;j>=0;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 ]