Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
projection.h

The projection for yield-stress rheologies e.g. the yield slip problem.

The projection for yield-stress rheologies e.g. the yield slip problem

#include "phi.h"
// p(x) = phi(|x|-a)*sgn(x)
struct projection {
Float operator() (const Float& x) const {
if (fabs(x) <= a) return 0;
return (x > 0) ? _phi(x-a) : -_phi(-x-a);
}
projection (Float a1, Float n=1, Float c=1, Float r=0)
: a(a1), _phi(n,c,r) {}
};
see the Float page for the full documentation
The Mossolov problem – the phi function.
Definition phi.h:25
Float operator()(const Float &x) const
Definition projection.h:28