sequential-integration 0.0.2

use core::marker::PhantomData;
use fehler::throws;
use mexprp::{Context, Expression};

use crate::{
    engine::{
        helper_equation_traits::{Bounds, EquationOfOneVariable, EquationOfTwoVariable},
        quadrature::GetQuadratureRange,
        utils, CalculationResult, CalculationStep,
    },
    errors::Error,
};

pub struct SecondIntegrator<G: GetQuadratureRange, E: EquationOfTwoVariable> {
    a_equation: Expression<f64>,
    b_equation: Expression<f64>,
    h: f64,
    equation: E,
    _p: PhantomData<G>,
}

impl<G: GetQuadratureRange, E: EquationOfTwoVariable> SecondIntegrator<G, E> {
    #[throws]
    pub fn new(a_equation: &str, b_equation: &str, h: f64, equation: E) -> Self {
        let a_equation = Expression::parse(a_equation)?;
        let b_equation = Expression::parse(b_equation)?;

        Self {
            a_equation,
            b_equation,
            h,
            equation,
            _p: PhantomData,
        }
    }
}

impl<G: GetQuadratureRange, E: EquationOfTwoVariable> EquationOfOneVariable
    for SecondIntegrator<G, E>
{
    #[throws]
    fn calculate(&self, x: CalculationStep, bounds: Bounds) -> CalculationResult {
        let mut context = Context::new();
        context.set_var("x", *x);

        let a = utils::calculate_expression_one_value_result(&context, &self.a_equation)?;
        let b = utils::calculate_expression_one_value_result(&context, &self.b_equation)?;

        let mut result = CalculationResult::new();
        let mut range = if let Some(range) = G::get_range_generator(a, b, self.h)? {
            range
        } else {
            return result;
        };

        loop {
            let step = range.next()?;
            result += self.equation.calculate(x, bounds, step, (a, b))?;

            if step.is_last() {
                break;
            }
        }

        result
    }
}