1207 …MatScalar       *pc, *pc1, *pc2, *pc3, *pc4, mul1, mul2, mul3, mul4, *pv, *rtmp1, *rtmp2, *rtmp3, …  in MatLUFactorNumeric_SeqAIJ_Inode()  local
1403           pc1 = rtmp1 + row;  in MatLUFactorNumeric_SeqAIJ_Inode()
1405           if (*pc1 != 0.0 || *pc2 != 0.0) {  in MatLUFactorNumeric_SeqAIJ_Inode()
1407             mul1 = *pc1 * (*pv);  in MatLUFactorNumeric_SeqAIJ_Inode()
1409             *pc1 = mul1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1428         pc1 = b->a + bi[i];  in MatLUFactorNumeric_SeqAIJ_Inode()
1433           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1434           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1437         pc1 = b->a + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1442           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1443           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1450         pc1  = b->a + bdiag[i]; /* Mark diagonal */  in MatLUFactorNumeric_SeqAIJ_Inode()
1451         *pc1 = 1.0 / sctx.pv;  in MatLUFactorNumeric_SeqAIJ_Inode()
1456           mul1 = (*pc2) * (*pc1);             /* *pc1=diag[i] is inverted! */  in MatLUFactorNumeric_SeqAIJ_Inode()
1540           pc1 = rtmp1 + row;  in MatLUFactorNumeric_SeqAIJ_Inode()
1543           if (*pc1 != 0.0 || *pc2 != 0.0 || *pc3 != 0.0) {  in MatLUFactorNumeric_SeqAIJ_Inode()
1545             mul1 = *pc1 * (*pv);  in MatLUFactorNumeric_SeqAIJ_Inode()
1548             *pc1 = mul1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1569         pc1 = b->a + bi[i];  in MatLUFactorNumeric_SeqAIJ_Inode()
1574           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1575           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1578         pc1 = b->a + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1583           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1584           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1591         pc1  = b->a + bdiag[i]; /* Mark diag[i] */  in MatLUFactorNumeric_SeqAIJ_Inode()
1592         *pc1 = 1.0 / sctx.pv;  in MatLUFactorNumeric_SeqAIJ_Inode()
1598           mul2 = (*pc2) * (*pc1);  in MatLUFactorNumeric_SeqAIJ_Inode()
1600           mul3 = (*pc3) * (*pc1);  in MatLUFactorNumeric_SeqAIJ_Inode()
1731           pc1 = rtmp1 + row;  in MatLUFactorNumeric_SeqAIJ_Inode()
1735           if (*pc1 != 0.0 || *pc2 != 0.0 || *pc3 != 0.0 || *pc4 != 0.0) {  in MatLUFactorNumeric_SeqAIJ_Inode()
1737             mul1 = *pc1 * (*pv);  in MatLUFactorNumeric_SeqAIJ_Inode()
1741             *pc1 = mul1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1764         pc1 = b->a + bi[i];  in MatLUFactorNumeric_SeqAIJ_Inode()
1769           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1770           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1773         pc1 = b->a + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ_Inode()
1778           pc1[j] = rtmp1[col];  in MatLUFactorNumeric_SeqAIJ_Inode()
1779           rs += PetscAbsScalar(pc1[j]);  in MatLUFactorNumeric_SeqAIJ_Inode()
1786         pc1  = b->a + bdiag[i]; /* Mark diag[i] */  in MatLUFactorNumeric_SeqAIJ_Inode()
1787         *pc1 = 1.0 / sctx.pv;  in MatLUFactorNumeric_SeqAIJ_Inode()
1794           mul2 = (*pc2) * (*pc1);  in MatLUFactorNumeric_SeqAIJ_Inode()
1796           mul3 = (*pc3) * (*pc1);  in MatLUFactorNumeric_SeqAIJ_Inode()
1798           mul4 = (*pc4) * (*pc1);  in MatLUFactorNumeric_SeqAIJ_Inode()