Ejemplo n.º 1
0
        /// <summary>Initializes a new instance of the <see cref="ClosedNewtonCotesFormula"/> class.
        /// </summary>
        /// <param name="rule">The Newton-Cotes type, i.e. the degree of the Newton-Cotes method.</param>
        /// <param name="exitCondition">The exit condition.</param>
        public ClosedNewtonCotesFormula(Rule rule, ExitCondition exitCondition)
            : base(BoundDescriptor.Closed, BoundDescriptor.Closed)
        {
            ExitCondition   = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
            NewtonCotesType = rule;

            m_Name = new IdentifierString(String.Format("Newton-Cotes {0} Integrator", rule));
        }
Ejemplo n.º 2
0
 /// <summary>Initializes a new instance of the <see cref="GaussHermiteIntegrator"/> class.
 /// </summary>
 /// <param name="initialOrder">The initial order of the Gauss-Hermite approach, i.e. the order in the first iteration step.</param>
 /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussHermiteIntegrator(int initialOrder, int orderStepSize, ExitCondition exitCondition)
     : base(BoundDescriptor.Create(BoundEvaluationType.Unbounded, Double.NegativeInfinity), BoundDescriptor.Create(BoundEvaluationType.Unbounded, Double.PositiveInfinity))
 {
     ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
     m_Name        = new IdentifierString("Gauss-Hermite Integrator");
     InitialOrder  = initialOrder;
     OrderStepSize = orderStepSize;
 }
 /// <summary>Initializes a new instance of the <see cref="GaussLegendreIntegrator"/> class.
 /// </summary>
 /// <param name="initialOrder">The initial order of the Gauss-Legendre approach, i.e. the order in the first iteration step.</param>
 /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussLegendreIntegrator(int initialOrder, int orderStepSize, ExitCondition exitCondition)
     : base(BoundDescriptor.Closed, BoundDescriptor.Closed)
 {
     ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
     m_Name        = new IdentifierString("Gauss-Legendre Integrator");
     InitialOrder  = initialOrder;
     OrderStepSize = orderStepSize;
 }
        /// <summary>Gets informations of the current object as a specific <see cref="T:Dodoni.BasicComponents.Containers.InfoOutput"/> instance.
        /// </summary>
        /// <param name="infoOutput">The <see cref="T:Dodoni.BasicComponents.Containers.InfoOutput"/> object which is to be filled with informations concering the current instance.</param>
        /// <param name="categoryName">The name of the category, i.e. all informations will be added to these category.</param>
        public override void FillInfoOutput(InfoOutput infoOutput, string categoryName = "General")
        {
            var infoOutputPackage = infoOutput.AcquirePackage(categoryName);

            infoOutputPackage.Add("Dimension", 1);
            infoOutputPackage.Add("Weight function is 1.0", true);

            LowerBoundDescriptor.FillInfoOutput(infoOutput, categoryName + ".LowerBoundDescriptor");
            UpperBoundDescriptor.FillInfoOutput(infoOutput, categoryName + ".UpperBoundDescriptor");
            ExitCondition.FillInfoOutput(infoOutput, categoryName + ".ExitCondition");
        }
Ejemplo n.º 5
0
        /// <summary>Gets informations of the current object as a specific <see cref="T:Dodoni.BasicComponents.Containers.InfoOutput"/> instance.
        /// </summary>
        /// <param name="infoOutput">The <see cref="T:Dodoni.BasicComponents.Containers.InfoOutput"/> object which is to be filled with informations concering the current instance.</param>
        /// <param name="categoryName">The name of the category, i.e. all informations will be added to these category.</param>
        public override void FillInfoOutput(InfoOutput infoOutput, string categoryName = "General")
        {
            var infoOutputPackage = infoOutput.AcquirePackage(categoryName);

            infoOutputPackage.Add("Dimension", 1);
            infoOutputPackage.Add("Weight function is 1.0", false);
            infoOutputPackage.Add("Initial order", InitialOrder);
            infoOutputPackage.Add("Order step size", OrderStepSize);

            LowerBoundDescriptor.FillInfoOutput(infoOutput, categoryName + ".LowerBoundDescriptor");
            UpperBoundDescriptor.FillInfoOutput(infoOutput, categoryName + ".UpperBoundDescriptor");
            ExitCondition.FillInfoOutput(infoOutput, categoryName + ".ExitCondition");
        }
        /// <summary>Initializes a new instance of the <see cref="GaussTschebyscheffIntegrator"/> class.
        /// </summary>
        /// <param name="initialOrder">The initial order of the Gauss-Tschebyscheff approach, i.e. the order in the first iteration step.</param>
        /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
        /// <param name="exitCondition">The exit condition.</param>
        public GaussTschebyscheffIntegrator(int initialOrder, int orderStepSize, ExitCondition exitCondition)
            : base(BoundDescriptor.Open, BoundDescriptor.Open)
        {
            ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));

            if (initialOrder < 2)
            {
                throw new ArgumentOutOfRangeException(String.Format(ExceptionMessages.ArgumentOutOfRangeGreaterEqual, "Initial order", 2));
            }
            InitialOrder = initialOrder;

            if (orderStepSize < 1)
            {
                throw new ArgumentOutOfRangeException(String.Format(ExceptionMessages.ArgumentOutOfRangeGreaterEqual, "Order step size", 1));
            }
            OrderStepSize = orderStepSize;

            m_Name = new IdentifierString("Gauss-Tschebyscheff Integrator");
        }
