wave_pressure.hh

#ifndef VTESTBED_CORE_WAVE_PRESSURE_HH
#define VTESTBED_CORE_WAVE_PRESSURE_HH

#include <vtestbed/core/grid.hh>
#include <vtestbed/core/math.hh>
#include <vtestbed/core/policy.hh>
#include <vtestbed/core/ship_hull_panel.hh>
#include <vtestbed/core/types.hh>

namespace vtb {

    namespace core {

        template <class T>
        class Wave_pressure_solver {

        public:
            using panel_type = Ship_hull_panel<T,3>;
            using panel_array = std::vector<panel_type>;
            using grid3 = Grid<T,3>;
            using grid4 = Grid<T,4>;
            using array3 = Array<T,3>;
            using array4 = Array<T,4>;
            using vec3 = Vector3<T>;

        private:
            T _density{1030};
            T _g{base::constants<T>::g()};
            bool _clip = true;

        public:

            virtual ~Wave_pressure_solver() = default;

            virtual void
            solve(
                const Grid<T,3>& wavy_surface_grid,
                Array<T,3> wavy_surface,
                const Grid<T,4>& velocity_potential_grid,
                Array<T,4> velocity_potential,
                panel_array& panels
            ) = 0;

            inline T density() const noexcept { return this->_density; }
            inline void density(T rhs) noexcept { this->_density = rhs; }
            inline T g() const noexcept { return this->_g; }
            inline void g(T rhs) noexcept { this->_g = rhs; }
            inline void clip(bool b) { this->_clip = b; }
            inline bool clip() const { return this->_clip; }

        };

        template <class T, Policy P>
        std::unique_ptr<Wave_pressure_solver<T>>
        make_wave_pressure_solver();

    }

}

#endif // vim:filetype=cpp