Lines Matching refs:PETSc

11 from petsc4py import PETSc
16 …sizes = PETSc.IS().load(PETSc.Viewer().createBinary(f"{prefix}{identifier}_sizes_{rank}_{size}.dat…
24   A = A.load(PETSc.Viewer().createBinary(f"{prefix}{identifier}.dat", "r", comm = PETSc.COMM_WORLD))
27 …aux_IS.load(PETSc.Viewer().createBinary(f"{prefix}{identifier}_is_{rank}_{size}.dat", "r", comm = …
29 …aux_Mat.load(PETSc.Viewer().createBinary(f"{prefix}{identifier}_aux_{rank}_{size}.dat", "r", comm …
32 size = PETSc.COMM_WORLD.getSize()
34 rank = PETSc.COMM_WORLD.getRank()
41 system_str = PETSc.Options().getString("system", "elasticity")
46   empty_A11 = PETSc.Options().getBool("empty_A11", False)
55 A[0] = PETSc.Mat().create(comm = PETSc.COMM_WORLD)
57 aux_IS[0] = PETSc.IS().create(comm = PETSc.COMM_SELF)
58 aux_Mat[0] = PETSc.Mat().create(comm = PETSc.COMM_SELF)
59 A[3] = PETSc.Mat().create(comm = PETSc.COMM_WORLD)
61 aux_IS[1] = PETSc.IS().create(comm = PETSc.COMM_SELF)
62 aux_Mat[1] = PETSc.Mat().create(comm = PETSc.COMM_SELF)
65 load_dir = PETSc.Options().getString("load_dir", "${DATAFILESPATH}/matrices/hpddm/GENEO")
74 A[2] = PETSc.Mat().create(comm = PETSc.COMM_WORLD)
82 A[2].load(PETSc.Viewer().createBinary(f"{load_dir}/{ 'B' if id_sys == 1 else 'A' }10.dat", "r", com…
84 A[1] = PETSc.Mat().createTranspose(A[2])
87 S = PETSc.Mat().createNest([[A[0], A[1]], [A[2], A[3] if not empty_A11 else None]])
89 ksp = PETSc.KSP().create(comm = PETSc.COMM_WORLD)
94 pc.setType(PETSc.PC.Type.FIELDSPLIT)
96 pc.setFieldSplitType(PETSc.PC.CompositeType.SCHUR)
98 pc.setFieldSplitSchurPreType(PETSc.PC.FieldSplitSchurPreType.SELF)
110 pc0.setType(PETSc.PC.Type.HPDDM)
118 pc1.setType(PETSc.PC.Type.HPDDM)
127 b.load(PETSc.Viewer().createBinary(f"{load_dir}/rhs_{ 'B' if id_sys == 1 else 'A' }.dat", "r", comm…