45 const Dune::FieldMatrix<Scalar, dim>& T2,
46 const Dune::FieldVector<Scalar, dim>&
normal)
49 Dune::FieldVector<Scalar, dim> tmp;
50 Dune::FieldVector<Scalar, dim> tmp2;
53 const Scalar alpha1 = tmp*
normal;
54 const Scalar alpha2 = tmp2*
normal;
57 const Scalar alphaAverage = 0.5*(alpha1 + alpha2);
59 Dune::FieldMatrix<Scalar, dim> T(0.0);
60 for (
int i = 0; i < dim; ++i)
62 for (
int j = 0; j < dim; ++j)
64 T[i][j] += 0.5*(T1[i][j] + T2[i][j]);
66 T[i][j] += alphaHarmonic - alphaAverage;
constexpr Scalar harmonicMean(Scalar x, Scalar y, Scalar wx=1.0, Scalar wy=1.0) noexcept
Calculate the (weighted) harmonic mean of two scalar values.
Definition math.hh:57
Scalar faceTensorAverage(const Scalar T1, const Scalar T2, const Dune::FieldVector< Scalar, dim > &normal)
Average of a discontinuous scalar field at discontinuity interface (for compatibility reasons with th...
Definition facetensoraverage.hh:29