/// <summary>
/// Gauss rules for \f$ \oint_{\frakI \cap K_j } \ldots \dS \f$ in the 3D case
/// </summary>
///
public static IQuadRuleFactory <QuadRule> SayeGaussRule_Volume3D(
LevelSetTracker.LevelSetData _lsData,
IRootFindingAlgorithm RootFinder)
{
ISayeGaussRule rule = new SayeFactory_Cube(
_lsData,
RootFinder,
SayeFactory_Cube.QuadratureMode.Volume);
return(new SayeGaussRuleFactory(rule));
}
public static SayeGaussComboRuleFactory SayeGaussRule_Combo3D(
LevelSetTracker.LevelSetData _lsData,
IRootFindingAlgorithm RootFinder
)
{
ISayeGaussComboRule rule = new SayeFactory_Cube(
_lsData,
RootFinder,
SayeFactory_Cube.QuadratureMode.Combo
);
return(new SayeGaussComboRuleFactory(
rule,
SayeGaussComboRuleFactory.Mode.CalculateOnceAfterInstantiation));
}
public static SayeGaussComboRuleFactory SayeGaussRule_Combo3D(
LevelSetTracker.LevelSetData lsData,
IRootFindingAlgorithm RootFinder
)
{
ISayeGaussComboRule rule = new SayeFactory_Cube(
lsData,
RootFinder,
SayeFactory_Cube.QuadratureMode.Combo
);
CellMask maxGrid = lsData.GridDat.Cells.GetCells4Refelement(rule.RefElement).Intersect(
lsData.Region.GetCutCellMask().ToGeometicalMask());
return(new SayeGaussComboRuleFactory(
rule,
maxGrid));
}