// Determine the largest abs numerical value for partial pivoting
Index diagind = iperm_c(jcol); // diagonal index
RealScalar pivmax = 0.0;
RealScalar pivmax = -1.0;
Index pivptr = nsupc;
Index diag = emptyIdxLU;
RealScalar rtemp;
if (lsub_ptr[isub] == diagind) diag = isub;
}
// Test for singularity (pivoting fails)
if ( pivmax == -1.0 ) {
perm_r[diagind] = jcol;
return (jcol+1);
// Test for singularity
if ( pivmax == 0.0 ) {
pivrow = lsub_ptr[pivptr];