| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex66.c | 70 MatNullSpace nullspace; in ComputeRHS() local
92 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
93 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
94 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
104 MatNullSpace nullspace; in ComputeJacobian() local
139 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeJacobian()
140 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeJacobian()
141 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeJacobian()
|
| H A D | ex50.c | 74 MatNullSpace nullspace; in ComputeRHS() local
96 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
97 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
98 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
108 MatNullSpace nullspace; in ComputeJacobian() local
183 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeJacobian()
184 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeJacobian()
185 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeJacobian()
|
| H A D | ex32.c | 108 MatNullSpace nullspace; in ComputeRHS() local
110 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
111 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
112 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
203 MatNullSpace nullspace; in ComputeMatrix() local
205 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeMatrix()
206 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeMatrix()
207 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeMatrix()
|
| H A D | ex29.c | 133 MatNullSpace nullspace; in ComputeRHS() local
135 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
136 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
137 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
263 MatNullSpace nullspace; in ComputeMatrix() local
265 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeMatrix()
266 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeMatrix()
267 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeMatrix()
|
| H A D | ex34.c | 101 MatNullSpace nullspace; in ComputeRHS() local
127 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
128 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
129 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
139 MatNullSpace nullspace; in ComputeMatrix() local
301 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeMatrix()
302 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeMatrix()
303 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeMatrix()
|
| H A D | ex36.cxx | 333 MatNullSpace nullspace; in ComputeRHS_MOAB() local
335 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS_MOAB()
336 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS_MOAB()
337 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS_MOAB()
442 MatNullSpace nullspace; in ComputeMatrix_MOAB() local
444 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeMatrix_MOAB()
445 PetscCall(MatSetNullSpace(jac, nullspace)); in ComputeMatrix_MOAB()
446 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeMatrix_MOAB()
|
| H A D | ex35.cxx | 370 MatNullSpace nullspace; in ComputeRHS() local
371 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS()
372 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS()
373 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS()
482 MatNullSpace nullspace; in ComputeMatrix() local
484 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeMatrix()
485 PetscCall(MatSetNullSpace(J, nullspace)); in ComputeMatrix()
486 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeMatrix()
|
| H A D | ex73.c | 354 MatNullSpace nullspace = NULL; in DMCreateMatrix_ShellDA() local
357 PetscCall(MatGetNullSpace(*A, &nullspace)); in DMCreateMatrix_ShellDA()
358 if (!nullspace) { in DMCreateMatrix_ShellDA()
360 PetscCall(MatNullSpaceCreate(comm, PETSC_TRUE, 0, 0, &nullspace)); in DMCreateMatrix_ShellDA()
361 PetscCall(MatSetNullSpace(*A, nullspace)); in DMCreateMatrix_ShellDA()
362 PetscCall(MatNullSpaceDestroy(&nullspace)); in DMCreateMatrix_ShellDA()
965 MatNullSpace nullspace; in ComputeRHS_DMDA() local
967 PetscCall(MatNullSpaceCreate(PETSC_COMM_WORLD, PETSC_TRUE, 0, 0, &nullspace)); in ComputeRHS_DMDA()
968 PetscCall(MatNullSpaceRemove(nullspace, b)); in ComputeRHS_DMDA()
969 PetscCall(MatNullSpaceDestroy(&nullspace)); in ComputeRHS_DMDA()
|
| /petsc/src/ksp/pc/impls/telescope/ |
| H A D | telescope_coarsedm.c | 246 …scopeSubNullSpaceCreate_CoarseDM(PC pc, PC_Telescope sred, MatNullSpace nullspace, MatNullSpace *s… in PCTelescopeSubNullSpaceCreate_CoarseDM() argument
258 PetscCall(MatNullSpaceGetVecs(nullspace, &has_const, &n, &vecs)); in PCTelescopeSubNullSpaceCreate_CoarseDM()
285 MatNullSpace nullspace, sub_nullspace; in PCTelescopeMatNullSpaceCreate_CoarseDM() local
286 PetscCall(MatGetNullSpace(B, &nullspace)); in PCTelescopeMatNullSpaceCreate_CoarseDM()
287 if (nullspace) { in PCTelescopeMatNullSpaceCreate_CoarseDM()
289 PetscCall(PCTelescopeSubNullSpaceCreate_CoarseDM(pc, sred, nullspace, &sub_nullspace)); in PCTelescopeMatNullSpaceCreate_CoarseDM()
294 if (nullspace->remove && !nullspace->rmctx) { in PCTelescopeMatNullSpaceCreate_CoarseDM()
295 PetscCall(MatNullSpaceSetFunction(sub_nullspace, nullspace->remove, context)); in PCTelescopeMatNullSpaceCreate_CoarseDM()
296 } else if (nullspace->remove && nullspace->rmctx) { in PCTelescopeMatNullSpaceCreate_CoarseDM()
303 PetscCall(MatNullSpaceSetFunction(sub_nullspace, nullspace->remove, context)); in PCTelescopeMatNullSpaceCreate_CoarseDM()
|
| H A D | telescope.c | 229 …copeSubNullSpaceCreate_Telescope(PC pc, PC_Telescope sred, MatNullSpace nullspace, MatNullSpace *s… in PCTelescopeSubNullSpaceCreate_Telescope() argument
239 PetscCall(MatNullSpaceGetVecs(nullspace, &has_const, &n, &vecs)); in PCTelescopeSubNullSpaceCreate_Telescope()
271 …PetscCheck(!nullspace->remove, PetscObjectComm((PetscObject)pc), PETSC_ERR_SUP, "Propagation of cu… in PCTelescopeSubNullSpaceCreate_Telescope()
272 …PetscCheck(!nullspace->rmctx, PetscObjectComm((PetscObject)pc), PETSC_ERR_SUP, "Propagation of cus… in PCTelescopeSubNullSpaceCreate_Telescope()
285 MatNullSpace nullspace, sub_nullspace; in PCTelescopeMatNullSpaceCreate_default() local
286 PetscCall(MatGetNullSpace(B, &nullspace)); in PCTelescopeMatNullSpaceCreate_default()
287 if (nullspace) { in PCTelescopeMatNullSpaceCreate_default()
289 PetscCall(PCTelescopeSubNullSpaceCreate_Telescope(pc, sred, nullspace, &sub_nullspace)); in PCTelescopeMatNullSpaceCreate_default()
|
| H A D | telescope_dmda.c | 781 …ubNullSpaceCreate_dmda_Telescope(PC pc, PC_Telescope sred, MatNullSpace nullspace, MatNullSpace *s… in PCTelescopeSubNullSpaceCreate_dmda_Telescope() argument
793 PetscCall(MatNullSpaceGetVecs(nullspace, &has_const, &n, &vecs)); in PCTelescopeSubNullSpaceCreate_dmda_Telescope()
830 …PetscCheck(!nullspace->remove, PetscObjectComm((PetscObject)pc), PETSC_ERR_SUP, "Propagation of cu… in PCTelescopeSubNullSpaceCreate_dmda_Telescope()
831 …PetscCheck(!nullspace->rmctx, PetscObjectComm((PetscObject)pc), PETSC_ERR_SUP, "Propagation of cus… in PCTelescopeSubNullSpaceCreate_dmda_Telescope()
843 MatNullSpace nullspace, sub_nullspace; in PCTelescopeMatNullSpaceCreate_dmda() local
844 PetscCall(MatGetNullSpace(B, &nullspace)); in PCTelescopeMatNullSpaceCreate_dmda()
845 if (nullspace) { in PCTelescopeMatNullSpaceCreate_dmda()
847 PetscCall(PCTelescopeSubNullSpaceCreate_dmda_Telescope(pc, sred, nullspace, &sub_nullspace)); in PCTelescopeMatNullSpaceCreate_dmda()
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex73_basic_neumann.out | 3 [size 4] DMCreateMatrix_ShellDA: operator does not have nullspace - attaching
6 [size 4] DMCreateMatrix_ShellDA: operator does not have nullspace - attaching
11 [size 1] DMCreateMatrix_ShellDA: operator does not have nullspace - attaching
|
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 40 …orrectDiscretePressure_Private(SNES snes, PetscInt pfield, MatNullSpace nullspace, Vec u, PetscCtx… in SNESCorrectDiscretePressure_Private() argument
53 …PetscCheck(nullspace, comm, PETSC_ERR_ARG_WRONG, "Cannot compute test without a Jacobian nullspace… in SNESCorrectDiscretePressure_Private()
56 PetscCall(MatNullSpaceGetVecs(nullspace, NULL, &Nv, &nullvecs)); in SNESCorrectDiscretePressure_Private()
116 MatNullSpace nullspace; in SNESConvergedCorrectPressure() local
122 PetscCall(MatGetNullSpace(J, &nullspace)); in SNESConvergedCorrectPressure()
123 PetscCall(MatNullSpaceGetVecs(nullspace, NULL, NULL, &nullvecs)); in SNESConvergedCorrectPressure()
130 MatNullSpace nullspace; in SNESConvergedCorrectPressure() local
135 PetscCall(MatGetNullSpace(J, &nullspace)); in SNESConvergedCorrectPressure()
136 PetscCall(SNESCorrectDiscretePressure_Private(snes, pfield, nullspace, u, ctx)); in SNESConvergedCorrectPressure()
901 MatNullSpace nullspace; in DMSNESCheckJacobian() local
[all …]
|
| /petsc/src/ts/utils/ |
| H A D | dmplexts.c | 365 MatNullSpace nullspace; in DMTSCheckJacobian() local
398 PetscCall(MatGetNullSpace(J, &nullspace)); in DMTSCheckJacobian()
399 if (nullspace) { in DMTSCheckJacobian()
401 PetscCall(MatNullSpaceTest(nullspace, J, &isNull)); in DMTSCheckJacobian()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex62.c | 341 …rorCode CreatePressureNullSpace(DM dm, PetscInt origField, PetscInt field, MatNullSpace *nullspace) in CreatePressureNullSpace() argument
357 PetscCall(MatNullSpaceCreate(PetscObjectComm((PetscObject)dm), PETSC_FALSE, 1, &vec, nullspace)); in CreatePressureNullSpace()
|
| H A D | ex17.c | 601 …rCode CreateElasticityNullSpace(DM dm, PetscInt origField, PetscInt field, MatNullSpace *nullspace) in CreateElasticityNullSpace() argument
604 PetscCall(DMPlexCreateRigidBody(dm, origField, nullspace)); in CreateElasticityNullSpace()
|
| H A D | ex69.c | 3196 …rorCode CreatePressureNullSpace(DM dm, PetscInt origField, PetscInt field, MatNullSpace *nullspace) in CreatePressureNullSpace() argument
3212 PetscCall(MatNullSpaceCreate(PetscObjectComm((PetscObject)dm), PETSC_FALSE, 1, &vec, nullspace)); in CreatePressureNullSpace()
3279 static PetscErrorCode CorrectDiscretePressure(DM dm, MatNullSpace nullspace, Vec u, AppCtx *user) in CorrectDiscretePressure() argument
3290 PetscCall(MatNullSpaceGetVecs(nullspace, NULL, NULL, &nullvecs)); in CorrectDiscretePressure()
3311 MatNullSpace nullspace; in SNESConvergenceCorrectPressure() local
3316 PetscCall(MatGetNullSpace(J, &nullspace)); in SNESConvergenceCorrectPressure()
3317 …PetscCheck(user->solType == SOLZERO || nullspace, PetscObjectComm((PetscObject)snes), PETSC_ERR_AR… in SNESConvergenceCorrectPressure()
3318 if (nullspace) PetscCall(CorrectDiscretePressure(dm, nullspace, u, user)); in SNESConvergenceCorrectPressure()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | KSP.pyx | 1757 the nullspace(Aᵀ), see `Mat.setNullSpace`.
1762 orthogonal to the nullspace(A). The algorithm is simply in each
1763 iteration of the Krylov method we remove the nullspace(A) from
1766 nullspace(A).
1769 systems. If nullspace(A) = nullspace(Aᵀ) (note symmetric
1771 preconditioning will work If nullspace(A) != nullspace(Aᵀ) then
|
| H A D | Mat.pyx | 3753 """Set the nullspace.
3765 """Return the nullspace.
3780 """Set the transpose nullspace.
3792 """Return the transpose nullspace.
3807 """Set the near-nullspace.
3819 """Return the near-nullspace.
6177 """Set the callback to remove the nullspace.
6243 """Return the callback to remove the nullspace.
|
| /petsc/doc/changes/ |
| H A D | 33.md | 27 - MatMult() no longer removes the nullspace set by MatSetNullSpace()
|
| H A D | 38.md | 90 - Removed PCBDDCSetNullSpace. Local nullspace information should now
|
| H A D | 314.md | 294 - Add argument to DM nullspace constructors
|
| /petsc/src/ts/tutorials/ |
| H A D | ex53.c | 1910 …rCode CreateElasticityNullSpace(DM dm, PetscInt origField, PetscInt field, MatNullSpace *nullspace) in CreateElasticityNullSpace() argument
1913 PetscCall(DMPlexCreateRigidBody(dm, origField, nullspace)); in CreateElasticityNullSpace()
|
| /petsc/doc/faq/ |
| H A D | index.md | 1260 ### How can one compute the nullspace of a sparse matrix with MUMPS?
1262 Assuming you have an existing matrix $A$ whose nullspace $V$ you want to find:
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 1526 any of the local problems is singular, the nullspace of the local
|