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