Lines Matching refs:rtmp
15 MatScalar *rtmp, *pc, *mwork, *pv; in MatLUFactorNumeric_SeqBAIJ_2() local
25 PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork)); in MatLUFactorNumeric_SeqBAIJ_2()
26 PetscCall(PetscArrayzero(rtmp, bs2 * n)); in MatLUFactorNumeric_SeqBAIJ_2()
34 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
39 for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
45 …for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, … in MatLUFactorNumeric_SeqBAIJ_2()
54 pc = rtmp + bs2 * row; in MatLUFactorNumeric_SeqBAIJ_2()
72 v = rtmp + 4 * pj[j]; in MatLUFactorNumeric_SeqBAIJ_2()
85 … for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
90 PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
102 … for (PetscInt j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
105 PetscCall(PetscFree2(rtmp, mwork)); in MatLUFactorNumeric_SeqBAIJ_2()
124 MatScalar *rtmp, *pc, *mwork, *pv; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
134 PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
135 PetscCall(PetscArrayzero(rtmp, bs2 * n)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
142 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
147 for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
153 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
160 pc = rtmp + bs2 * row; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
178 v = rtmp + 4 * pj[j]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
191 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
196 PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
209 for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
211 PetscCall(PetscFree2(rtmp, mwork)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
233 MatScalar *pv, *v, *rtmp, m1, m2, m3, m4, *pc, *w, *x, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
247 PetscCall(PetscMalloc1(4 * (n + 1), &rtmp)); in MatILUFactorNumeric_SeqBAIJ_2_inplace()
253 x = rtmp + 4 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
262 x = rtmp + 4 * ic[ajtmpold[j]]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
271 pc = rtmp + 4 * row; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
294 x = rtmp + 4 * pj[j]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
310 x = rtmp + 4 * pj[j]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
323 PetscCall(PetscFree(rtmp)); in MatILUFactorNumeric_SeqBAIJ_2_inplace()
344 MatScalar *pv, *v, *rtmp, *pc, *w, *x; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() local
356 PetscCall(PetscMalloc1(4 * (n + 1), &rtmp)); in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
361 x = rtmp + 4 * ajtmp[j]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
369 x = rtmp + 4 * ajtmpold[j]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
378 pc = rtmp + 4 * row; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
401 x = rtmp + 4 * pj[j]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
417 x = rtmp + 4 * pj[j]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
436 PetscCall(PetscFree(rtmp)); in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
458 MatScalar *rtmp, *pc, multiplier, *pv; in MatLUFactorNumeric_SeqBAIJ_1() local
490 PetscCall(PetscMalloc1(n + 1, &rtmp)); in MatLUFactorNumeric_SeqBAIJ_1()
500 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqBAIJ_1()
505 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatLUFactorNumeric_SeqBAIJ_1()
511 for (j = 0; j < nz; j++) rtmp[ics[ajtmp[j]]] = v[j]; in MatLUFactorNumeric_SeqBAIJ_1()
514 rtmp[i] += sctx.shift_amount; /* shift the diagonal of the matrix */ in MatLUFactorNumeric_SeqBAIJ_1()
521 pc = rtmp + row; in MatLUFactorNumeric_SeqBAIJ_1()
530 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqBAIJ_1()
543 pv[j] = rtmp[pj[j]]; in MatLUFactorNumeric_SeqBAIJ_1()
552 pv[j] = rtmp[pj[j]]; in MatLUFactorNumeric_SeqBAIJ_1()
557 sctx.pv = rtmp[i]; in MatLUFactorNumeric_SeqBAIJ_1()
560 rtmp[i] = sctx.pv; /* sctx.pv might be updated in the case of MAT_SHIFT_INBLOCKS */ in MatLUFactorNumeric_SeqBAIJ_1()
564 *pv = (PetscScalar)1.0 / rtmp[i]; in MatLUFactorNumeric_SeqBAIJ_1()
582 PetscCall(PetscFree(rtmp)); in MatLUFactorNumeric_SeqBAIJ_1()
622 MatScalar *pv, *v, *rtmp, multiplier, *pc; in MatILUFactorNumeric_SeqBAIJ_1_inplace() local
633 PetscCall(PetscMalloc1(n + 1, &rtmp)); in MatILUFactorNumeric_SeqBAIJ_1_inplace()
638 for (j = 0; j < nz; j++) rtmp[ajtmp[j]] = 0.0; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
644 for (j = 0; j < nz; j++) rtmp[ic[ajtmpold[j]]] = v[j]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
648 pc = rtmp + row; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
655 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
664 for (j = 0; j < nz; j++) pv[j] = rtmp[pj[j]]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
671 PetscCall(PetscFree(rtmp)); in MatILUFactorNumeric_SeqBAIJ_1_inplace()
756 MatScalar *rtmp, *ba = b->a, *bval, *aa = a->a, dk, uikdi; in MatCholeskyFactorNumeric_SeqBAIJ_N() local
776 PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
783 rtmp[i] = 0.0; in MatCholeskyFactorNumeric_SeqBAIJ_N()
796 rtmp[col] = aa[j]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
802 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqBAIJ_N()
805 dk = rtmp[k]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
821 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
839 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqBAIJ_N()
857 ba[j] = rtmp[col]; in MatCholeskyFactorNumeric_SeqBAIJ_N()
858 rtmp[col] = 0.0; in MatCholeskyFactorNumeric_SeqBAIJ_N()
868 PetscCall(PetscFree3(rtmp, il, jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N()
893 MatScalar *rtmp, *ba = b->a, *aa = a->a, dk, uikdi, *aval, *bval; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering() local
901 PetscCall(PetscMalloc3(am, &rtmp, am, &il, am, &jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
906 rtmp[i] = 0.0; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
922 rtmp[*acol++] = *aval++; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
928 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
931 dk = rtmp[k]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
948 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
966 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
985 *bval++ = rtmp[*bcol]; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
986 rtmp[*bcol++] = 0.0; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
996 PetscCall(PetscFree3(rtmp, il, jl)); in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()