mpi/superlu_dist/superlu_dist.c

100   superlu_dist_options_t options;  member
198     if (lu->options.SolveInitialized) {  in MatDestroy_SuperLU_DIST()
200 …cStackCallExternalVoid("SuperLU_DIST:SolveFinalize", sSolveFinalize(&lu->options, &lu->sSOLVEstruc…  in MatDestroy_SuperLU_DIST()
203 …PetscStackCallExternalVoid("SuperLU_DIST:SolveFinalize", SolveFinalize(&lu->options, &lu->SOLVEstr…  in MatDestroy_SuperLU_DIST()
285 …PetscCheck(lu->options.Fact == FACTORED, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "SuperLU_DIST optio…  in MatSolve_SuperLU_DIST()
290   if (lu->options.SolveInitialized && !lu->matsolve_iscalled) {  in MatSolve_SuperLU_DIST()
293 …cStackCallExternalVoid("SuperLU_DIST:SolveFinalize", sSolveFinalize(&lu->options, &lu->sSOLVEstruc…  in MatSolve_SuperLU_DIST()
296 …PetscStackCallExternalVoid("SuperLU_DIST:SolveFinalize", SolveFinalize(&lu->options, &lu->SOLVEstr…  in MatSolve_SuperLU_DIST()
297     lu->options.SolveInitialized = NO;  in MatSolve_SuperLU_DIST()
314 …sion) PetscStackCallExternalVoid("SuperLU_DIST:pgssvx3d", psgssvx3d(&lu->options, &lu->A_sup, &lu-…  in MatSolve_SuperLU_DIST()
317 …PetscStackCallExternalVoid("SuperLU_DIST:pgssvx3d", pgssvx3d(&lu->options, &lu->A_sup, &lu->ScaleP…  in MatSolve_SuperLU_DIST()
323 …if (lu->singleprecision) PetscStackCallExternalVoid("SuperLU_DIST:pgssvx", psgssvx(&lu->options, &…  in MatSolve_SuperLU_DIST()
326 …PetscStackCallExternalVoid("SuperLU_DIST:pgssvx", pgssvx(&lu->options, &lu->A_sup, &lu->ScalePerms…  in MatSolve_SuperLU_DIST()
329 …if (lu->options.PrintStat) PetscStackCallExternalVoid("SuperLU_DIST:PStatPrint", PStatPrint(&lu->o…  in MatSolve_SuperLU_DIST()
358 …PetscCheck(lu->options.Fact == FACTORED, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "SuperLU_DIST optio…  in MatMatSolve_SuperLU_DIST()
366   if (lu->options.SolveInitialized && !lu->matmatsolve_iscalled) {  in MatMatSolve_SuperLU_DIST()
371 …PetscStackCallExternalVoid("SuperLU_DIST:SolveFinalize", SolveFinalize(&lu->options, &lu->SOLVEstr…  in MatMatSolve_SuperLU_DIST()
372     lu->options.SolveInitialized = NO;  in MatMatSolve_SuperLU_DIST()
382 …if (lu->use3d) PetscStackCallExternalVoid("SuperLU_DIST:pgssvx3d", pgssvx3d(&lu->options, &lu->A_s…  in MatMatSolve_SuperLU_DIST()
385 …PetscStackCallExternalVoid("SuperLU_DIST:pgssvx", pgssvx(&lu->options, &lu->A_sup, &lu->ScalePerms…  in MatMatSolve_SuperLU_DIST()
389 …if (lu->options.PrintStat) PetscStackCallExternalVoid("SuperLU_DIST:PStatPrint", PStatPrint(&lu->o…  in MatMatSolve_SuperLU_DIST()
416 …PetscCheck(lu->options.RowPerm == NOROWPERM, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Must set NOROW…  in MatGetInertia_SuperLU_DIST()
473   if (lu->options.Fact == DOFACT) { /* first numeric factorization */  in MatLUFactorNumeric_SuperLU_DIST()
481       lu->options.Fact = SamePattern_SameRowPerm; /* matrix has similar numerical values */  in MatLUFactorNumeric_SuperLU_DIST()
488 …PetscStackCallExternalVoid("SuperLU_DIST:SolveFinalize", sSolveFinalize(&lu->options, &lu->sSOLVEs…  in MatLUFactorNumeric_SuperLU_DIST()
493 …PetscStackCallExternalVoid("SuperLU_DIST:SolveFinalize", SolveFinalize(&lu->options, &lu->SOLVEstr…  in MatLUFactorNumeric_SuperLU_DIST()
503       lu->options.Fact = SamePattern;  in MatLUFactorNumeric_SuperLU_DIST()
511       lu->options.Fact = DOFACT;  in MatLUFactorNumeric_SuperLU_DIST()
535   if (lu->options.Fact == DOFACT) {  in MatLUFactorNumeric_SuperLU_DIST()
549 …sion) PetscStackCallExternalVoid("SuperLU_DIST:pgssvx3d", psgssvx3d(&lu->options, &lu->A_sup, &lu-…  in MatLUFactorNumeric_SuperLU_DIST()
552 …PetscStackCallExternalVoid("SuperLU_DIST:pgssvx3d", pgssvx3d(&lu->options, &lu->A_sup, &lu->ScaleP…  in MatLUFactorNumeric_SuperLU_DIST()
557 …if (lu->singleprecision) PetscStackCallExternalVoid("SuperLU_DIST:pgssvx", psgssvx(&lu->options, &…  in MatLUFactorNumeric_SuperLU_DIST()
560 …PetscStackCallExternalVoid("SuperLU_DIST:pgssvx", pgssvx(&lu->options, &lu->A_sup, &lu->ScalePerms…  in MatLUFactorNumeric_SuperLU_DIST()
577 …if (lu->options.PrintStat) PetscStackCallExternalVoid("SuperLU_DIST:PStatPrint", PStatPrint(&lu->o…  in MatLUFactorNumeric_SuperLU_DIST()
581 …lu->options.Fact = FACTORED; /* The factored form of A is supplied. Local option used by this func…  in MatLUFactorNumeric_SuperLU_DIST()
604 …PetscCall(PetscOptionsBool("-mat_superlu_dist_equil", "Equilibrate matrix", "None", lu->options.Eq…  in MatLUFactorSymbolic_SuperLU_DIST()
605   if (set && !flg) lu->options.Equil = NO;  in MatLUFactorSymbolic_SuperLU_DIST()
611       lu->options.RowPerm = NOROWPERM;  in MatLUFactorSymbolic_SuperLU_DIST()
614       lu->options.RowPerm = LargeDiag_MC64;  in MatLUFactorSymbolic_SuperLU_DIST()
617       lu->options.RowPerm = LargeDiag_AWPM;  in MatLUFactorSymbolic_SuperLU_DIST()
620       lu->options.RowPerm = MY_PERMR;  in MatLUFactorSymbolic_SuperLU_DIST()
631       lu->options.ColPerm = NATURAL;  in MatLUFactorSymbolic_SuperLU_DIST()
634       lu->options.ColPerm = MMD_AT_PLUS_A;  in MatLUFactorSymbolic_SuperLU_DIST()
637       lu->options.ColPerm = MMD_ATA;  in MatLUFactorSymbolic_SuperLU_DIST()
640       lu->options.ColPerm = METIS_AT_PLUS_A;  in MatLUFactorSymbolic_SuperLU_DIST()
643       lu->options.ColPerm = PARMETIS; /* only works for np>1 */  in MatLUFactorSymbolic_SuperLU_DIST()
650   lu->options.ReplaceTinyPivot = NO;  in MatLUFactorSymbolic_SuperLU_DIST()
651 …"-mat_superlu_dist_replacetinypivot", "Replace tiny pivots", "None", lu->options.ReplaceTinyPivot …  in MatLUFactorSymbolic_SuperLU_DIST()
652   if (set && flg) lu->options.ReplaceTinyPivot = YES;  in MatLUFactorSymbolic_SuperLU_DIST()
654   lu->options.ParSymbFact = NO;  in MatLUFactorSymbolic_SuperLU_DIST()
658     lu->options.ParSymbFact = YES;  in MatLUFactorSymbolic_SuperLU_DIST()
659 …lu->options.ColPerm     = PARMETIS; /* in v2.2, PARMETIS is forced for ParSymbFact regardless of u…  in MatLUFactorSymbolic_SuperLU_DIST()
681   lu->options.IterRefine = NOREFINE;  in MatLUFactorSymbolic_SuperLU_DIST()
682 …("-mat_superlu_dist_iterrefine", "Use iterative refinement", "None", lu->options.IterRefine == NOR…  in MatLUFactorSymbolic_SuperLU_DIST()
683   if (set && flg) lu->options.IterRefine = SLU_DOUBLE;  in MatLUFactorSymbolic_SuperLU_DIST()
685   if (PetscLogPrintInfo) lu->options.PrintStat = YES;  in MatLUFactorSymbolic_SuperLU_DIST()
686   else lu->options.PrintStat = NO;  in MatLUFactorSymbolic_SuperLU_DIST()
688 …nt factorization information", "None", (PetscBool)lu->options.PrintStat, (PetscBool *)&lu->options…  in MatLUFactorSymbolic_SuperLU_DIST()
691   lu->options.superlu_acc_offload = 1;  in MatLUFactorSymbolic_SuperLU_DIST()
692 …orization onto the GPUs", "None", (PetscBool)lu->options.superlu_acc_offload, (PetscBool *)&lu->op…  in MatLUFactorSymbolic_SuperLU_DIST()
812   superlu_dist_options_t options;  in MatView_Info_SuperLU_DIST()  local
818   options = lu->options;  in MatView_Info_SuperLU_DIST()
827 …PetscCall(PetscViewerASCIIPrintf(viewer, "  Equilibrate matrix %s \n", PetscBools[options.Equil !=…  in MatView_Info_SuperLU_DIST()
828 …PetscCall(PetscViewerASCIIPrintf(viewer, "  Replace tiny pivots %s \n", PetscBools[options.Replace…  in MatView_Info_SuperLU_DIST()
829 …PetscCall(PetscViewerASCIIPrintf(viewer, "  Use iterative refinement %s \n", PetscBools[options.It…  in MatView_Info_SuperLU_DIST()
832   switch (options.RowPerm) {  in MatView_Info_SuperLU_DIST()
849   switch (options.ColPerm) {  in MatView_Info_SuperLU_DIST()
870 …SCIIPrintf(viewer, "  Parallel symbolic factorization %s \n", PetscBools[options.ParSymbFact != NO…  in MatView_Info_SuperLU_DIST()
905   superlu_dist_options_t options;  in MatGetFactor_aij_superlu_dist()  local
934   set_default_options_dist(&options);  in MatGetFactor_aij_superlu_dist()
943     options.SymPattern             = YES;  in MatGetFactor_aij_superlu_dist()
954   lu->options              = options;  in MatGetFactor_aij_superlu_dist()
955   lu->options.Fact         = DOFACT;  in MatGetFactor_aij_superlu_dist()