using System;
|
|
namespace IStation.Numerics.Integration.GaussRule
|
{
|
/// <summary>
|
/// Creates a Gauss-Kronrod point.
|
/// </summary>
|
internal static class GaussKronrodPointFactory
|
{
|
[ThreadStatic]
|
private static GaussPointPair gaussKronrodPoint;
|
|
/// <summary>
|
/// Getter for the GaussKronrodPoint.
|
/// </summary>
|
/// <param name="order">Defines an Nth order Gauss-Kronrod rule. Precomputed Gauss-Kronrod abscissas/weights for orders 15, 21, 31, 41, 51, 61 are used, otherwise they're calculated on the fly.</param>
|
/// <returns>Object containing the non-negative abscissas/weights, and order.</returns>
|
public static GaussPointPair GetGaussPoint(int order)
|
{
|
// Try to get the GaussKronrodPoint from the cached static field.
|
bool gaussKronrodPointIsCached = gaussKronrodPoint != null && gaussKronrodPoint.Order == order;
|
if (!gaussKronrodPointIsCached)
|
{
|
// Try to find the GaussKronrodPoint in the precomputed dictionary.
|
if (!GaussKronrodPoint.PreComputed.TryGetValue(order, out gaussKronrodPoint))
|
{
|
gaussKronrodPoint = GaussKronrodPoint.Generate(order, 1E-10);
|
}
|
}
|
|
return gaussKronrodPoint;
|
}
|
}
|
}
|
