using System;

namespace Science.Mathematics.VectorCalculus
{
    public class IntegralOf1Form
    {
        private Function.ToLastType<double, double, double, double> Pdx, Qdy, Rdz;
        private Path p;
        public IntegralOf1Form(Function.ToLastType<double, double, double, double> dx,
            Function.ToLastType<double, double, double, double> dy,
            Function.ToLastType<double, double, double, double> dz, Path c)
        {
            Pdx = dx;
            Qdy = dy;
            Rdz = dz;
            p = c;
        }
        public void Compute()
        {
            Function.ToLastType<double[],double[]> f
                = new Function.ToLastType<double[],double[]>(func);
            LineIntegral obj = new LineIntegral(f, p);
            obj.Compute();
            res = obj.Result;
        }
        private double[] func(double[] x)
        {
            double[] r = new double[3];
            r[0] = Pdx(x[0],x[1],x[2]);
            r[1] = Qdy(x[0],x[1],x[2]);
            r[2] = Rdz(x[0],x[1],x[2]);
            return r;
        }
        private double res;
        public double Result
        {
            get { return res; }
        }
    }    
}