1 #include <../src/ksp/pc/impls/factor/factor.h> /*I "petscpc.h"  I*/
3 PetscErrorCode PCFactorSetUpMatSolverType_Factor(PC pc)  in PCFactorSetUpMatSolverType_Factor()  argument
5   PC_Factor *icc = (PC_Factor *)pc->data;  in PCFactorSetUpMatSolverType_Factor()
8 …PetscCheck(pc->pmat, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONGSTATE, "You can only cal…  in PCFactorSetUpMatSolverType_Factor()
9   if (!icc->fact) PetscCall(MatGetFactor(pc->pmat, icc->solvertype, icc->factortype, &icc->fact));  in PCFactorSetUpMatSolverType_Factor()
10 …PetscCheck(icc->fact, PetscObjectComm((PetscObject)pc->pmat), PETSC_ERR_SUP, "MatFactor type %s no…  in PCFactorSetUpMatSolverType_Factor()
11              MatFactorTypes[icc->factortype], ((PetscObject)pc->pmat)->type_name, icc->solvertype);  in PCFactorSetUpMatSolverType_Factor()
15 PetscErrorCode PCFactorSetZeroPivot_Factor(PC pc, PetscReal z)  in PCFactorSetZeroPivot_Factor()  argument
17   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorSetZeroPivot_Factor()
24 PetscErrorCode PCFactorSetShiftType_Factor(PC pc, MatFactorShiftType shifttype)  in PCFactorSetShiftType_Factor()  argument
26   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetShiftType_Factor()
37 PetscErrorCode PCFactorSetShiftAmount_Factor(PC pc, PetscReal shiftamount)  in PCFactorSetShiftAmount_Factor()  argument
39   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetShiftAmount_Factor()
47 PetscErrorCode PCFactorSetDropTolerance_Factor(PC pc, PetscReal dt, PetscReal dtcol, PetscInt dtcou…  in PCFactorSetDropTolerance_Factor()  argument
49   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorSetDropTolerance_Factor()
52 …pc->setupcalled && (!ilu->info.usedt || ilu->info.dt != dt || ilu->info.dtcol != dtcol || ilu->inf…  in PCFactorSetDropTolerance_Factor()
61 PetscErrorCode PCFactorSetFill_Factor(PC pc, PetscReal fill)  in PCFactorSetFill_Factor()  argument
63   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetFill_Factor()
70 PetscErrorCode PCFactorSetMatOrderingType_Factor(PC pc, MatOrderingType ordering)  in PCFactorSetMatOrderingType_Factor()  argument
72   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetMatOrderingType_Factor()
76   if (!pc->setupcalled) {  in PCFactorSetMatOrderingType_Factor()
81 …PetscCheck(flg, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONGSTATE, "Cannot change orderin…  in PCFactorSetMatOrderingType_Factor()
86 PetscErrorCode PCFactorGetLevels_Factor(PC pc, PetscInt *levels)  in PCFactorGetLevels_Factor()  argument
88   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorGetLevels_Factor()
95 PetscErrorCode PCFactorGetZeroPivot_Factor(PC pc, PetscReal *pivot)  in PCFactorGetZeroPivot_Factor()  argument
97   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorGetZeroPivot_Factor()
104 PetscErrorCode PCFactorGetShiftAmount_Factor(PC pc, PetscReal *shift)  in PCFactorGetShiftAmount_Factor()  argument
106   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorGetShiftAmount_Factor()
113 PetscErrorCode PCFactorGetShiftType_Factor(PC pc, MatFactorShiftType *type)  in PCFactorGetShiftType_Factor()  argument
115   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorGetShiftType_Factor()
122 PetscErrorCode PCFactorSetLevels_Factor(PC pc, PetscInt levels)  in PCFactorSetLevels_Factor()  argument
124   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorSetLevels_Factor()
127   if (!pc->setupcalled) ilu->info.levels = levels;  in PCFactorSetLevels_Factor()
129     PetscUseTypeMethod(pc, reset);  /* remove previous factored matrices */  in PCFactorSetLevels_Factor()
130 … pc->setupcalled  = PETSC_FALSE; /* force a complete rebuild of preconditioner factored matrices */  in PCFactorSetLevels_Factor()
132 …} else PetscCheck(!ilu->info.usedt, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONGSTATE, "C…  in PCFactorSetLevels_Factor()
136 PetscErrorCode PCFactorSetAllowDiagonalFill_Factor(PC pc, PetscBool flg)  in PCFactorSetAllowDiagonalFill_Factor()  argument
138   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetAllowDiagonalFill_Factor()
145 PetscErrorCode PCFactorGetAllowDiagonalFill_Factor(PC pc, PetscBool *flg)  in PCFactorGetAllowDiagonalFill_Factor()  argument
147   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorGetAllowDiagonalFill_Factor()
154 PetscErrorCode PCFactorSetPivotInBlocks_Factor(PC pc, PetscBool pivot)  in PCFactorSetPivotInBlocks_Factor()  argument
156   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetPivotInBlocks_Factor()
163 PetscErrorCode PCFactorGetMatrix_Factor(PC pc, Mat *mat)  in PCFactorGetMatrix_Factor()  argument
165   PC_Factor *ilu = (PC_Factor *)pc->data;  in PCFactorGetMatrix_Factor()
168 …PetscCheck(ilu->fact, PetscObjectComm((PetscObject)pc), PETSC_ERR_ORDER, "Matrix not yet factored;…  in PCFactorGetMatrix_Factor()
176 PetscErrorCode PCFactorSetMatSolverType_Factor(PC pc, MatSolverType stype)  in PCFactorSetMatSolverType_Factor()  argument
178   PC_Factor *lu = (PC_Factor *)pc->data;  in PCFactorSetMatSolverType_Factor()
187 …, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_WRONGSTATE, "Cannot change solver matrix package…  in PCFactorSetMatSolverType_Factor()
195 PetscErrorCode PCFactorGetMatSolverType_Factor(PC pc, MatSolverType *stype)  in PCFactorGetMatSolverType_Factor()  argument
197   PC_Factor *lu = (PC_Factor *)pc->data;  in PCFactorGetMatSolverType_Factor()
204 PetscErrorCode PCFactorSetColumnPivot_Factor(PC pc, PetscReal dtcol)  in PCFactorSetColumnPivot_Factor()  argument
206   PC_Factor *dir = (PC_Factor *)pc->data;  in PCFactorSetColumnPivot_Factor()
209 …PetscCheck(dtcol >= 0.0 && dtcol <= 1.0, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_OUTOFRANG…  in PCFactorSetColumnPivot_Factor()
214 PetscErrorCode PCSetFromOptions_Factor(PC pc, PetscOptionItems PetscOptionsObject)  in PCSetFromOptions_Factor()  argument
216   PC_Factor        *factor = (PC_Factor *)pc->data;  in PCSetFromOptions_Factor()
224   PetscCall(PCFactorGetUseInPlace(pc, &inplace));  in PCSetFromOptions_Factor()
226   if (set) PetscCall(PCFactorSetUseInPlace(pc, flg));  in PCSetFromOptions_Factor()
230   if (flg) PetscCall(PCFactorSetShiftType(pc, (MatFactorShiftType)etmp));  in PCSetFromOptions_Factor()
237   if (set) PetscCall(PCFactorSetPivotInBlocks(pc, flg));  in PCSetFromOptions_Factor()
240   if (set) PetscCall(PCFactorSetReuseFill(pc, flg));  in PCSetFromOptions_Factor()
242   if (set) PetscCall(PCFactorSetReuseOrdering(pc, flg));  in PCSetFromOptions_Factor()
246   if (flg) PetscCall(PCFactorSetMatSolverType(pc, solvertype));  in PCSetFromOptions_Factor()
247   PetscCall(PCFactorSetDefaultOrdering_Factor(pc));  in PCSetFromOptions_Factor()
250   if (flg) PetscCall(PCFactorSetMatOrderingType(pc, tname));  in PCSetFromOptions_Factor()
256 PetscErrorCode PCView_Factor(PC pc, PetscViewer viewer)  in PCView_Factor()  argument
258   PC_Factor        *factor = (PC_Factor *)pc->data;  in PCView_Factor()