28   const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0;  in test_3d()  local
41   PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz));  in test_3d()
48         const PetscScalar x = (Lx * i) / M;  in test_3d()
51 …for (c = 0; c < dof; ++c) va[k][j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(…  in test_3d()
71   const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0;  in test_2d()  local
83   PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz));  in test_2d()
89       const PetscScalar x = (Lx * i) / M;  in test_2d()
91 …for (c = 0; c < dof; ++c) va[j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y -…  in test_2d()
110   const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0;  in test_3d_compat()  local
123   PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz));  in test_3d_compat()
134         const PetscScalar x = (Lx * i) / M;  in test_3d_compat()
137 …va[k][j][i]         = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2) + Pe…  in test_3d_compat()
138 …for (c = 0; c < dof; ++c) vVectora[k][j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScal…  in test_3d_compat()
162   const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0;  in test_2d_compat()  local
174   PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz));  in test_2d_compat()
184       const PetscScalar x = (Lx * i) / M;  in test_2d_compat()
186       va[j][i]            = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2);  in test_2d_compat()
187 …for (c = 0; c < dof; ++c) vVectora[j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarI…  in test_2d_compat()