opencl/derivative.hh

#ifndef VTESTBED_OPENCL_DERIVATIVE_HH
#define VTESTBED_OPENCL_DERIVATIVE_HH

#include <fstream>
#include <sstream>

#include <vtestbed/base/blitz.hh>
#include <vtestbed/core/grid.hh>
#include <vtestbed/core/types.hh>
#include <vtestbed/opencl/opencl.hh>
#include <vtestbed/opencl/pipeline.hh>

namespace vtb {

    namespace opencl {

        template <class T, int N>
        class Derivative: public Context_base {

        public:
            using grid_type = vtb::core::Grid<T,N>;
            using buffer_type = Buffer<T>;

        private:
            clx::kernel _derivative, _local_derivative;

        public:
            using Context_base::context;

            void context(Context* rhs) override;

            void solve(const buffer_type& func,
                       const grid_type& grid,
                       int dimension,
                       buffer_type& result);

            inline void operator()(
                const buffer_type& func,
                const grid_type& grid,
                int dimension,
                buffer_type& result
            ) { this->solve(func, grid, dimension, result); }

        };

    }

}

#endif // vim:filetype=cpp