The projection for yield-stress rheology – its derivative.
The projection for yield-stress rheology – its derivative
if (fabs(x) <=
a)
return 0;
if (
n == 1)
return 1/(
c +
r);
if (
r == 0)
return pow(fabs(x)-
a,-1+1/
n)/(
n*pow(
c,1/
n));
return 1/(
r +
n*
c*pow(
_phi(fabs(x)-
a),-1+
n));
}
:
a(a1),
n(n1),
c(c1),
r(r1),
_phi(n1,c1,r1) {}
};
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