227   MatScalar       *rtmp, *pc, multiplier, *pv;  in MatLUFactorNumeric_SeqAIJ()  local
262   PetscCall(PetscMalloc1(n + 1, &rtmp));  in MatLUFactorNumeric_SeqAIJ()
272       for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0;  in MatLUFactorNumeric_SeqAIJ()
277       for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0;  in MatLUFactorNumeric_SeqAIJ()
283       for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j];  in MatLUFactorNumeric_SeqAIJ()
285       rtmp[i] += sctx.shift_amount; /* shift the diagonal of the matrix */  in MatLUFactorNumeric_SeqAIJ()
292         pc = rtmp + row;  in MatLUFactorNumeric_SeqAIJ()
302           for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j];  in MatLUFactorNumeric_SeqAIJ()
315         pv[j] = rtmp[pj[j]];  in MatLUFactorNumeric_SeqAIJ()
324         pv[j] = rtmp[pj[j]];  in MatLUFactorNumeric_SeqAIJ()
329       sctx.pv = rtmp[i];  in MatLUFactorNumeric_SeqAIJ()
332       rtmp[i] = sctx.pv;        /* sctx.pv might be updated in the case of MAT_SHIFT_INBLOCKS */  in MatLUFactorNumeric_SeqAIJ()
336       *pv = 1.0 / rtmp[i];  in MatLUFactorNumeric_SeqAIJ()
357   PetscCall(PetscFree(rtmp));  in MatLUFactorNumeric_SeqAIJ()
406   MatScalar       *pv, *rtmp, *pc, multiplier, d;  in MatLUFactorNumeric_SeqAIJ_inplace()  local
441   PetscCall(PetscMalloc1(n + 1, &rtmp));  in MatLUFactorNumeric_SeqAIJ_inplace()
449       for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0;  in MatLUFactorNumeric_SeqAIJ_inplace()
455       for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j];  in MatLUFactorNumeric_SeqAIJ_inplace()
456       rtmp[ics[r[i]]] += sctx.shift_amount; /* shift the diagonal of the matrix */  in MatLUFactorNumeric_SeqAIJ_inplace()
460         pc = rtmp + row;  in MatLUFactorNumeric_SeqAIJ_inplace()
467           for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j];  in MatLUFactorNumeric_SeqAIJ_inplace()
479         pv[j] = rtmp[pj[j]];  in MatLUFactorNumeric_SeqAIJ_inplace()
506   PetscCall(PetscFree(rtmp));  in MatLUFactorNumeric_SeqAIJ_inplace()
563   PetscScalar    *rtmp, *pc, multiplier, d;  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()  local
600   PetscCall(PetscMalloc1(n + 1, &rtmp));  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
601   PetscCall(PetscArrayzero(rtmp, n + 1));  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
646         rtmp[ajtmp[j]] = v[j];  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
649       rtmp[r[i]] += sctx.shift_amount; /* shift the diagonal of the matrix */  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
653         pc = rtmp + row;  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
661           for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j];  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
674         pv[j] = rtmp[pj[j]];  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
702   PetscCall(PetscFree(rtmp));  in MatLUFactorNumeric_SeqAIJ_InplaceWithPerm()
1713   MatScalar       *rtmp, *ba = b->a, *bval, dk, uikdi;  in MatCholeskyFactorNumeric_SeqAIJ()  local
1751   PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &c2r));  in MatCholeskyFactorNumeric_SeqAIJ()
1763       for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0;  in MatCholeskyFactorNumeric_SeqAIJ()
1772           rtmp[col] = aa[j];  in MatCholeskyFactorNumeric_SeqAIJ()
1777       rtmp[k] += sctx.shift_amount; /* shift the diagonal of the matrix */  in MatCholeskyFactorNumeric_SeqAIJ()
1780       dk = rtmp[k];  in MatCholeskyFactorNumeric_SeqAIJ()
1796           for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j];  in MatCholeskyFactorNumeric_SeqAIJ()
1813           ba[j] = rtmp[col];  in MatCholeskyFactorNumeric_SeqAIJ()
1835   PetscCall(PetscFree3(rtmp, il, c2r));  in MatCholeskyFactorNumeric_SeqAIJ()
1880   MatScalar       *rtmp, *ba = b->a, *bval, dk, uikdi;  in MatCholeskyFactorNumeric_SeqAIJ_inplace()  local
1915   PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &jl));  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1927       for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0;  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1936           rtmp[col] = aa[j];  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1941       if (sctx.nshift) rtmp[k] += sctx.shift_amount;  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1944       dk = rtmp[k];  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1960           for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j];  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1976       for (j = 0; j < nz; j++) rs += PetscAbsScalar(rtmp[bcol[j]]);  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
1991           ba[j] = rtmp[col];  in MatCholeskyFactorNumeric_SeqAIJ_inplace()
2002   PetscCall(PetscFree3(rtmp, il, jl));  in MatCholeskyFactorNumeric_SeqAIJ_inplace()