Lines Matching refs:rtmp
567 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
574 PetscCall(PetscCalloc1(bs2 * mbs, &rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
635 rtmp_ptr = rtmp + vj * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
640 PetscCall(PetscArraycpy(dk, rtmp + k * bs2, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
668 rtmp_ptr = rtmp + bj[j] * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
698 rtmp_ptr = rtmp + vj * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
713 PetscCall(PetscFree(rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
740 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
746 PetscCall(PetscCalloc1(bs2 * mbs, &rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
767 rtmp_ptr = rtmp + vj * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
772 PetscCall(PetscArraycpy(dk, rtmp + k * bs2, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
800 rtmp_ptr = rtmp + bj[j] * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
830 rtmp_ptr = rtmp + vj * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
845 PetscCall(PetscFree(rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
873 MatScalar *u, *diag, *rtmp, *rtmp_ptr, dk[4], uik[4]; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
889 PetscCall(PetscCalloc1(4 * mbs, &rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
943 rtmp_ptr = rtmp + vj * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
948 PetscCall(PetscArraycpy(dk, rtmp + k * 4, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
982 rtmp_ptr = rtmp + bj[j] * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1014 rtmp_ptr = rtmp + vj * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1029 PetscCall(PetscFree(rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1052 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering() local
1068 PetscCall(PetscCalloc1(4 * mbs, &rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1085 rtmp_ptr = rtmp + vj * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1090 PetscCall(PetscArraycpy(dk, rtmp + k * 4, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1124 rtmp_ptr = rtmp + bj[j] * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1156 rtmp_ptr = rtmp + vj * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1171 PetscCall(PetscFree(rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1196 MatScalar *rtmp, *ba = b->a, *bval, *aa, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1228 PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1234 rtmp[i] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1245 rtmp[col] = aa[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1250 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1253 dk = rtmp[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1269 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1289 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1307 ba[j] = rtmp[col]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1308 rtmp[col] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1318 PetscCall(PetscFree3(rtmp, il, jl)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1353 MatScalar *rtmp, *ba = b->a, *bval, *aa = a->a, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1361 PetscCall(PetscMalloc3(mbs, &rtmp, mbs, &il, mbs, &c2r)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1369 PetscCall(PetscArrayzero(rtmp, mbs)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1374 rtmp[i] += -PetscRealPart(d); /* diagonal entry */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1379 rtmp[i] += PetscAbsScalar(v[j]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1380 rtmp[ajtmp[j]] += PetscAbsScalar(v[j]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1382 if (PetscRealPart(rtmp[i]) > sctx.shift_top) sctx.shift_top = PetscRealPart(rtmp[i]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1404 for (j = 0; j < nz; j++) rtmp[bjtmp[j]] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1412 rtmp[col] = aa[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1416 rtmp[k] += sctx.shift_amount; /* shift the diagonal of the matrix */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1419 dk = rtmp[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1435 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1452 ba[j] = rtmp[col]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1472 PetscCall(PetscFree3(rtmp, il, c2r)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1505 MatScalar *rtmp, *ba = b->a, *aa = a->a, dk, uikdi, *aval, *bval; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1522 PetscCall(PetscMalloc1(mbs, &rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1529 rtmp[i] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1540 rtmp[*acol++] = *aval++; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1545 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1548 dk = rtmp[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1566 while (nz--) rtmp[*bcol++] += uikdi * (*bval++); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1585 rs += PetscAbsScalar(rtmp[*bcol]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1604 *bval++ = rtmp[*bcol]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1605 rtmp[*bcol++] = 0.0; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1615 PetscCall(PetscFree(rtmp)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()