15 …rorCode DMSwarmCeedContextCreate(DM dm_swarm, const char *ceed_resource, DMSwarmCeedContext *ctx) {  in DMSwarmCeedContextCreate()  argument
23   PetscCall(PetscNew(ctx));  in DMSwarmCeedContextCreate()
27   CeedInit(ceed_resource, &(*ctx)->ceed);  in DMSwarmCeedContextCreate()
32     PetscCall(CreateBasisFromPlex((*ctx)->ceed, dm_mesh, NULL, 0, 0, 0, bp_data, &basis_u));  in DMSwarmCeedContextCreate()
33     PetscCall(CreateBasisFromPlex((*ctx)->ceed, dm_coord, NULL, 0, 0, 0, bp_data, &basis_x));  in DMSwarmCeedContextCreate()
34     PetscCall(CreateRestrictionFromPlex((*ctx)->ceed, dm_mesh, 0, NULL, 0, &elem_restr_u_mesh));  in DMSwarmCeedContextCreate()
35     PetscCall(CreateRestrictionFromPlex((*ctx)->ceed, dm_coord, 0, NULL, 0, &elem_restr_x_mesh));  in DMSwarmCeedContextCreate()
38     CeedElemRestrictionCreateVector(elem_restr_u_mesh, &(*ctx)->u_mesh, NULL);  in DMSwarmCeedContextCreate()
39     CeedElemRestrictionCreateVector(elem_restr_u_mesh, &(*ctx)->v_mesh, NULL);  in DMSwarmCeedContextCreate()
63 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, num_comp, num_points * num_c…  in DMSwarmCeedContextCreate()
65 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, dim, num_points * dim, CEED_…  in DMSwarmCeedContextCreate()
67 …CeedElemRestrictionCreateAtPoints((*ctx)->ceed, num_elem, num_points, 1, num_points, CEED_MEM_HOST…  in DMSwarmCeedContextCreate()
71     CeedElemRestrictionCreateVector(elem_restr_u_points, &(*ctx)->u_points, NULL);  in DMSwarmCeedContextCreate()
79       CeedVectorCreate((*ctx)->ceed, num_points * dim, &x_ref_points);  in DMSwarmCeedContextCreate()
99         CeedVectorCreate((*ctx)->ceed, len, &x_coord);  in DMSwarmCeedContextCreate()
107       CeedQFunctionCreateInterior((*ctx)->ceed, 1, SetupMass, SetupMass_loc, &qf_setup);  in DMSwarmCeedContextCreate()
112 …CeedOperatorCreateAtPoints((*ctx)->ceed, qf_setup, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &op_s…  in DMSwarmCeedContextCreate()
131   PetscCall(DMSetApplicationContext(dm_mesh, (void *)(*ctx)));  in DMSwarmCeedContextCreate()
139 …CeedQFunctionCreateIdentity((*ctx)->ceed, num_comp, CEED_EVAL_INTERP, CEED_EVAL_NONE, &qf_mesh_to_…  in DMSwarmCeedContextCreate()
141 …CeedOperatorCreateAtPoints((*ctx)->ceed, qf_mesh_to_points, NULL, NULL, &(*ctx)->op_mesh_to_points…  in DMSwarmCeedContextCreate()
142 …CeedOperatorSetField((*ctx)->op_mesh_to_points, "input", elem_restr_u_mesh, basis_u, CEED_VECTOR_A…  in DMSwarmCeedContextCreate()
143 …CeedOperatorSetField((*ctx)->op_mesh_to_points, "output", elem_restr_u_points, CEED_BASIS_NONE, CE…  in DMSwarmCeedContextCreate()
144     CeedOperatorAtPointsSetPoints((*ctx)->op_mesh_to_points, elem_restr_x_points, x_ref_points);  in DMSwarmCeedContextCreate()
156     CeedQFunctionCreateInterior((*ctx)->ceed, 1, Mass, Mass_loc, &qf_pts_to_mesh);  in DMSwarmCeedContextCreate()
162     CeedQFunctionContextCreate((*ctx)->ceed, &qf_ctx);  in DMSwarmCeedContextCreate()
167 …CeedOperatorCreateAtPoints((*ctx)->ceed, qf_pts_to_mesh, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE,…  in DMSwarmCeedContextCreate()
168 …CeedOperatorSetField((*ctx)->op_points_to_mesh, "q data", elem_restr_q_data_points, CEED_BASIS_NON…  in DMSwarmCeedContextCreate()
169 …CeedOperatorSetField((*ctx)->op_points_to_mesh, "u", elem_restr_u_points, CEED_BASIS_NONE, CEED_VE…  in DMSwarmCeedContextCreate()
170 …CeedOperatorSetField((*ctx)->op_points_to_mesh, "v", elem_restr_u_mesh, basis_u, CEED_VECTOR_ACTIV…  in DMSwarmCeedContextCreate()
171     CeedOperatorAtPointsSetPoints((*ctx)->op_points_to_mesh, elem_restr_x_points, x_ref_points);  in DMSwarmCeedContextCreate()
184     CeedQFunctionCreateInterior((*ctx)->ceed, 1, Mass, Mass_loc, &qf_mass);  in DMSwarmCeedContextCreate()
190     CeedQFunctionContextCreate((*ctx)->ceed, &ctx_mass);  in DMSwarmCeedContextCreate()
195 …CeedOperatorCreateAtPoints((*ctx)->ceed, qf_mass, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &(*ctx…  in DMSwarmCeedContextCreate()
196 …CeedOperatorSetField((*ctx)->op_mass, "q data", elem_restr_q_data_points, CEED_BASIS_NONE, q_data_…  in DMSwarmCeedContextCreate()
197     CeedOperatorSetField((*ctx)->op_mass, "u", elem_restr_u_mesh, basis_u, CEED_VECTOR_ACTIVE);  in DMSwarmCeedContextCreate()
198     CeedOperatorSetField((*ctx)->op_mass, "v", elem_restr_u_mesh, basis_u, CEED_VECTOR_ACTIVE);  in DMSwarmCeedContextCreate()
199     CeedOperatorAtPointsSetPoints((*ctx)->op_mass, elem_restr_x_points, x_ref_points);  in DMSwarmCeedContextCreate()
219 PetscErrorCode DMSwarmCeedContextDestroy(DMSwarmCeedContext *ctx) {  in DMSwarmCeedContextDestroy()  argument
221   CeedDestroy(&(*ctx)->ceed);  in DMSwarmCeedContextDestroy()
222   CeedVectorDestroy(&(*ctx)->u_mesh);  in DMSwarmCeedContextDestroy()
223   CeedVectorDestroy(&(*ctx)->v_mesh);  in DMSwarmCeedContextDestroy()
224   CeedVectorDestroy(&(*ctx)->u_points);  in DMSwarmCeedContextDestroy()
225   CeedOperatorDestroy(&(*ctx)->op_mesh_to_points);  in DMSwarmCeedContextDestroy()
226   CeedOperatorDestroy(&(*ctx)->op_points_to_mesh);  in DMSwarmCeedContextDestroy()
227   CeedOperatorDestroy(&(*ctx)->op_mass);  in DMSwarmCeedContextDestroy()
228   PetscCall(PetscFree(*ctx));  in DMSwarmCeedContextDestroy()