72   PetscInt         d, q, npoints_q;  in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()  local
78   for (q = 0; q < npoints_q; q++) {  in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
79     for (d = 0; d < dim; d++) xi[dim * q + d] = quadrature_xi[dim * q + d];  in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
91   PetscInt           nel, npe, e, q, k, d;  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()  local
132   for (q = 0; q < npoints_q; q++) {  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
133     PetscCall(PetscMalloc1(npe, &basis[q]));  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
137       basis[q][0] = 0.5 * (1.0 - xi[dim * q + 0]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
138       basis[q][1] = 0.5 * (1.0 + xi[dim * q + 0]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
141       basis[q][0] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
142       basis[q][1] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
143       basis[q][2] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
144       basis[q][3] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
148 … basis[q][0] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
149 … basis[q][1] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
150 … basis[q][2] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
151 … basis[q][3] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
152 … basis[q][4] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
153 … basis[q][5] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
154 … basis[q][6] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
155 … basis[q][7] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]);  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
179     for (q = 0; q < npoints_q; q++) {  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
182         for (d = 0; d < dim; d++) swarm_coor[dim * pcnt + d] += basis[q][k] * elcoor[dim * k + d];  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
195   for (q = 0; q < npoints_q; q++) PetscCall(PetscFree(basis[q]));  in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()