1285 …r *interp_1d, const CeedScalar *grad_1d, const CeedScalar *q_weight_1d, CeedInt P_1d, CeedInt Q_1d,  in CeedBuildMassLaplace()  argument
1287   for (CeedInt i = 0; i < P_1d; i++) {  in CeedBuildMassLaplace()
1288     for (CeedInt j = 0; j < P_1d; j++) {  in CeedBuildMassLaplace()
1290 … (CeedInt k = 0; k < Q_1d; k++) sum += interp_1d[k * P_1d + i] * q_weight_1d[k] * interp_1d[k * P_…  in CeedBuildMassLaplace()
1291       mass[i + j * P_1d] = sum;  in CeedBuildMassLaplace()
1295   for (CeedInt i = 0; i < P_1d; i++) {  in CeedBuildMassLaplace()
1296     for (CeedInt j = 0; j < P_1d; j++) {  in CeedBuildMassLaplace()
1299 …for (CeedInt k = 0; k < Q_1d; k++) sum += grad_1d[k * P_1d + i] * q_weight_1d[k] * grad_1d[k * P_1…  in CeedBuildMassLaplace()
1300       laplace[i + j * P_1d] = sum;  in CeedBuildMassLaplace()
1304   for (CeedInt i = 0; i < P_1d; i++) laplace[i + P_1d * i] += perturbation;  in CeedBuildMassLaplace()
2907 …CeedInt              num_input_fields, P_1d, Q_1d, num_nodes, num_qpts, dim, num_comp = 1, num_ele…  in CeedOperatorCreateFDMElementInverse()  local
2960   CeedCall(CeedBasisGetNumNodes1D(basis, &P_1d));  in CeedOperatorCreateFDMElementInverse()
2971   CeedCall(CeedCalloc(P_1d * P_1d, &mass));  in CeedOperatorCreateFDMElementInverse()
2972   CeedCall(CeedCalloc(P_1d * P_1d, &laplace));  in CeedOperatorCreateFDMElementInverse()
2973   CeedCall(CeedCalloc(P_1d * P_1d, &x));  in CeedOperatorCreateFDMElementInverse()
2974   CeedCall(CeedCalloc(P_1d * P_1d, &fdm_interp));  in CeedOperatorCreateFDMElementInverse()
2975   CeedCall(CeedCalloc(P_1d, &lambda));  in CeedOperatorCreateFDMElementInverse()
2980   CeedCall(CeedBuildMassLaplace(interp_1d, grad_1d, q_weight_1d, P_1d, Q_1d, dim, mass, laplace));  in CeedOperatorCreateFDMElementInverse()
2983   CeedCall(CeedSimultaneousDiagonalization(ceed, laplace, mass, x, lambda, P_1d));  in CeedOperatorCreateFDMElementInverse()
2986   for (CeedInt i = 0; i < P_1d; i++) {  in CeedOperatorCreateFDMElementInverse()
2987     for (CeedInt j = 0; j < P_1d; j++) fdm_interp[i + j * P_1d] = x[j + i * P_1d];  in CeedOperatorCreateFDMElementInverse()
3046             CeedInt i = (n / CeedIntPow(P_1d, d)) % P_1d;  in CeedOperatorCreateFDMElementInverse()
3073     CeedCall(CeedCalloc(P_1d * P_1d, &grad_dummy));  in CeedOperatorCreateFDMElementInverse()
3074     CeedCall(CeedCalloc(P_1d, &q_ref_dummy));  in CeedOperatorCreateFDMElementInverse()
3075     CeedCall(CeedCalloc(P_1d, &q_weight_dummy));  in CeedOperatorCreateFDMElementInverse()
3076 …CeedCall(CeedBasisCreateTensorH1(ceed_parent, dim, num_comp, P_1d, P_1d, fdm_interp, grad_dummy, q…  in CeedOperatorCreateFDMElementInverse()