ship_hull_panel.hh

#ifndef VTESTBED_CORE_SHIP_HULL_PANEL_HH
#define VTESTBED_CORE_SHIP_HULL_PANEL_HH

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

namespace vtb {

    namespace core {

        template <class T, int N>
        class alignas(16) Ship_hull_panel {

        public:
            static constexpr const int
                Centre = 0,
                Normal = 1,
                Force = 2,
                Polyhedron = 3;

        public:
            using vec = blitz::TinyVector<T,N>;

        private:
            vec _centre{T{}};
            vec _normal{T{}};
            vec _force{T{}};
            T _area{T{}};
            T _volume{T{}};
            union {
                struct { bool _wetted = false; bool _waterline = false; };
                T padding[1];
            };

        public:

            Ship_hull_panel(const vec& centre, vec n, T area, T volume);
            Ship_hull_panel() = default;
            ~Ship_hull_panel() = default;
            Ship_hull_panel& operator=(const Ship_hull_panel&) = default;

            inline const vec& centre() const { return this->_centre; }
            inline void centre(const vec& rhs) { this->_centre = rhs; }

            inline const vec& normal() const { return this->_normal; }
            inline void normal(const vec& rhs) { this->_normal = rhs; }

            inline T area() const { return this->_area; }
            inline void area(T rhs) { this->_area = rhs; }

            inline T volume() const { return this->_volume; }
            inline void volume(T rhs) { this->_volume = rhs; }

            inline const vec& force() const { return this->_force; }
            inline vec& force() { return this->_force; }

            inline bool wetted() const { return this->_wetted; }
            inline void wetted(bool rhs) { this->_wetted = rhs; }

            inline bool waterline() const { return this->_waterline; }
            inline void waterline(bool rhs) { this->_waterline = rhs; }

        };

        template <class T>
        Grid<T,3>
        clamp(const Grid<T,3>& grid_xyz, const std::vector<Ship_hull_panel<T,3>>& panels);

    }

}

#endif // vim:filetype=cpp