Rheolef  7.1
an efficient C++ finite element environment
eta.h
Go to the documentation of this file.
1 struct eta {
26  Float operator() (const Float& z) const {
27  check_macro(z != 0 || p > 2, "eta: division by zero (HINT: check mesh)");
28  return pow(z, (p-2)/2);
29  }
30  Float derivative (const Float& z) const {
31  check_macro(z != 0 || p > 4, "eta': division by zero (HINT: check mesh)");
32  return 0.5*(p-2)*pow(z, (p-4)/2);
33  }
34  eta (const Float& q) : p(q) {}
36 };
check_macro
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
eta::p
Float p
Definition: eta.h:35
eta
Definition: eta.h:25
rheolef::pow
space_mult_list< T, M > pow(const space_basic< T, M > &X, size_t n)
Definition: space_mult.h:120
p
Definition: sphere.icc:25
eta::eta
eta(const Float &q)
Definition: eta.h:34
Float
see the Float page for the full documentation
eta::derivative
Float derivative(const Float &z) const
Definition: eta.h:30
eta::operator()
Float operator()(const Float &z) const
Definition: eta.h:26