Lines Matching refs:PetscCallA

25   PetscCallA(PetscInitialize(PETSC_NULL_CHARACTER, ierr))
26   PetscCallA(DMCreate(PETSC_COMM_WORLD, dm, ierr))
27   PetscCallA(DMSetType(dm, DMPLEX, ierr))
28   PetscCallA(DMSetFromOptions(dm, ierr))
29   PetscCallA(DMGetDimension(dm, dim, ierr))
30   PetscCallA(DMViewFromOptions(dm, PETSC_NULL_OBJECT, '-dm_view', ierr))
33 …PetscCallA(PetscFECreateLagrange(PETSC_COMM_SELF, dim, Nc, PETSC_FALSE, degree, PETSC_DETERMINE, f…
34   PetscCallA(DMSetField(dm, field, PETSC_NULL_DMLABEL, PetscObjectCast(fe), ierr))
35   PetscCallA(DMCreateDS(dm, ierr))
36   PetscCallA(PetscFEDestroy(fe, ierr))
39   PetscCallA(DMCreate(PETSC_COMM_WORLD, sw, ierr))
40   PetscCallA(DMSetType(sw, DMSWARM, ierr))
41   PetscCallA(DMSetDimension(sw, dim, ierr))
42   PetscCallA(DMSwarmSetType(sw, DMSWARM_PIC, ierr))
43   PetscCallA(DMSwarmSetCellDM(sw, dm, ierr))
44   PetscCallA(DMSwarmRegisterPetscDatatypeField(sw, 'w_q', Nc, PETSC_SCALAR, ierr))
45   PetscCallA(DMSwarmFinalizeFieldRegister(sw, ierr))
46   PetscCallA(DMSwarmSetLocalSizes(sw, Np, zero, ierr))
47   PetscCallA(DMSetFromOptions(sw, ierr))
48   PetscCallA(DMSwarmGetField(sw, 'w_q', bs, dtype, wq, ierr))
49   PetscCallA(DMSwarmGetField(sw, 'DMSwarmPIC_coor', bs, dtype, coords, ierr))
55   PetscCallA(DMSwarmRestoreField(sw, 'DMSwarmPIC_coor', bs, dtype, coords, ierr))
56   PetscCallA(DMSwarmRestoreField(sw, 'w_q', bs, dtype, wq, ierr))
57   PetscCallA(DMSwarmMigrate(sw, removePoints, ierr))
58   PetscCallA(DMSwarmVectorDefineField(sw, 'w_q', ierr))
59   PetscCallA(DMViewFromOptions(sw, PETSC_NULL_OBJECT, '-swarm_view', ierr))
63   PetscCallA(DMCreateMassMatrix(sw, dm, M_p, ierr))
64   PetscCallA(DMCreateGlobalVector(dm, rho, ierr))
65   PetscCallA(DMSwarmCreateGlobalVectorFromField(sw, 'w_q', f, ierr))
66   PetscCallA(MatMultTranspose(M_p, f, rho, ierr))
69   PetscCallA(DMSetOutputSequenceNumber(dm, timestep, time, ierr))
70   PetscCallA(PetscObjectViewFromOptions(PetscObjectCast(rho), PETSC_NULL_OBJECT, '-rho_view', ierr))
74   PetscCallA(DMCreateMassMatrix(dm, dm, M, ierr))
75   PetscCallA(DMCreateGlobalVector(dm, rhs, ierr))
77     PetscCallA(MatMult(M, rho, rhs, ierr)) ! this is what you would do for and FE solve
79     PetscCallA(VecCopy(rho, rhs, ierr)) ! Identity: M^1 M rho
81   PetscCallA(KSPCreate(PETSC_COMM_WORLD, ksp, ierr))
82   PetscCallA(KSPSetOptionsPrefix(ksp, 'ftop_', ierr))
83   PetscCallA(KSPSetFromOptions(ksp, ierr))
84   PetscCallA(KSPSetOperators(ksp, M_p, M_p, ierr))
85   PetscCallA(KSPSolveTranspose(ksp, rhs, f, ierr))
86   PetscCallA(KSPDestroy(ksp, ierr))
87   PetscCallA(VecDestroy(rhs, ierr))
88   PetscCallA(MatDestroy(M_p, ierr))
89   PetscCallA(MatDestroy(M, ierr))
92   PetscCallA(DMSetOutputSequenceNumber(sw, timestep, time, ierr))
93 …PetscCallA(PetscObjectViewFromOptions(PetscObjectCast(f), PETSC_NULL_OBJECT, '-weights_view', ierr…
94   PetscCallA(VecNorm(f, NORM_1, norm, ierr))
97   PetscCallA(DMSwarmDestroyGlobalVectorFromField(sw, 'w_q', f, ierr))
98   PetscCallA(VecDestroy(rho, ierr))
99   PetscCallA(DMDestroy(sw, ierr))
100   PetscCallA(DMDestroy(dm, ierr))
102   PetscCallA(PetscFinalize(ierr))