Lines Matching refs:sctx

1198   FactorShiftCtx  sctx;  in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()  local
1202 PetscCall(PetscMemzero(&sctx, sizeof(FactorShiftCtx))); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1231 sctx.newshift = PETSC_FALSE; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1250 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1294 sctx.rs = rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1295 sctx.pv = dk; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1296 PetscCall(MatPivotCheck(C, A, info, &sctx, k)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1297 if (sctx.newshift) break; /* sctx.shift_amount is updated */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1298 dk = sctx.pv; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1317 } while (sctx.newshift); in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1332 if (sctx.nshift) { in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1334 …"number of shiftnz tries %" PetscInt_FMT ", shift_amount %g\n", sctx.nshift, (double)sctx.shift_am… in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1336 …"number of shiftpd tries %" PetscInt_FMT ", shift_amount %g\n", sctx.nshift, (double)sctx.shift_am… in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1354 FactorShiftCtx sctx; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1364 PetscCall(PetscMemzero(&sctx, sizeof(FactorShiftCtx))); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1367 sctx.shift_top = info->zeropivot; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1382 if (PetscRealPart(rtmp[i]) > sctx.shift_top) sctx.shift_top = PetscRealPart(rtmp[i]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1384 sctx.shift_top *= 1.1; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1385 sctx.nshift_max = 5; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1386 sctx.shift_lo = 0.; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1387 sctx.shift_hi = 1.; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1395 sctx.newshift = PETSC_FALSE; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1416 rtmp[k] += sctx.shift_amount; /* shift the diagonal of the matrix */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1462 sctx.rs = rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1463 sctx.pv = dk; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1464 PetscCall(MatPivotCheck(B, A, info, &sctx, k)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1465 if (sctx.newshift) break; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1466 dk = sctx.pv; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1470 } while (sctx.newshift); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1487 if (sctx.nshift) { in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1489 … up by %e fraction top_value %e\n", sctx.nshift, (double)sctx.shift_amount, (double)sctx.shift_fra… in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1491 …number of shift_nz tries %" PetscInt_FMT ", shift_amount %g\n", sctx.nshift, (double)sctx.shift_am… in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1493 …ber of shift_inblocks applied %" PetscInt_FMT ", each shift_amount %g\n", sctx.nshift, (double)inf… in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1507 FactorShiftCtx sctx; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1511 PetscCall(PetscMemzero(&sctx, sizeof(FactorShiftCtx))); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1526 sctx.newshift = PETSC_FALSE; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1545 if (sctx.nshift) rtmp[k] += sctx.shift_amount; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1590 sctx.rs = rs; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1591 sctx.pv = dk; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1592 PetscCall(MatPivotCheck(C, A, info, &sctx, k)); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1593 if (sctx.newshift) break; /* sctx.shift_amount is updated */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1594 dk = sctx.pv; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1614 } while (sctx.newshift); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1628 if (sctx.nshift) { in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1630 …"number of shiftnz tries %" PetscInt_FMT ", shift_amount %g\n", sctx.nshift, (double)sctx.shift_am… in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1632 …"number of shiftpd tries %" PetscInt_FMT ", shift_amount %g\n", sctx.nshift, (double)sctx.shift_am… in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()