Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
d_projection_dx.h
Go to the documentation of this file.
1
25#include "projection.h"
27 Float operator() (const Float& x) const {
28 if (fabs(x) <= a) return 0;
29 if (n == 1) return 1/(c + r);
30 if (r == 0) return pow(fabs(x)-a,-1+1/n)/(n*pow(c,1/n));
31 return 1/(r + n*c*pow(_phi(fabs(x)-a),-1+n));
32 }
33 d_projection_dx (Float a1, Float n1=1, Float c1=1, Float r1=0)
34 : a(a1), n(n1), c(c1), r(r1), _phi(n1,c1,r1) {}
37};
see the Float page for the full documentation
The projection for yield-stress rheologies e.g. the yield slip problem.
Float operator()(const Float &x) const
d_projection_dx(Float a1, Float n1=1, Float c1=1, Float r1=0)
Definition phi.h:25