| /libCEED/examples/petsc/ |
| H A D | bpssphere.c | 58 KSP ksp; in main() local
202 PetscCall(KSPCreate(comm, &ksp)); in main()
205 PetscCall(KSPGetPC(ksp, &pc)); in main()
217 PetscCall(KSPSetType(ksp, KSPCG)); in main()
218 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in main()
219 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in main()
221 PetscCall(KSPSetFromOptions(ksp)); in main()
222 PetscCall(KSPSetOperators(ksp, mat_O, mat_O)); in main()
226 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, 1)); in main()
228 PetscCall(KSPSolve(ksp, rhs, X)); in main()
[all …]
|
| H A D | bpsswarm.c | 61 KSP ksp; in main() local
299 PetscCall(KSPCreate(comm, &ksp)); in main()
302 PetscCall(KSPGetPC(ksp, &pc)); in main()
309 PetscCall(KSPSetType(ksp, KSPCG)); in main()
310 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in main()
311 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in main()
313 PetscCall(KSPSetFromOptions(ksp)); in main()
314 PetscCall(KSPSetOperators(ksp, mat_O, mat_O)); in main()
318 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, 1)); in main()
320 PetscCall(KSPSolve(ksp, rhs, X)); in main()
[all …]
|
| H A D | bps.c | 60 KSP ksp; in RunWithDM() local
192 PetscCall(KSPCreate(rp->comm, &ksp)); in RunWithDM()
195 PetscCall(KSPGetPC(ksp, &pc)); in RunWithDM()
207 PetscCall(KSPSetType(ksp, KSPCG)); in RunWithDM()
208 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in RunWithDM()
209 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in RunWithDM()
211 PetscCall(KSPSetOperators(ksp, mat_O, mat_O)); in RunWithDM()
214 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, 1)); in RunWithDM()
216 PetscCall(KSPSolve(ksp, rhs, X)); in RunWithDM()
221 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, rp->ksp_max_it_clip[0])); in RunWithDM()
[all …]
|
| H A D | multigrid.c | 59 KSP ksp; in main() local
379 // Set up KSP in main()
380 PetscCall(KSPCreate(comm, &ksp)); in main()
382 PetscCall(KSPSetType(ksp, KSPCG)); in main()
383 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in main()
384 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in main()
386 PetscCall(KSPSetFromOptions(ksp)); in main()
387 PetscCall(KSPSetOperators(ksp, mat_O[fine_level], mat_O[fine_level])); in main()
390 PetscCall(KSPGetPC(ksp, &pc)); in main()
399 KSP smoother; in main()
[all …]
|
| H A D | bpsraw.c | 350 KSP ksp; in main() local
688 PetscCall(KSPCreate(comm, &ksp)); in main()
691 PetscCall(KSPGetPC(ksp, &pc)); in main()
698 PetscCall(KSPSetType(ksp, KSPCG)); in main()
699 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in main()
700 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in main()
702 PetscCall(KSPSetOperators(ksp, mat, mat)); in main()
704 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, 1)); in main()
706 PetscCall(KSPSolve(ksp, rhs, X)); in main()
711 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, ksp_max_it_clip[0])); in main()
[all …]
|
| /libCEED/examples/fluids/src/ |
| H A D | grid_anisotropy_tensor.c | 23 KSP ksp; in GridAnisotropyTensorProjectionSetupApply() local
75 { // -- Setup KSP for L^2 projection in GridAnisotropyTensorProjectionSetupApply()
80 PetscCall(KSPCreate(comm, &ksp)); in GridAnisotropyTensorProjectionSetupApply()
81 PetscCall(KSPSetOptionsPrefix(ksp, "grid_anisotropy_tensor_projection_")); in GridAnisotropyTensorProjectionSetupApply()
84 PetscCall(KSPGetPC(ksp, &pc)); in GridAnisotropyTensorProjectionSetupApply()
87 PetscCall(KSPSetType(ksp, KSPCG)); in GridAnisotropyTensorProjectionSetupApply()
88 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in GridAnisotropyTensorProjectionSetupApply()
89 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in GridAnisotropyTensorProjectionSetupApply()
91 PetscCall(KSPSetFromOptions_WithMatCeed(ksp, mat_mass)); in GridAnisotropyTensorProjectionSetupApply()
104 PetscCall(KSPSolve(ksp, Grid_Anisotropy, Grid_Anisotropy)); in GridAnisotropyTensorProjectionSetupApply()
[all …]
|
| H A D | petsc_ops.c | 302 * @brief Return Mats for KSP solve
304 …* Uses command-line flag with `ksp`'s prefix to determine if mat_ceed should be used directly or w…
307 …d`, then `Pmat` is either assembled or uses `mat_ceed` based on the preconditioner choice in `ksp`.
309 * @param[in] ksp `KSP` object for used for solving
315 PetscErrorCode CreateSolveOperatorsFromMatCeed(KSP ksp, Mat mat_ceed, PetscBool assemble, Mat *Amat… in CreateSolveOperatorsFromMatCeed() argument
322 PetscCall(KSPGetOptionsPrefix(ksp, &ksp_prefix)); in CreateSolveOperatorsFromMatCeed()
324 …PetscCall(PetscOptionsBool("-matceed_assemble_amat", "Assemble the A matrix for KSP solve", NULL, … in CreateSolveOperatorsFromMatCeed()
344 PetscCall(KSPGetPC(ksp, &pc)); in CreateSolveOperatorsFromMatCeed()
364 * @param[in] ksp `KSP` of the solve
367 PetscErrorCode KSPSetFromOptions_WithMatCeed(KSP ksp, Mat mat_ceed) { in KSPSetFromOptions_WithMatCeed() argument
[all …]
|
| H A D | velocity_gradient_projection.c | 106 { // -- Setup KSP for L^2 projection with lumped mass operator in VelocityGradientProjectionSetup()
112 PetscCall(KSPCreate(comm, &grad_velo_proj->ksp)); in VelocityGradientProjectionSetup()
113 PetscCall(KSPSetOptionsPrefix(grad_velo_proj->ksp, "velocity_gradient_projection_")); in VelocityGradientProjectionSetup()
116 PetscCall(KSPGetPC(grad_velo_proj->ksp, &pc)); in VelocityGradientProjectionSetup()
119 PetscCall(KSPSetType(grad_velo_proj->ksp, KSPPREONLY)); in VelocityGradientProjectionSetup()
121 PetscCall(KSPSetFromOptions_WithMatCeed(grad_velo_proj->ksp, mat_mass)); in VelocityGradientProjectionSetup()
143 PetscCall(KSPSolve(grad_velo_proj->ksp, VelocityGradient, VelocityGradient)); in VelocityGradientProjectionApply()
|
| H A D | differential_filter.c | 85 { // Setup LHS Operator and KSP for the differential filtering solve in DifferentialFilterCreateOperators()
164 PetscCall(KSPCreate(PetscObjectComm((PetscObject)dm_filter), &diff_filter->ksp)); in DifferentialFilterCreateOperators()
165 PetscCall(KSPSetOptionsPrefix(diff_filter->ksp, "diff_filter_")); in DifferentialFilterCreateOperators()
168 PetscCall(KSPGetPC(diff_filter->ksp, &pc)); in DifferentialFilterCreateOperators()
171 PetscCall(KSPSetType(diff_filter->ksp, KSPCG)); in DifferentialFilterCreateOperators()
172 PetscCall(KSPSetNormType(diff_filter->ksp, KSP_NORM_NATURAL)); in DifferentialFilterCreateOperators()
173 … PetscCall(KSPSetTolerances(diff_filter->ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in DifferentialFilterCreateOperators()
175 PetscCall(KSPSetFromOptions_WithMatCeed(diff_filter->ksp, mat_lhs)); in DifferentialFilterCreateOperators()
306 PetscCall(KSPSolve(diff_filter->ksp, RHS, Filtered_Solution)); in DifferentialFilterApply()
336 PetscCall(KSPDestroy(&diff_filter->ksp)); in DifferentialFilterDataDestroy()
|
| H A D | turb_spanstats.c | 327 { // Setup KSP for L^2 projection in SetupL2ProjectionStats()
329 KSP ksp; in SetupL2ProjectionStats() local
333 PetscCall(KSPCreate(PetscObjectComm((PetscObject)user->spanstats.dm), &ksp)); in SetupL2ProjectionStats()
334 PetscCall(KSPSetOptionsPrefix(ksp, "turbulence_spanstats_")); in SetupL2ProjectionStats()
337 PetscCall(KSPGetPC(ksp, &pc)); in SetupL2ProjectionStats()
340 PetscCall(KSPSetType(ksp, KSPCG)); in SetupL2ProjectionStats()
341 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in SetupL2ProjectionStats()
342 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in SetupL2ProjectionStats()
344 PetscCall(KSPSetFromOptions_WithMatCeed(ksp, mat_mass)); in SetupL2ProjectionStats()
345 user->spanstats.ksp = ksp; in SetupL2ProjectionStats()
[all …]
|
| H A D | setupts.c | 332 KSP ksp; in TSSolve_NS() local
336 PetscCall(SNESGetKSP(snes, &ksp)); in TSSolve_NS()
337 … PetscCall(CreateSolveOperatorsFromMatCeed(ksp, user->mat_ijacobian, PETSC_FALSE, &Amat, &Pmat)); in TSSolve_NS()
|
| H A D | setuplibceed.c | 16 // @brief Create CeedOperator for unstabilized mass KSP for explicit timestepping
56 // @brief Create KSP to solve the inverse mass operator for explicit time stepping schemes
66 { // -- Setup KSP for mass operator in CreateKSPMass()
|
| H A D | misc.c | 285 PetscCall(KSPDestroy(&context->ksp)); in NodalProjectionDataDestroy()
|
| /libCEED/examples/fluids/include/ |
| H A D | petsc_ops.h | 41 PetscErrorCode CreateSolveOperatorsFromMatCeed(KSP ksp, Mat mat_ceed, PetscBool assemble, Mat *Amat…
42 PetscErrorCode KSPSetFromOptions_WithMatCeed(KSP ksp, Mat mat_ceed);
|
| /libCEED/examples/petsc/src/ |
| H A D | swarmutils.c | 501 KSP ksp; in DMSwarmProjectFromSwarmToCells() local
527 // Setup KSP in DMSwarmProjectFromSwarmToCells()
531 PetscCall(KSPCreate(comm, &ksp)); in DMSwarmProjectFromSwarmToCells()
532 PetscCall(KSPGetPC(ksp, &pc)); in DMSwarmProjectFromSwarmToCells()
535 PetscCall(KSPSetType(ksp, KSPCG)); in DMSwarmProjectFromSwarmToCells()
536 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in DMSwarmProjectFromSwarmToCells()
537 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in DMSwarmProjectFromSwarmToCells()
538 PetscCall(KSPSetOperators(ksp, M, M)); in DMSwarmProjectFromSwarmToCells()
539 PetscCall(KSPSetFromOptions(ksp)); in DMSwarmProjectFromSwarmToCells()
540 PetscCall(PetscObjectSetName((PetscObject)ksp, "Swarm-to-Mesh Projection")); in DMSwarmProjectFromSwarmToCells()
[all …]
|
| /libCEED/examples/solids/ |
| H A D | elasticity.c | 471 // Setup KSP in main()
475 KSP ksp; in main() local
477 // -- KSP in main()
478 PetscCall(SNESGetKSP(snes, &ksp)); in main()
479 PetscCall(KSPSetType(ksp, KSPCG)); in main()
480 PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL)); in main()
481 PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT)); in main()
482 PetscCall(KSPSetOptionsPrefix(ksp, "outer_")); in main()
485 PetscCall(KSPGetPC(ksp, &pc)); in main()
504 KSP ksp_smoother, ksp_est; in main()
[all …]
|
| /libCEED/examples/fluids/ |
| H A D | navierstokes.h | 143 KSP ksp; // For the L^2 projection solve member
154 KSP ksp; member
163 KSP ksp; member
183 KSP mass_ksp;
|
| H A D | advection.yaml | 15 ksp:
|
| /libCEED/benchmarks/ |
| H A D | postprocess_base.py | 94 elif 'Total KSP Iterations' in line:
|
| /libCEED/examples/fluids/problems/ |
| H A D | advection.c | 18 // @brief Create CeedOperator for stabilized mass KSP for explicit timestepping
|
| H A D | newtonian.c | 155 // @brief Create CeedOperator for stabilized mass KSP for explicit timestepping
|