37 PetscErrorCode CeedDataDestroy(CeedInt level, CeedData data) {  in CeedDataDestroy()  argument
41   CeedVectorDestroy(&data->x_ceed);  in CeedDataDestroy()
42   CeedVectorDestroy(&data->y_ceed);  in CeedDataDestroy()
43   CeedVectorDestroy(&data->geo_data);  in CeedDataDestroy()
44   for (CeedInt i = 0; i < SOLIDS_MAX_NUMBER_FIELDS; i++) CeedVectorDestroy(&data->stored_fields[i]);  in CeedDataDestroy()
45   CeedVectorDestroy(&data->geo_data_diagnostic);  in CeedDataDestroy()
46   CeedVectorDestroy(&data->true_soln);  in CeedDataDestroy()
48   CeedElemRestrictionDestroy(&data->elem_restr_x);  in CeedDataDestroy()
49   CeedElemRestrictionDestroy(&data->elem_restr_u);  in CeedDataDestroy()
50   CeedElemRestrictionDestroy(&data->elem_restr_geo_data_i);  in CeedDataDestroy()
51 …for (CeedInt i = 0; i < SOLIDS_MAX_NUMBER_FIELDS; i++) CeedElemRestrictionDestroy(&data->elem_rest…  in CeedDataDestroy()
52   CeedElemRestrictionDestroy(&data->elem_restr_energy);  in CeedDataDestroy()
53   CeedElemRestrictionDestroy(&data->elem_restr_diagnostic);  in CeedDataDestroy()
54   CeedElemRestrictionDestroy(&data->elem_restr_geo_data_diagnostic_i);  in CeedDataDestroy()
56   CeedBasisDestroy(&data->basis_x);  in CeedDataDestroy()
57   CeedBasisDestroy(&data->basis_u);  in CeedDataDestroy()
58   CeedBasisDestroy(&data->basis_energy);  in CeedDataDestroy()
59   CeedBasisDestroy(&data->basis_diagnostic);  in CeedDataDestroy()
61   CeedQFunctionDestroy(&data->qf_residual);  in CeedDataDestroy()
62   CeedQFunctionDestroy(&data->qf_jacobian);  in CeedDataDestroy()
63   CeedQFunctionDestroy(&data->qf_energy);  in CeedDataDestroy()
64   CeedQFunctionDestroy(&data->qf_diagnostic);  in CeedDataDestroy()
66   CeedOperatorDestroy(&data->op_residual);  in CeedDataDestroy()
67   CeedOperatorDestroy(&data->op_jacobian);  in CeedDataDestroy()
68   CeedOperatorDestroy(&data->op_energy);  in CeedDataDestroy()
69   CeedOperatorDestroy(&data->op_diagnostic);  in CeedDataDestroy()
71   CeedBasisDestroy(&data->basis_c_to_f);  in CeedDataDestroy()
72   CeedOperatorDestroy(&data->op_prolong);  in CeedDataDestroy()
73   CeedOperatorDestroy(&data->op_restrict);  in CeedDataDestroy()
75   PetscCall(PetscFree(data));  in CeedDataDestroy()
127 …                                  CeedVector force_ceed, CeedVector neumann_ceed, CeedData *data) {  in SetupLibceedFineLevel()  argument
150 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_x, 2, Q, problem_data.quadrature_mode, &data[f…  in SetupLibceedFineLevel()
152 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_u, P, Q, problem_data.quadrature_mode, &data[f…  in SetupLibceedFineLevel()
154 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_e, P, Q, problem_data.quadrature_mode, &data[f…  in SetupLibceedFineLevel()
156 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_u, P, P, CEED_GAUSS_LOBATTO, &data[fine_level]…  in SetupLibceedFineLevel()
165   PetscCall(CreateRestrictionFromPlex(ceed, dm_coord, 0, 0, 0, &(data[fine_level]->elem_restr_x)));  in SetupLibceedFineLevel()
167   PetscCall(CreateRestrictionFromPlex(ceed, dm, 0, 0, 0, &data[fine_level]->elem_restr_u));  in SetupLibceedFineLevel()
169 …PetscCall(CreateRestrictionFromPlex(ceed, dm_energy, 0, 0, 0, &data[fine_level]->elem_restr_energy…  in SetupLibceedFineLevel()
171 …PetscCall(CreateRestrictionFromPlex(ceed, dm_diagnostic, 0, 0, 0, &data[fine_level]->elem_restr_di…  in SetupLibceedFineLevel()
176   CeedBasisGetNumQuadraturePoints(data[fine_level]->basis_u, &num_qpts);  in SetupLibceedFineLevel()
179                                    &data[fine_level]->elem_restr_geo_data_i);  in SetupLibceedFineLevel()
183 …                           CEED_STRIDES_BACKEND, &data[fine_level]->elem_restr_stored_fields_i[i]);  in SetupLibceedFineLevel()
187                                    &data[fine_level]->elem_restr_geo_data_diagnostic_i);  in SetupLibceedFineLevel()
195   CeedElemRestrictionCreateVector(data[fine_level]->elem_restr_x, &x_coord, NULL);  in SetupLibceedFineLevel()
203   CeedVectorCreate(ceed, U_loc_size, &data[fine_level]->x_ceed);  in SetupLibceedFineLevel()
204   CeedVectorCreate(ceed, U_loc_size, &data[fine_level]->y_ceed);  in SetupLibceedFineLevel()
206   CeedVectorCreate(ceed, num_elem * num_qpts * q_data_size, &data[fine_level]->geo_data);  in SetupLibceedFineLevel()
209 …CeedVectorCreate(ceed, num_elem * num_qpts * problem_data.field_sizes[i], &data[fine_level]->store…  in SetupLibceedFineLevel()
212 …CeedVectorCreate(ceed, num_elem * P * P * P * q_data_size, &data[fine_level]->geo_data_diagnostic);  in SetupLibceedFineLevel()
226 …CeedOperatorSetField(op_setup_geo, "dx", data[fine_level]->elem_restr_x, data[fine_level]->basis_x…  in SetupLibceedFineLevel()
227 …CeedOperatorSetField(op_setup_geo, "weight", CEED_ELEMRESTRICTION_NONE, data[fine_level]->basis_x,…  in SetupLibceedFineLevel()
228 …CeedOperatorSetField(op_setup_geo, "qdata", data[fine_level]->elem_restr_geo_data_i, CEED_BASIS_NO…  in SetupLibceedFineLevel()
230   CeedOperatorApply(op_setup_geo, x_coord, data[fine_level]->geo_data, CEED_REQUEST_IMMEDIATE);  in SetupLibceedFineLevel()
251 …CeedOperatorSetField(op_residual, "du", data[fine_level]->elem_restr_u, data[fine_level]->basis_u,…  in SetupLibceedFineLevel()
252 …CeedOperatorSetField(op_residual, "qdata", data[fine_level]->elem_restr_geo_data_i, CEED_BASIS_NON…  in SetupLibceedFineLevel()
253 …CeedOperatorSetField(op_residual, "dv", data[fine_level]->elem_restr_u, data[fine_level]->basis_u,…  in SetupLibceedFineLevel()
255 …CeedOperatorSetField(op_residual, problem_data.field_names[i], data[fine_level]->elem_restr_stored…  in SetupLibceedFineLevel()
256                          data[fine_level]->stored_fields[i]);  in SetupLibceedFineLevel()
259   data[fine_level]->qf_residual = qf_residual;  in SetupLibceedFineLevel()
260   data[fine_level]->op_residual = op_residual;  in SetupLibceedFineLevel()
278 …CeedOperatorSetField(op_jacobian, "delta du", data[fine_level]->elem_restr_u, data[fine_level]->ba…  in SetupLibceedFineLevel()
279 …CeedOperatorSetField(op_jacobian, "qdata", data[fine_level]->elem_restr_geo_data_i, CEED_BASIS_NON…  in SetupLibceedFineLevel()
280 …CeedOperatorSetField(op_jacobian, "delta dv", data[fine_level]->elem_restr_u, data[fine_level]->ba…  in SetupLibceedFineLevel()
282 …CeedOperatorSetField(op_jacobian, problem_data.field_names[i], data[fine_level]->elem_restr_stored…  in SetupLibceedFineLevel()
283                          data[fine_level]->stored_fields[i]);  in SetupLibceedFineLevel()
286   data[fine_level]->qf_jacobian = qf_jacobian;  in SetupLibceedFineLevel()
287   data[fine_level]->op_jacobian = op_jacobian;  in SetupLibceedFineLevel()
366 …CeedOperatorSetField(op_setup_force, "x", data[fine_level]->elem_restr_x, data[fine_level]->basis_…  in SetupLibceedFineLevel()
367 …eedOperatorSetField(op_setup_force, "qdata", data[fine_level]->elem_restr_geo_data_i, CEED_BASIS_N…  in SetupLibceedFineLevel()
368 …CeedOperatorSetField(op_setup_force, "force", data[fine_level]->elem_restr_u, data[fine_level]->ba…  in SetupLibceedFineLevel()
390     CeedVectorCreate(ceed, U_loc_size, &(data[fine_level]->true_soln));  in SetupLibceedFineLevel()
399 …CeedOperatorSetField(op_true, "x", data[fine_level]->elem_restr_x, basis_x_true, CEED_VECTOR_ACTIV…  in SetupLibceedFineLevel()
400 …CeedOperatorSetField(op_true, "true solution", data[fine_level]->elem_restr_u, CEED_BASIS_NONE, CE…  in SetupLibceedFineLevel()
402     CeedOperatorApply(op_true, x_coord, data[fine_level]->true_soln, CEED_REQUEST_IMMEDIATE);  in SetupLibceedFineLevel()
404     CeedElemRestrictionCreateVector(data[fine_level]->elem_restr_u, &mult_vec, NULL);  in SetupLibceedFineLevel()
406     CeedElemRestrictionGetMultiplicity(data[fine_level]->elem_restr_u, mult_vec);  in SetupLibceedFineLevel()
408     CeedVectorGetArray(data[fine_level]->true_soln, CEED_MEM_HOST, &true_array);  in SetupLibceedFineLevel()
411     CeedVectorRestoreArray(data[fine_level]->true_soln, &true_array);  in SetupLibceedFineLevel()
433 …CeedOperatorSetField(op_energy, "du", data[fine_level]->elem_restr_u, data[fine_level]->basis_u, C…  in SetupLibceedFineLevel()
434 …CeedOperatorSetField(op_energy, "qdata", data[fine_level]->elem_restr_geo_data_i, CEED_BASIS_NONE,…  in SetupLibceedFineLevel()
435 …CeedOperatorSetField(op_energy, "energy", data[fine_level]->elem_restr_energy, data[fine_level]->b…  in SetupLibceedFineLevel()
437   data[fine_level]->qf_energy = qf_energy;  in SetupLibceedFineLevel()
438   data[fine_level]->op_energy = op_energy;  in SetupLibceedFineLevel()
456 …CeedOperatorSetField(op_setup_geo, "dx", data[fine_level]->elem_restr_x, basis_x, CEED_VECTOR_ACTI…  in SetupLibceedFineLevel()
458 …CeedOperatorSetField(op_setup_geo, "qdata", data[fine_level]->elem_restr_geo_data_diagnostic_i, CE…  in SetupLibceedFineLevel()
460 …CeedOperatorApply(op_setup_geo, x_coord, data[fine_level]->geo_data_diagnostic, CEED_REQUEST_IMMED…  in SetupLibceedFineLevel()
476 …CeedOperatorSetField(op_diagnostic, "u", data[fine_level]->elem_restr_u, data[fine_level]->basis_d…  in SetupLibceedFineLevel()
477 …CeedOperatorSetField(op_diagnostic, "du", data[fine_level]->elem_restr_u, data[fine_level]->basis_…  in SetupLibceedFineLevel()
478 …CeedOperatorSetField(op_diagnostic, "qdata", data[fine_level]->elem_restr_geo_data_diagnostic_i, C…  in SetupLibceedFineLevel()
479                        data[fine_level]->geo_data_diagnostic);  in SetupLibceedFineLevel()
480 …CeedOperatorSetField(op_diagnostic, "diagnostic values", data[fine_level]->elem_restr_diagnostic, …  in SetupLibceedFineLevel()
482   data[fine_level]->qf_diagnostic = qf_diagnostic;  in SetupLibceedFineLevel()
483   data[fine_level]->op_diagnostic = op_diagnostic;  in SetupLibceedFineLevel()
496                                  PetscInt U_loc_size, CeedVector fine_mult, CeedData *data) {  in SetupLibceedLevel()  argument
511   PetscCall(CreateRestrictionFromPlex(ceed, dm, 0, 0, 0, &data[level]->elem_restr_u));  in SetupLibceedLevel()
517 …CeedBasisCreateTensorH1Lagrange(ceed, dim, num_comp_u, P, Q, problem_data.quadrature_mode, &data[l…  in SetupLibceedLevel()
522   CeedVectorCreate(ceed, U_loc_size, &data[level]->x_ceed);  in SetupLibceedLevel()
523   CeedVectorCreate(ceed, U_loc_size, &data[level]->y_ceed);  in SetupLibceedLevel()
530 …CeedOperatorMultigridLevelCreate(data[level + 1]->op_jacobian, fine_mult, data[level]->elem_restr_…  in SetupLibceedLevel()
534   data[level]->op_jacobian     = op_jacobian;  in SetupLibceedLevel()
535   data[level + 1]->op_prolong  = op_prolong;  in SetupLibceedLevel()
536   data[level + 1]->op_restrict = op_restrict;  in SetupLibceedLevel()