Muscat.FE.IntegrationRules module

CheckIntegrity(GUI: bool = False)

class ElementaryQuadrature(elementType: ElementType, points: ArrayLike, weights: ArrayLike)

Bases: object

CopyFor(elementType: ElementType) → ElementaryQuadrature

GetNumberOfPoints()

elementType: ElementType

points: ndarray

weights: ndarray

GenerateTrapezoidalOrderRuleOfOrder(i)

GetConstructorTrapezoidalOrderRuleOfOrder(i)

GetIntegrationRule(ruleName: str | None = None, rule: MeshQuadrature | None = None) → MeshQuadrature

Return the integration rule for the mesh using the name or an integration rule. In the case where both are None then LagrangeIsoParamQuadrature is returned

Parameters:

• ruleName (str, optional) – name of the integration rule, by default None

• rule (MeshQuadrature, optional) – integration rule in the case ruleName is None, by default None

Returns:

the integration rule

Return type:

MeshQuadrature

GetIntegrationRuleByName(ruleName: str) → MeshQuadrature

class IntegrationRulesFactory

Bases: Factory

class MeshQuadrature(name: str)

Bases: object

AddElementaryQuadrature(elementaryQuadrature: ElementaryQuadrature)

elementQuadrature: Dict[ElementType, ElementaryQuadrature]

items() → Generator[ED.ElementType, ElementaryQuadrature]

keys() → ElementType

name: str

values() → Generator[ElementaryQuadrature]

RegisterIntegrationRule(ir: MeshQuadrature)

TensorProd(p1: ArrayLike, w1: ArrayLike, p2: ArrayLike, w2: ArrayLike, p3: ArrayLike | None = None, w3: ArrayLike | None = None)

TensorProdHomogeneous(dim: int, p: ArrayLike, w: ArrayLike)

TensorProductGauss(dim: int, nbPoints: int = 2) → Tuple[np.ndarray, np.ndarray]

Tensor product of gauss point in dimension (dim) using nbPoints (number of points) per dimension

Parameters:

• dim (int) – the dimensionality of the tensor product

• nbPoints (int, optional) – Number of points per dimension , by default 2

Returns:

Integration Rule (points and weights)

Return type:

IntegrationRuleType