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