Rheolef  7.1
an efficient C++ finite element environment
projection.h

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) {}
};
projection::a
Float a
Definition: projection.h:34
phi
Definition: phi.h:25
mkgeo_ball.c
c
Definition: mkgeo_ball.sh:153
phi.h
The Mossolov problem – the phi function.
projection::projection
projection(Float a1, Float n=1, Float c=1, Float r=0)
Definition: projection.h:32
a
Definition: diffusion_isotropic.h:25
projection::_phi
phi _phi
Definition: projection.h:35
projection::operator()
Float operator()(const Float &x) const
Definition: projection.h:28
Float
see the Float page for the full documentation
mkgeo_ball.n
n
Definition: mkgeo_ball.sh:150
projection
Definition: projection.h:27