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 * FEI2dQuadLin :: giveIntegrationRule(int order) { IntegrationRule *iRule = new GaussIntegrationRule(1, NULL); int points = iRule->getRequiredNumberOfIntegrationPoints(_Square, order + 2); iRule->SetUpPointsOnSquare(points, _Unknown); return iRule; }
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; }