INNER CODE UNIT · JavaScript
big
lh3/yak · misc/test-qv.js:18
big = Math.abs(a[j][k]);
ir = j; ic = k;
}
} else if (ipiv[k] > 1) return -2; // singular matrix
}
}
}
++ipiv[ic];
if (ir != ic) {
for (l = 0; l < n; ++l) tmp = a[ir][l], a[ir][l] = a[ic][l], a[ic][l] = tmp;
if (b) for (l = 0; l < m; ++l) tmp = b[ir][l], b[ir][l] = b[ic][l], b[ic][l] = tmp;
}
xr[i] = ir; xc[i] = ic;
if (a[ic][ic] == 0) return -3; // singular matrix
var pivinv = 1. / a[ic][ic];
a[ic][ic] = 1.;
for (l = 0; l < n; ++l) a[ic][l] *= pivinv;
if (b) for (l = 0; l < m; ++l) b[ic][l] *= pivinv;