/*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));
}
}