Proto_Operator.H

#pragma once

#include "Proto_BoxData.H"
//#include "AMRState.H"

namespace Operator {

    using namespace Proto;    

    /*
    template<unsigned int C>
    inline void faceDivergence(PatchState<C, 0>& a_div,
            const PatchState<0,C>& a_flux,
            double a_dx);
    */

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> convolve(
            const BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void convolve(
            BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_2nd);


    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> deconvolve(
            const BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void deconvolve(
            BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_2nd);


    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> convolveFace(
            const BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_2nd,
            int a_dir);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void convolveFace(
            BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_2nd,
            int a_dir);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> deconvolveFace(
            const BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_2nd,
            int a_dir);
   

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void deconvolveFace(
            BoxData<T,C,D,E>& a_ctr,
            const BoxData<T,C,D,E>& a_avg,
            const BoxData<T,C,D,E>& a_2nd,
            int a_dir);
   

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> cellProduct(
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void cellProduct(
            BoxData<T,C,D,E>& a_product,
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> faceProduct(
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd,
            int a_dir);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void faceProduct(BoxData<T,C,D,E>& a_dest,
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd,
            int a_dir);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> cellQuotient(
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void cellQuotient(
            BoxData<T,C,D,E>& a_quotient,
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd);
    

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> faceQuotient(
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd,
            int a_dir);


    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void faceQuotient(
            BoxData<T,C,D,E>& a_quotient,
            const BoxData<T,C,D,E>& a_pavg,
            const BoxData<T,C,D,E>& a_qavg,
            const BoxData<T,C,D,E>& a_p2nd,
            const BoxData<T,C,D,E>& a_q2nd,
            int a_dir);

    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> cellToFace(
            const BoxData<T,C,D,E>& a_cell,
            int a_dir);
    
    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void cellToFace(
            BoxData<T,C,D,E>& a_face,
            const BoxData<T,C,D,E>& a_cell,
            int a_dir);
    
    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline BoxData<T,C,D,E> diffCellToFace(
            const BoxData<T,C,D,E>& a_data,
            int a_dir);
    
    template<class T, unsigned int C, unsigned char D, unsigned char E>
    inline void diffCellToFace(
            BoxData<T,C,D,E>& a_deriv,
            const BoxData<T,C,D,E>& a_data,
            int a_dir);
}
#include "implem/Proto_OperatorI.H"