if( j == ks -1 ) val = b.A[j][0];
else if( j > ks - 1 ) {
int jm = Math.min(n-1, j);
DoubleVector u = getColumn(ks,jm,p[j]);
DoubleVector v = b.getColumn(ks,jm,0);
val = v.innerProduct(u) / u.norm2();
}
else { // ks > j
for( k = j+1; k < ks; k++ ) // make a copy of A[j][]
xxx[k] = A[j][p[k]];
val = b.A[j][0];