IntegrationRule *
FEI3dWedgeLin :: giveBoundaryIntegrationRule(int order, int boundary)
{
IntegrationRule *iRule = new GaussIntegrationRule(1, NULL);
if ( boundary <= 2 ) {
int points = iRule->getRequiredNumberOfIntegrationPoints(_Triangle, order + 0);
iRule->SetUpPointsOnTriangle(points, _Unknown);
} else {
int points = iRule->getRequiredNumberOfIntegrationPoints(_Square, order + 2);
iRule->SetUpPointsOnSquare(points, _Unknown);
}
return iRule;
}
IntegrationRule *FEI2dLineLin :: giveIntegrationRule(int order)
{
IntegrationRule *iRule = new GaussIntegrationRule(1, NULL);
int points = iRule->getRequiredNumberOfIntegrationPoints(_Line, order + 0);
iRule->SetUpPointsOnLine(points, _Unknown);
return iRule;
}
IntegrationRule *
Tria1PlateSubSoil :: GetSurfaceIntegrationRule(int approxOrder)
{
IntegrationRule *iRule = new GaussIntegrationRule(1, this, 1, 1);
int npoints = iRule->getRequiredNumberOfIntegrationPoints(_Triangle, approxOrder);
iRule->SetUpPointsOnTriangle(npoints, _Unknown);
return iRule;
}
IntegrationRule *
QBrick1_ht :: GetSurfaceIntegrationRule(int approxOrder)
{
IntegrationRule *iRule = new GaussIntegrationRule(1, this, 1, 1);
int npoints = iRule->getRequiredNumberOfIntegrationPoints(_Square, approxOrder);
iRule->SetUpPointsOnSquare(npoints, _Unknown);
return iRule;
}
IntegrationRule *
FEI3dTrQuad :: giveIntegrationRule(int order)
{
IntegrationRule *iRule = new GaussIntegrationRule(1, NULL);
int points = iRule->getRequiredNumberOfIntegrationPoints(_Triangle, order);
iRule->SetUpPointsOnTriangle(points, _Unknown);
return iRule;
}