Rheolef  7.2
an efficient C++ finite element environment
p_laplacian_circle.h
Go to the documentation of this file.
1 struct u_exact {
26  Float operator() (const point& x) const {
27  return (1 - pow(norm2(x), p/(2*p-2)))/((p/(p-1))*pow(2.,1/(p-1)));
28  }
29  u_exact (Float q) : p(q) {}
30  protected: Float p;
31 };
32 struct grad_u {
33  point operator() (const point& x) const {
34  return - (pow(norm2(x), p/(2*p-2) - 1)/pow(2.,1/(p-1)))*x;
35  }
36  grad_u (Float q) : p(q) {}
37  protected: Float p;
38 };
see the Float page for the full documentation
see the point page for the full documentation
T norm2(const vec< T, M > &x)
norm2(x): see the expression page for the full documentation
Definition: vec.h:379
space_mult_list< T, M > pow(const space_basic< T, M > &X, size_t n)
Definition: space_mult.h:120
point operator()(const point &x) const
grad_u(Float q)
Definition: sphere.icc:25
point operator()(const point &x) const
u_exact(Float q)