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