Rheolef  7.2
an efficient C++ finite element environment
vector_projection.h

The projection for yield-stress rheologies – vector-valued case for the Mossolov problem

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