117   for (d = 0; d < dim; ++d) u[d] = PetscSqr(x[d]) - (d > 0 ? 2.0 * x[d - 1] * x[d] : 0.0);  in quadratic_u()
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()
367   f0[0] -= sum / M - 4 * PetscSqr(PETSC_PI) * kappa * sum * t;  in f0_trig_linear_p()
458     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_u()
462 … (m % 2 == 1) F2 += (8.0 / PetscSqr(m * PETSC_PI)) * PetscCosReal(0.5 * m * PETSC_PI * zstar) * (1…  in terzaghi_2d_u()
495     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_eps()
499 …* L)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * (1.0 - PetscExpReal(-PetscSqr(m * PETSC_PI) * t…  in terzaghi_2d_eps()
531     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_p()
537 …* PETSC_PI)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * PetscExpReal(-PetscSqr(m * PETSC_PI) * t…  in terzaghi_2d_p()
570     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_u_t()
574 …if (m % 2 == 1) F2_t += (2.0 * c / PetscSqr(L)) * PetscCosReal(0.5 * m * PETSC_PI * zstar) * Petsc…  in terzaghi_2d_u_t()
607     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_eps_t()
611 …(L * L * L)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * PetscExpReal(-PetscSqr(m * PETSC_PI) * t…  in terzaghi_2d_eps_t()
638     u[0] = -((P_0 * eta) / (G * S)) * PetscSqr(0 * PETSC_PI) * c / PetscSqr(2.0 * L); /* Pa / s */  in terzaghi_2d_p_t()
655     PetscReal   tstar = PetscRealPart(c * time) / PetscSqr(2.0 * L); /* - */  in terzaghi_2d_p_t()
661 …= 1) F1_t += ((-m * PETSC_PI * c) / PetscSqr(L)) * PetscSinReal(0.5 * m * PETSC_PI * zstar) * Pets…  in terzaghi_2d_p_t()
1126     PetscReal   tstar  = PetscRealPart(c * time) / PetscSqr(R_0); /* - */  in cryer_3d_u()
1132 …const PetscReal E_n = PetscRealPart(PetscSqr(1 - nu) * PetscSqr(1 + nu_u) * x_n - 18.0 * (1 + nu) …  in cryer_3d_u()
1136 …A_n += PetscRealPart((12.0 * (1.0 + nu) * (nu_u - nu)) / ((1.0 - 2.0 * nu) * E_n * PetscSqr(R_star…  in cryer_3d_u()
1176     PetscReal tstar  = PetscRealPart(c * time) / PetscSqr(R_0); /* - */  in cryer_3d_eps()
1182 …const PetscReal E_n = PetscRealPart(PetscSqr(1 - nu) * PetscSqr(1 + nu_u) * x_n - 18.0 * (1 + nu) …  in cryer_3d_eps()
1184 ….0 * nu) * E_n * PetscSqr(R_star) * x_n * PetscSinReal(PetscSqrtReal(x_n))) * (3.0 * (nu_u - nu) *…  in cryer_3d_eps()
1189 …const PetscReal E_n = PetscRealPart(PetscSqr(1 - nu) * PetscSqr(1 + nu_u) * x_n - 18.0 * (1 + nu) …  in cryer_3d_eps()
1191 …cRealPart((12.0 * (1.0 + nu) * (nu_u - nu)) / ((1.0 - 2.0 * nu) * E_n * PetscSqr(R_star) * x_n * P…  in cryer_3d_eps()
1228     PetscReal t_star = PetscRealPart(c * time) / PetscSqr(R_0);  in cryer_3d_p()
1233 …const PetscReal E_n = PetscRealPart(PetscSqr(1 - nu) * PetscSqr(1 + nu_u) * x_n - 18.0 * (1 + nu) …  in cryer_3d_p()
1235 …A_x += PetscRealPart(((18.0 * PetscSqr(nu_u - nu)) / (eta * E_n)) * (PetscSinReal(R_star * PetscSq…  in cryer_3d_p()
1586       a1 = PetscSqr(n * PETSC_PI) - k * PETSC_PI;  in cryerZeros()
1587       a2 = PetscSqr(n * PETSC_PI) + k * PETSC_PI;  in cryerZeros()
1670       ctx->t_r = PetscSqr(ctx->xmax[0] - ctx->xmin[0]) / c;  in SetupParameters()
1673       ctx->t_r = PetscSqr(2.0 * (ctx->xmax[1] - ctx->xmin[1])) / c;  in SetupParameters()
1676       ctx->t_r = PetscSqr(2.0 * (ctx->xmax[1] - ctx->xmin[1])) / c;  in SetupParameters()
1679       ctx->t_r = PetscSqr(ctx->xmax[1]) / c;  in SetupParameters()