Ejemplo n.º 7
0
 /// <summary>Initializes a new instance of the <see cref="GaussKronrodPatterson255ConstAbscissaIntegrator"/> class.
 /// </summary>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussKronrodPatterson255ConstAbscissaIntegrator(ExitCondition exitCondition)
     : base(OneDimNumericalIntegrator.BoundDescriptor.Closed, OneDimNumericalIntegrator.BoundDescriptor.Closed)
 {
     ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
     m_Name        = new IdentifierString("Gauss-Kronrod-Patterson 255 const abscissa Integrator");
 }
Ejemplo n.º 8
0
 /// <summary>Initializes a new instance of the <see cref="RombergIntegrator"/> class.
 /// </summary>
 /// <param name="exitCondition">The exit condition.</param>
 public RombergIntegrator(ExitCondition exitCondition)
     : base(BoundDescriptor.Closed, BoundDescriptor.Closed)
 {
     ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
     m_Name        = new IdentifierString("Romberg Integrator");
 }
 /// <summary>Initializes a new instance of the <see cref="GaussKronrodPatterson255Integrator"/> class.
 /// </summary>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussKronrodPatterson255Integrator(ExitCondition exitCondition)
     : base(BoundDescriptor.Closed, BoundDescriptor.Closed)
 {
     ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));
     m_Name        = new IdentifierString("Gauss-Kronrod-Patterson 255 Integrator");
 }
Ejemplo n.º 10
0
 /// <summary>Initializes a new instance of the <see cref="GaussLaguerreIntegrator"/> class.
 /// </summary>
 /// <param name="alpha">The parameter \alpha of the Guass-Laguerre integrator.</param>
 /// <param name="initialOrder">The initial order of the Gauss-Laguerre approach, i.e. the order in the first iteration step.</param>
 /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussLaguerreIntegrator(double alpha, int initialOrder, int orderStepSize, ExitCondition exitCondition)
     : this(alpha, initialOrder, orderStepSize, exitCondition, alphaIsZero : false)
 {
     m_Weights   = new Dictionary <int, double[]>();
     m_Abscissas = new Dictionary <int, double[]>();
 }
Ejemplo n.º 11
0
        /// <summary>Initializes a new instance of the <see cref="GaussLaguerreIntegrator"/> class.
        /// </summary>
        /// <param name="alpha">The parameter \alpha of the Guass-Laguerre integrator.</param>
        /// <param name="initialOrder">The initial order of the Gauss-Laguerre approach, i.e. the order in the first iteration step.</param>
        /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
        /// <param name="exitCondition">The exit condition.</param>
        /// <param name="alphaIsZero">A value indicating whether the parameter <paramref name="alpha"/> is 0.0.</param>
        private GaussLaguerreIntegrator(double alpha, int initialOrder, int orderStepSize, ExitCondition exitCondition, bool alphaIsZero)
            : base(OneDimNumericalIntegrator.BoundDescriptor.Create(OneDimNumericalIntegrator.BoundEvaluationType.Closed, 0.0), OneDimNumericalIntegrator.BoundDescriptor.Create(OneDimNumericalIntegrator.BoundEvaluationType.Unbounded, Double.PositiveInfinity))
        {
            ExitCondition = exitCondition ?? throw new ArgumentNullException(nameof(exitCondition));

            if (Double.IsNaN(alpha) == true)
            {
                throw new ArgumentOutOfRangeException("alpha");
            }
            Alpha         = alpha;
            m_Name        = new IdentifierString(String.Format("Gauss-Laguerre Integrator; alpha: {0}", alpha));
            InitialOrder  = initialOrder;
            OrderStepSize = orderStepSize;
            m_AlphaIsZero = alphaIsZero;
        }
Ejemplo n.º 12
0
 /// <summary>Initializes a new instance of the <see cref="GaussLaguerreIntegrator"/> class with \alpha = 0.0.
 /// </summary>
 /// <param name="initialOrder">The initial order of the Gauss-Laguerre approach, i.e. the order in the first iteration step.</param>
 /// <param name="orderStepSize">The step size of the order, i.e. in each iteration step the order will be increased by the specified number.</param>
 /// <param name="exitCondition">The exit condition.</param>
 public GaussLaguerreIntegrator(int initialOrder, int orderStepSize, ExitCondition exitCondition)
     : this(0.0, initialOrder, orderStepSize, exitCondition, alphaIsZero : true)
 {
 }