Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
mosolov_exact_circle.h
Go to the documentation of this file.
1
25struct u {
26 Float operator() (const point& x) const {
27 return (pow(1-Bi,1+1/n) - pow(max(Float(0),norm(x)-Bi),1+1/n))/(1+1/n);
28 }
29 u (Float Bi1, Float n1) : Bi(Bi1), n(n1) {}
30 protected: Float Bi, n;
31};
32struct grad_u {
33 point operator() (const point& x) const {
34 Float r = norm(x);
35 return (r <= Bi) ? point(0,0) : -pow(r-Bi, 1/n)*(x/r);
36 }
37 grad_u (Float Bi1, Float n1) : Bi(Bi1), n(n1) {}
38 protected: Float Bi, n;
39};
40struct sigma {
41 point operator() (const point& x) const { return -x; }
42 sigma (Float=0, Float=0) {}
43};
see the Float page for the full documentation
see the point page for the full documentation
point operator()(const point &x) const
grad_u(Float Bi1, Float n1)
sigma(Float=0, Float=0)
point operator()(const point &x) const
Definition leveque.h:25
point operator()(const point &x) const
Definition leveque.h:26
u(Float Bi1, Float n1)