prescribedgradientbcweakdirichlet.C

Go to the documentation of this file.

/*
 * prescribedgradientbcweakdirichlet.C
 *
 *  Created on: May 19, 2014
 *      Author: svennine
 */
 
#include "prescribedgradientbcweakdirichlet.h"
 
#include "classfactory.h"
 
namespace oofem {
REGISTER_BoundaryCondition(PrescribedGradientBCWeakDirichlet);
 
PrescribedGradientBCWeakDirichlet :: PrescribedGradientBCWeakDirichlet(int n, Domain *d) :
    PrescribedGradientBCWeak(n, d)
{
    // TODO Auto-generated constructor stub
}
 
PrescribedGradientBCWeakDirichlet :: ~PrescribedGradientBCWeakDirichlet()
{
    // TODO Auto-generated destructor stub
}
 
void PrescribedGradientBCWeakDirichlet :: initializeFrom(InputRecord &ir)
{
    PrescribedGradientBCWeak :: initializeFrom(ir);
    mMeshIsPeriodic = false;
}
 
void PrescribedGradientBCWeakDirichlet :: postInitialize()
{
    bool enforceCornerPeriodicity = false;
    int numSides = 4;
    createTractionMesh(enforceCornerPeriodicity, numSides);
}
 
void PrescribedGradientBCWeakDirichlet :: giveBoundaryCoordVector(FloatArray &oX, const FloatArray &iPos) const
{
    oX = Vec2(
        iPos [ 0 ] - mCenterCoord [ 0 ], iPos [ 1 ] - mCenterCoord [ 1 ]
    );
}
 
void PrescribedGradientBCWeakDirichlet :: checkIfCorner(bool &oIsCorner, bool &oDuplicatable, const FloatArray &iPos, const double &iNodeDistTol) const
{
    oIsCorner = false;
    oDuplicatable = false;
 
    FloatArray cornerPos = mLC;
    if ( distance(iPos, cornerPos) < iNodeDistTol ) {
        oIsCorner = true;
        oDuplicatable = true;
    }
 
    cornerPos = Vec2(
        mUC [ 0 ], mLC [ 1 ]
    );
    if ( distance(iPos, cornerPos) < iNodeDistTol ) {
        oIsCorner = true;
        oDuplicatable = true;
    }
 
    cornerPos = Vec2(
        mUC [ 0 ], mUC [ 1 ]
    );
    if ( distance(iPos, cornerPos) < iNodeDistTol ) {
        oIsCorner = true;
        oDuplicatable = true;
    }
 
    cornerPos = Vec2(
        mLC [ 0 ], mUC [ 1 ]
    );
    if ( distance(iPos, cornerPos) < iNodeDistTol ) {
        oIsCorner = true;
        oDuplicatable = true;
    }
}
} /* namespace oofem */