125   PetscScalar   *u, *uj, xj, xi;  in FVSample_3WaySplit()  local
135   PetscCall(PetscMalloc1(dof, &uj));  in FVSample_3WaySplit()
147 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj));  in FVSample_3WaySplit()
148         for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N;  in FVSample_3WaySplit()
156 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj));  in FVSample_3WaySplit()
157         for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N;  in FVSample_3WaySplit()
165 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj));  in FVSample_3WaySplit()
166         for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N;  in FVSample_3WaySplit()
174 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj));  in FVSample_3WaySplit()
175         for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N;  in FVSample_3WaySplit()
183 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj));  in FVSample_3WaySplit()
184         for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N;  in FVSample_3WaySplit()
189   PetscCall(PetscFree(uj));  in FVSample_3WaySplit()