System.hpp

/*
 * System.hpp
 *
 *  Created on: Oct 5, 2015
 *      Author: i-bird
 */

#ifndef OPENFPM_NUMERICS_SRC_FINITEDIFFERENCE_SYSTEM_HPP_
#define OPENFPM_NUMERICS_SRC_FINITEDIFFERENCE_SYSTEM_HPP_


template<unsigned int dim, unsigned int nvf, unsigned int ncf, typename ... eqs>
class System
{
    // Define the number of constant fields
    typedef num_cfields boost::mpl::int_<nf>;

    // Define the number of variable fields
    typedef num_vfields boost::mpl::int_<nf>;

    // set of equations as boost::mpl::vector
    typedef eqs_v make_vactor<eqs>;

    template<unsigned int ord=EQS_FIELDS> void value(const grid_key_dx_dist<dim> & pos)
    {
        if (EQS_FIELDS)
            value_f(pos);
        else
            value_s(pos);
    }

    template<unsigned int eq_id> void value_s(grid_key_dx_dist<dim> & it)
    {
        boost::mpl::at<eqs_v,boost::mpl::int_<eq_id>>::type eq;

        eq.value_s(it);
    }

    template<unsigned int eq_id> void value_f(grid_key_dx_dist<dim> & it)
    {
        boost::mpl::at<eqs_v,boost::mpl::int_<eq_id>>::type eq;

        eq.value_f(it);
    }
};


#endif /* OPENFPM_NUMERICS_SRC_FINITEDIFFERENCE_SYSTEM_HPP_ */