147 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_quadratic_linear_u()  argument
157   for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * t;  in f0_quadratic_linear_u()
158   f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * t;  in f0_quadratic_linear_u()
161 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_quadratic_linear_p()  argument
169   f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0;  in f0_quadratic_linear_p()
170   f0[0] += u_t ? u_t[uOff[2]] / M : 0.0;  in f0_quadratic_linear_p()
171   f0[0] -= sum / M;  in f0_quadratic_linear_p()
229 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_quadratic_trig_u()  argument
239   for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * PetscCosReal(t);  in f0_quadratic_trig_u()
240   f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * PetscCosReal(t);  in f0_quadratic_trig_u()
243 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_quadratic_trig_p()  argument
252   f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0;  in f0_quadratic_trig_p()
253   f0[0] += u_t ? u_t[uOff[2]] / M : 0.0;  in f0_quadratic_trig_p()
254   f0[0] += PetscSinReal(t) * sum / M;  in f0_quadratic_trig_p()
342 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_trig_linear_u()  argument
352 …for (d = 0; d < dim - 1; ++d) f0[d] += PetscSqr(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[d]…  in f0_trig_linear_u()
353 …f0[dim - 1] += PetscSqr(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[dim - 1]) * (2. * G + lamb…  in f0_trig_linear_u()
356 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_trig_linear_p()  argument
365   f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0;  in f0_trig_linear_p()
366   f0[0] += u_t ? u_t[uOff[2]] / M : 0.0;  in f0_trig_linear_p()
367   f0[0] -= sum / M - 4 * PetscSqr(PETSC_PI) * kappa * sum * t;  in f0_trig_linear_p()
1243 …l x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_terzaghi_bd_u()  argument
1247   f0[0] = 0.0;  in f0_terzaghi_bd_u()
1248   f0[1] = P;  in f0_terzaghi_bd_u()
1255 …l x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])
1345     f0[1] += sigma_zz;
1347     f0[1] -= sigma_zz;
1352 …l x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_cryer_bd_u()  argument
1357   for (d = 0; d < dim; ++d) f0[d] = -P_0 * n[d];  in f0_cryer_bd_u()
1362 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_epsilon()  argument
1366   for (d = 0; d < dim; ++d) f0[0] += u_x[d * dim + d];  in f0_epsilon()
1367   f0[0] -= u[uOff[1]];  in f0_epsilon()
1371 …eal t, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[])  in f0_p()  argument
1376   f0[0] += alpha * u_t[uOff[1]];  in f0_p()
1377   f0[0] += u_t[uOff[2]] / M;  in f0_p()
1378   if (f0[0] != f0[0]) abort();  in f0_p()