128 static inline PetscReal FIX_C_3d(PetscScalar C00, PetscScalar C01, PetscScalar C02, PetscScalar C11…  in FIX_C_3d()  argument
132   PetscBool twod = (PetscBool)(C02 == 0 && C12 == 0 && C22 == 0);  in FIX_C_3d()
133   Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs);  in FIX_C_3d()
178   const PetscScalar  C02      = split ? a[aOff[C_FIELD_ID] + 2] : u[uOff[C_FIELD_ID] + 2];  in C_0_3d()  local
182   const PetscScalar  norm     = NORM2C_3d(C00, C01, C02, C11, C12, C22) + eps;  in C_0_3d()
223   const PetscScalar C02    = split ? a[aOff[C_FIELD_ID] + 2] : u[uOff[C_FIELD_ID] + 2];  in JC_0_c0c0_3d()  local
227   const PetscScalar norm   = NORM2C_3d(C00, C01, C02, C11, C12, C22) + eps;  in JC_0_c0c0_3d()
231   const PetscScalar dC[]   = {2 * C00, 4 * C01, 4 * C02, 2 * C11, 4 * C12, 2 * C22};  in JC_0_c0c0_3d()
357   const PetscScalar  C02   = u[uOff[C_FIELD_ID] + 2];  in P_1_3d()  local
361   const PetscScalar  C20   = C02;  in P_1_3d()
366   const PetscScalar  s     = fix_c ? FIX_C_3d(C00, C01, C02, C11, C12, C22) : 0.0;  in P_1_3d()
368   f1[0] = -((C00 + s) * gradp[0] + C01 * gradp[1] + C02 * gradp[2]);  in P_1_3d()
394   const PetscScalar  C02   = a[aOff[C_FIELD_ID] + 2];  in P_1_aux_3d()  local
398   const PetscScalar  C20   = C02;  in P_1_aux_3d()
403   const PetscScalar  s     = fix_c ? FIX_C_3d(C00, C01, C02, C11, C12, C22) : 0.0;  in P_1_aux_3d()
405   f1[0] = (C00 + s) * gradp[0] + C01 * gradp[1] + C02 * gradp[2];  in P_1_aux_3d()
432   const PetscScalar C02   = u[uOff[C_FIELD_ID] + 2];  in JP_1_p1p1_3d()  local
436   const PetscScalar C20   = C02;  in JP_1_p1p1_3d()
440   const PetscScalar s     = fix_c ? FIX_C_3d(C00, C01, C02, C11, C12, C22) : 0.0;  in JP_1_p1p1_3d()
444   J[2] = -C02;  in JP_1_p1p1_3d()
474   const PetscScalar C02   = a[aOff[C_FIELD_ID] + 2];  in JP_1_p1p1_aux_3d()  local
478   const PetscScalar C20   = C02;  in JP_1_p1p1_aux_3d()
482   const PetscScalar s     = fix_c ? FIX_C_3d(C00, C01, C02, C11, C12, C22) : 0.0;  in JP_1_p1p1_aux_3d()
486   J[2] = C02;  in JP_1_p1p1_aux_3d()
639     const PetscScalar  C02     = u[uOff[C_FIELD_ID] + 2];  in energy()  local
643     const PetscScalar  C20     = C02;  in energy()
653     const PetscScalar  normC   = NORM2C_3d(C00, C01, C02, C11, C12, C22) + eps;  in energy()
658 …const PetscScalar t1 = gradp[0] * ((C00 + r) * gradp[0] + C01 * gradp[1] + C02 * gradp[2]) + gradp…  in energy()
670   PetscScalar     C00, C01, C02 = 0.0, C11, C12 = 0.0, C22 = 0.0;  in ellipticity_fail_private()  local
680     C02 = u[uOff[C_FIELD_ID] + 2];  in ellipticity_fail_private()
685   Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs);  in ellipticity_fail_private()
796     const PetscScalar C02 = u[uOff[C_FIELD_ID] + 2];  in flux()  local
800     const PetscScalar C20 = C02;  in flux()
804     f[0] = -C00 * gradp[0] - C01 * gradp[1] - C02 * gradp[2];  in flux()
822     const PetscScalar C02 = u[uOff[C_FIELD_ID] + 2];  in normc()  local
827     f[0] = PetscSqrtReal(PetscRealPart(NORM2C_3d(C00, C01, C02, C11, C12, C22)));  in normc()
836   PetscScalar C00, C01, C02 = 0.0, C11, C12 = 0.0, C22 = 0.0;  in eigsc()  local
844     C02 = u[uOff[C_FIELD_ID] + 2];  in eigsc()
849   Eigenvalues_Sym3x3(C00, C01, C02, C11, C12, C22, eigs);  in eigsc()