void elasticitySolver::setLagrangeMultipliers(int phys, double tau,
const SVector3 &d, int tag,
simpleFunction<double> *f)
{
LagrangeMultiplierField field;
field._tau = tau;
field._tag = tag;
field._f = f;
field._d = d.unit();
field.g = new groupOfElements(
_dim - 1, phys); // LM applied on entity of dimension = dim-1!
LagrangeMultiplierFields.push_back(field);
LagrangeMultiplierSpaces.push_back(
new ScalarLagrangeFunctionSpaceOfElement(tag));
}
