59   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSetUp_SPAI()  local
65   if (ispai->sp) {  in PCSetUp_SPAI()
66     PetscCall(ConvertMatToMatrix(ispai->comm_spai, pc->pmat, pc->pmat, &ispai->B));  in PCSetUp_SPAI()
70     PetscCall(ConvertMatToMatrix(ispai->comm_spai, pc->pmat, AT, &ispai->B));  in PCSetUp_SPAI()
90 …nal(bspai, ispai->B, &ispai->M, stdout, ispai->epsilon, ispai->nbsteps, ispai->max, ispai->maxnew,…  in PCSetUp_SPAI()
92   PetscCall(ConvertMatrixToMat(PetscObjectComm((PetscObject)pc), ispai->M, &ispai->PM));  in PCSetUp_SPAI()
95   sp_free_matrix(ispai->B);  in PCSetUp_SPAI()
96   sp_free_matrix(ispai->M);  in PCSetUp_SPAI()
102   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCApply_SPAI()  local
106   PetscCall(MatMult(ispai->PM, xx, y));  in PCApply_SPAI()
112   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCMatApply_SPAI()  local
116   PetscCall(MatMatMult(ispai->PM, X, MAT_REUSE_MATRIX, PETSC_CURRENT, &Y));  in PCMatApply_SPAI()
122   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCDestroy_SPAI()  local
125   PetscCall(MatDestroy(&ispai->PM));  in PCDestroy_SPAI()
126   PetscCallMPI(MPI_Comm_free(&ispai->comm_spai));  in PCDestroy_SPAI()
141   PC_SPAI  *ispai = (PC_SPAI *)pc->data;  in PCView_SPAI()  local
147     PetscCall(PetscViewerASCIIPrintf(viewer, "    epsilon %g\n", ispai->epsilon));  in PCView_SPAI()
148     PetscCall(PetscViewerASCIIPrintf(viewer, "    nbsteps %d\n", ispai->nbsteps));  in PCView_SPAI()
149     PetscCall(PetscViewerASCIIPrintf(viewer, "    max %d\n", ispai->max));  in PCView_SPAI()
150     PetscCall(PetscViewerASCIIPrintf(viewer, "    maxnew %d\n", ispai->maxnew));  in PCView_SPAI()
151     PetscCall(PetscViewerASCIIPrintf(viewer, "    block_size %d\n", ispai->block_size));  in PCView_SPAI()
152     PetscCall(PetscViewerASCIIPrintf(viewer, "    cache_size %d\n", ispai->cache_size));  in PCView_SPAI()
153     PetscCall(PetscViewerASCIIPrintf(viewer, "    verbose %d\n", ispai->verbose));  in PCView_SPAI()
154     PetscCall(PetscViewerASCIIPrintf(viewer, "    sp %d\n", ispai->sp));  in PCView_SPAI()
161   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetEpsilon_SPAI()  local
164   ispai->epsilon = (double)epsilon1;  in PCSPAISetEpsilon_SPAI()
170   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetNBSteps_SPAI()  local
173   ispai->nbsteps = (int)nbsteps1;  in PCSPAISetNBSteps_SPAI()
180   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetMax_SPAI()  local
183   ispai->max = (int)max1;  in PCSPAISetMax_SPAI()
189   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetMaxNew_SPAI()  local
192   ispai->maxnew = (int)maxnew1;  in PCSPAISetMaxNew_SPAI()
198   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetBlockSize_SPAI()  local
201   ispai->block_size = (int)block_size1;  in PCSPAISetBlockSize_SPAI()
207   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetCacheSize_SPAI()  local
210   ispai->cache_size = (int)cache_size;  in PCSPAISetCacheSize_SPAI()
216   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetVerbose_SPAI()  local
219   ispai->verbose = (int)verbose;  in PCSPAISetVerbose_SPAI()
225   PC_SPAI *ispai = (PC_SPAI *)pc->data;  in PCSPAISetSp_SPAI()  local
228   ispai->sp = (int)sp;  in PCSPAISetSp_SPAI()
435   PC_SPAI  *ispai = (PC_SPAI *)pc->data;  in PCSetFromOptions_SPAI()  local
442 …PetscCall(PetscOptionsReal("-pc_spai_epsilon", "", "PCSPAISetEpsilon", ispai->epsilon, &epsilon1, …  in PCSetFromOptions_SPAI()
444 …PetscCall(PetscOptionsInt("-pc_spai_nbsteps", "", "PCSPAISetNBSteps", ispai->nbsteps, &nbsteps1, &…  in PCSetFromOptions_SPAI()
447   PetscCall(PetscOptionsInt("-pc_spai_max", "", "PCSPAISetMax", ispai->max, &max1, &flg));  in PCSetFromOptions_SPAI()
449 …PetscCall(PetscOptionsInt("-pc_spai_maxnew", "", "PCSPAISetMaxNew", ispai->maxnew, &maxnew1, &flg)…  in PCSetFromOptions_SPAI()
451 …PetscCall(PetscOptionsInt("-pc_spai_block_size", "", "PCSPAISetBlockSize", ispai->block_size, &blo…  in PCSetFromOptions_SPAI()
453 …PetscCall(PetscOptionsInt("-pc_spai_cache_size", "", "PCSPAISetCacheSize", ispai->cache_size, &cac…  in PCSetFromOptions_SPAI()
455 …PetscCall(PetscOptionsInt("-pc_spai_verbose", "", "PCSPAISetVerbose", ispai->verbose, &verbose, &f…  in PCSetFromOptions_SPAI()
457   PetscCall(PetscOptionsInt("-pc_spai_sp", "", "PCSPAISetSp", ispai->sp, &sp, &flg));  in PCSetFromOptions_SPAI()
488   PC_SPAI *ispai;  in PCCreate_SPAI()  local
491   PetscCall(PetscNew(&ispai));  in PCCreate_SPAI()
492   pc->data = ispai;  in PCCreate_SPAI()
502   ispai->epsilon    = .4;  in PCCreate_SPAI()
503   ispai->nbsteps    = 5;  in PCCreate_SPAI()
504   ispai->max        = 5000;  in PCCreate_SPAI()
505   ispai->maxnew     = 5;  in PCCreate_SPAI()
506   ispai->block_size = 1;  in PCCreate_SPAI()
507   ispai->cache_size = 5;  in PCCreate_SPAI()
508   ispai->verbose    = 0;  in PCCreate_SPAI()
510   ispai->sp = 1;  in PCCreate_SPAI()
511   PetscCallMPI(MPI_Comm_dup(PetscObjectComm((PetscObject)pc), &ispai->comm_spai));  in PCCreate_SPAI()