Lines Matching refs:pv
15 MatScalar *rtmp, *pc, *mwork, *pv; in MatLUFactorNumeric_SeqBAIJ_2() local
62 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_2()
64 PetscCall(PetscKernel_A_gets_A_times_B_2(pc, pv, mwork)); in MatLUFactorNumeric_SeqBAIJ_2()
67 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_2()
73 PetscCall(PetscKernel_A_gets_A_minus_B_times_C_2(v, pc, pv)); in MatLUFactorNumeric_SeqBAIJ_2()
74 pv += 4; in MatLUFactorNumeric_SeqBAIJ_2()
82 pv = b->a + bs2 * bi[i]; 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()
88 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_2()
90 PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2)); in MatLUFactorNumeric_SeqBAIJ_2()
94 PetscCall(PetscKernel_A_gets_inverse_A_2(pv, shift, allowzeropivot, &zeropivotdetected)); in MatLUFactorNumeric_SeqBAIJ_2()
99 pv = b->a + bs2 * (bdiag[i + 1] + 1); 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()
124 MatScalar *rtmp, *pc, *mwork, *pv; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
168 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
170 PetscCall(PetscKernel_A_gets_A_times_B_2(pc, pv, mwork)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
173 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
179 PetscCall(PetscKernel_A_gets_A_minus_B_times_C_2(v, pc, pv)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
180 pv += 4; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
188 pv = b->a + bs2 * bi[i]; 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()
194 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
196 PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
197 PetscCall(PetscKernel_A_gets_inverse_A_2(pv, shift, allowzeropivot, &zeropivotdetected)); in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
206 pv = b->a + bs2 * (bdiag[i + 1] + 1); 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()
233 MatScalar *pv, *v, *rtmp, m1, m2, m3, m4, *pc, *w, *x, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_2_inplace() local
277 pv = ba + 4 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
279 x1 = pv[0]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
280 x2 = pv[1]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
281 x3 = pv[2]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
282 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
288 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
290 x1 = pv[0]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
291 x2 = pv[1]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
292 x3 = pv[2]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
293 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
299 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
306 pv = ba + 4 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
311 pv[0] = x[0]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
312 pv[1] = x[1]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
313 pv[2] = x[2]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
314 pv[3] = x[3]; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
315 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_inplace()
344 MatScalar *pv, *v, *rtmp, *pc, *w, *x; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace() local
384 pv = ba + 4 * diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
386 x1 = pv[0]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
387 x2 = pv[1]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
388 x3 = pv[2]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
389 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
395 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
397 x1 = pv[0]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
398 x2 = pv[1]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
399 x3 = pv[2]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
400 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
406 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
413 pv = ba + 4 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
418 pv[0] = x[0]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
419 pv[1] = x[1]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
420 pv[2] = x[2]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
421 pv[3] = x[3]; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
428 pv += 4; in MatILUFactorNumeric_SeqBAIJ_2_NaturalOrdering_inplace()
458 MatScalar *rtmp, *pc, multiplier, *pv; in MatLUFactorNumeric_SeqBAIJ_1() local
523 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_1()
524 multiplier = *pc * (*pv); in MatLUFactorNumeric_SeqBAIJ_1()
528 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_1()
530 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatLUFactorNumeric_SeqBAIJ_1()
539 pv = b->a + bi[i]; in MatLUFactorNumeric_SeqBAIJ_1()
543 pv[j] = rtmp[pj[j]]; in MatLUFactorNumeric_SeqBAIJ_1()
544 rs += PetscAbsScalar(pv[j]); in MatLUFactorNumeric_SeqBAIJ_1()
548 pv = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_1()
552 pv[j] = rtmp[pj[j]]; in MatLUFactorNumeric_SeqBAIJ_1()
553 rs += PetscAbsScalar(pv[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()
563 pv = b->a + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_1()
564 *pv = (PetscScalar)1.0 / rtmp[i]; in MatLUFactorNumeric_SeqBAIJ_1()
622 MatScalar *pv, *v, *rtmp, multiplier, *pc; in MatILUFactorNumeric_SeqBAIJ_1_inplace() local
650 pv = ba + diag_offset[row]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
652 multiplier = *pc * *pv++; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
655 for (j = 0; j < nz; j++) rtmp[pj[j]] -= multiplier * pv[j]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
661 pv = ba + bi[i]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
664 for (j = 0; j < nz; j++) pv[j] = rtmp[pj[j]]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
667 …PetscCheck(pv[diag] != 0.0, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Zero pivot: row in original … in MatILUFactorNumeric_SeqBAIJ_1_inplace()
668 pv[diag] = 1.0 / pv[diag]; in MatILUFactorNumeric_SeqBAIJ_1_inplace()
845 sctx.pv = dk; in MatCholeskyFactorNumeric_SeqBAIJ_N()
848 dk = sctx.pv; in MatCholeskyFactorNumeric_SeqBAIJ_N()
972 sctx.pv = dk; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()
975 dk = sctx.pv; in MatCholeskyFactorNumeric_SeqBAIJ_N_NaturalOrdering()