void
PhaseFieldElement :: computeNd_matrixAt(const FloatArray &lCoords, FloatMatrix &N)
{
NLStructuralElement *el = this->giveElement();
FloatArray Nvec;
el->giveInterpolation( )->evalN( Nvec, lCoords, FEIElementGeometryWrapper( el ) );
N.resize(1, Nvec.giveSize());
N.beNMatrixOf(Nvec,1);
}
void
PhaseFieldElement :: computeBd_matrixAt(GaussPoint *gp, FloatMatrix &answer, int li, int ui)
{
// Returns the [numSpaceDim x nDofs] gradient matrix {B_d} of the receiver,
// evaluated at gp.
NLStructuralElement *el = this->giveElement();
FloatMatrix dNdx;
el->giveInterpolation( )->evaldNdx( dNdx, *gp->giveNaturalCoordinates( ), FEIElementGeometryWrapper( el ) );
answer.beTranspositionOf( dNdx );
}
double
PhaseFieldElement :: computeDamageAt(GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN)
{
// d = N_d * a_d
NLStructuralElement *el = this->giveElement();
FloatArray dVec;
computeDamageUnknowns(dVec, valueMode, stepN);
FloatArray Nvec;
el->giveInterpolation()->evalN(Nvec, *gp->giveNaturalCoordinates(), FEIElementGeometryWrapper(el));
return Nvec.dotProduct(dVec);
}