using System;

namespace Science.Mathematics.VectorCalculus
{
    public class IntegralOf2Form
    {
        private Function.ToLastType<double, double, double, double> Fdxdy, Gdydz, Hdzdx;
        private Surface d;
        public IntegralOf2Form(Function.ToLastType<double, double, double, double> dxdy,
            Function.ToLastType<double, double, double, double> dydz,
            Function.ToLastType<double, double, double, double> dzdx, Surface s)
        {
            Fdxdy = dxdy;
            Gdydz = dydz;
            Hdzdx = dzdx;
            d = s;
        }
        public void Compute()
        {
            Function.ToLastType<double[],double[]> f
                = new Function.ToLastType<double[],double[]>(func);
            SurfaceIntegral obj = new SurfaceIntegral(f, d);
            obj.Compute();
            res = obj.Result;
        }
        private double[] func(double[] x)
        {
            double[] r = new double[3];
            r[0] = Gdydz(x[0],x[1],x[2]);
            r[1] = Hdzdx(x[0],x[1],x[2]);
            r[2] = Fdxdy(x[0],x[1],x[2]);
            return r;
        }
        private double res;
        public double Result
        {
            get { return res; }
        }
    }    
}