levelSetOctreeKernel.hpp

/*---------------------------------------------------------------------------*\
 *
 *  bitpit
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of bitpit.
 *
 *  bitpit is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License v3 (LGPL)
 *  as published by the Free Software Foundation.
 *
 *  bitpit is distributed in the hope that it will be useful, but WITHOUT
 *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 *  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
 *  License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with bitpit. If not, see <http://www.gnu.org/licenses/>.
 *
\*---------------------------------------------------------------------------*/
 
# ifndef __BITPIT_LEVELSET_OCTREE_KERNEL_HPP__
# define __BITPIT_LEVELSET_OCTREE_KERNEL_HPP__
 
#include "levelSetKernel.hpp"
 
#include "bitpit_voloctree.hpp"
 
namespace bitpit{
 
class LevelSetOctreeKernel : public LevelSetCachedKernel {
 
    private:
    std::size_t                                 m_cellCentroidCacheId;
 
    std::vector<double>                         m_octantTangentRadii ;    
    std::vector<double>                         m_octantBoundingRadii ;   
 
    public:
    template<typename value_t>
    using CellSparseCacheContainer = std::unordered_map<long, value_t>;
    template<typename value_t>
    using CellDenseCacheContainer = bitpit::PiercedStorage<value_t, long>;
 
    LevelSetOctreeKernel( VolOctree &patch, LevelSetFillIn fillIn );
 
    VolOctree *                                 getMesh() const override;
 
    double                                      getOctantTangentRadius(int level) const;
    double                                      getOctantBoundingRadius(int level) const;
 
    std::array<double, 3>                       computeCellCentroid(long) const override;
    double                                      computeCellTangentRadius(long) const override;
    double                                      computeCellBoundingRadius(long) const override;
 
};
 
// Typdefs for compatibility with older versions
typedef LevelSetOctreeKernel LevelSetOctree;
 
}
 
// Include template implementations
#include "levelSetOctreeKernel.tpp"
 
#endif