levelSetCartesianKernel.cpp

/*---------------------------------------------------------------------------*\
 *
 *  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/>.
 *
\*---------------------------------------------------------------------------*/
 
# include "bitpit_volcartesian.hpp"
 
# include "levelSetKernel.hpp"
# include "levelSetCartesianKernel.hpp"
 
namespace bitpit {


LevelSetCartesianKernel::LevelSetCartesianKernel(VolCartesian &patch, LevelSetFillIn fillIn ): LevelSetCachedKernel( &patch, fillIn ){
 
    BITPIT_UNUSED(fillIn);
 
    // Get mesh information
    const VolCartesian *mesh = getMesh();
    int dimension = mesh->getDimension();
    std::array<double, 3> spacing = mesh->getSpacing();
 
    // Initialize bounding and tangent radii
    m_cellTangentRadius  = std::numeric_limits<double>::max();
    m_cellBoundingRadius = 0.;
    for(int i=0; i<dimension; ++i){
        m_cellTangentRadius   = std::min(m_cellTangentRadius, spacing[i]);
        m_cellBoundingRadius += spacing[i] * spacing[i];
    }
    m_cellTangentRadius  = 0.5 * m_cellTangentRadius;
    m_cellBoundingRadius = 0.5 * std::sqrt(m_cellBoundingRadius);
 
}

VolCartesian * LevelSetCartesianKernel::getMesh() const{
 
    return static_cast<VolCartesian *>(LevelSetKernel::getMesh()) ;
 
}

double LevelSetCartesianKernel::getCellTangentRadius( ) const {
 
    return m_cellTangentRadius;
 
}

double LevelSetCartesianKernel::getCellBoundingRadius( ) const {
 
    return m_cellBoundingRadius;
 
}

std::array<double, 3> LevelSetCartesianKernel::computeCellCentroid( long id ) const {
 
    const VolCartesian *mesh = getMesh();
 
    return mesh->evalCellCentroid( id );
}

double LevelSetCartesianKernel::computeCellTangentRadius( long id ) const {
 
    BITPIT_UNUSED( id );
 
    return getCellTangentRadius();
 
}

double LevelSetCartesianKernel::computeCellBoundingRadius( long id ) const {
 
    BITPIT_UNUSED( id );
 
    return getCellBoundingRadius();
 
}
 
}