core/hull.hh

#ifndef VTESTBED_CORE_HULL_HH
#define VTESTBED_CORE_HULL_HH

#include <vtestbed/base/constants.hh>
#include <vtestbed/geometry/polyhedron.hh>

namespace vtb {

    namespace core {

        template <class T>
        class Hull: public vtb::geometry::Polyhedron<T,3> {

        private:
            using base_type = vtb::geometry::Polyhedron<T,3>;

        public:
            using typename base_type::scalar_type;
            using typename base_type::vertex_type;
            using typename base_type::face_type;
            using typename base_type::vertex_array;
            using typename base_type::face_array;
            using typename base_type::size_type;
            using typename base_type::index_type;
            using typename base_type::box_type;
            using typename base_type::triangle_type;
            using typename base_type::triangle_array;
            using typename base_type::bsp_tree;

        private:
            scalar_type _mass = 0;

        public:

            using base_type::base_type;
            using base_type::operator=;

            inline Hull(base_type&& rhs): base_type(std::move(rhs)) {}
            inline Hull(const base_type& rhs): base_type(rhs) {}

            inline scalar_type mass() const { return this->_mass; }

            inline void mass(scalar_type rhs) {
                if (rhs < 0) { throw std::invalid_argument("bad mass"); }
                this->_mass = rhs;
            }

            inline void
            displacement(scalar_type rhs) {
                this->mass(vtb::base::constants<T>::water_density()*rhs);
            }

            inline scalar_type
            displacement() const {
                return this->mass() / vtb::base::constants<T>::water_density();
            }

            inline void
            draught(scalar_type rhs) {
                if (rhs < 0) { throw std::invalid_argument("bad draught"); }
                auto level = this->bounds(2).min() + rhs;
                this->displacement(this->signed_volume_below(2, level));
            }

        };

    }

}

#endif // vim:filetype=cpp