Rheolef  7.2
an efficient C++ finite element environment
yield_slip_circle.h
Go to the documentation of this file.
1 struct u {
26  Float operator() (const point& x) const { return (1-norm2(x))/4 + us; }
27  u (Float S, Float n, Float Cf) : us(pow(max(Float(0),(0.5-S)/Cf), 1/n)) {}
28  protected: Float us;
29 };
30 struct grad_u {
31  point operator() (const point& x) const { return -x/2; }
32  grad_u (Float S, Float n, Float Cf) {}
33 };
34 struct lambda {
35  Float operator() (const point& x) const { return 1./2; }
36  lambda (Float S, Float n, Float Cf) {}
37 };
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 S, Float n, Float Cf)
lambda(Float S, Float n, Float Cf)
Float operator()(const point &x) const
Definition: leveque.h:25
Float us
point operator()(const point &x) const
Definition: leveque.h:26
u(Float S, Float n, Float Cf)