levelSetProxyObject.tpp

/*---------------------------------------------------------------------------*\
 *
 *  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_PROXY_OBJECT_TPP__
# define __BITPIT_LEVELSET_PROXY_OBJECT_TPP__
 
namespace bitpit {


template<typename SourceLevelSetObject, typename BaseLevelSetObject>
LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::LevelSetProxyObject(int id) : BaseLevelSetObject(id){
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
bool LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::isPrimary() const{
    return false;
} 

 * A cell can be either in the narrow band or in the bulk. It will be considered inside the narrow
 * band if one of the following conditions holds:
 *  - its distance from the surface is less than the narrow band size;
 *  - it intersects the zero-levelset iso-surface (intersections are checked using the
 *    FAST_GUARANTEE_FALSE criterium);
 *  - one of its neighbors intersects the zero-levelset iso-surface.
 */
template<typename SourceLevelSetObject, typename BaseLevelSetObject>
void LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::fillCellLocationCache()
{
    // Mesh information
    const VolumeKernel &mesh = *(this->getKernel()->getMesh()) ;
    VolumeKernel::CellConstIterator cellBegin = mesh.cellConstBegin();
    VolumeKernel::CellConstIterator cellEnd   = mesh.cellConstEnd();
 
    // Get cell zone cache
    typedef typename BaseLevelSetObject::CellCacheCollection::template ValueCache<char> BaseZoneCache;
    BaseZoneCache *locationCache = this->template getCellCache<char>(this->m_cellLocationCacheId);
 
    // Get zone information from reference object
    for (VolumeKernel::CellConstIterator cellItr = cellBegin; cellItr != cellEnd; ++cellItr) {
        long cellId = cellItr.getId();
        LevelSetCellLocation cellLocation = getCellReferenceObject(cellId)->getCellLocation(cellId);
        locationCache->insertEntry(cellId, static_cast<char>(cellLocation));
    }
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
void LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::fillCellLocationCache(const std::vector<adaption::Info> &adaptionData)
{
    // Get cell location cache
    typedef typename BaseLevelSetObject::CellCacheCollection::template ValueCache<char> BaseZoneCache;
    BaseZoneCache *locationCache = this->template getCellCache<char>(this->m_cellLocationCacheId);
 
    // Get location information from reference object
    for (const adaption::Info &adaptionInfo : adaptionData) {
        if (adaptionInfo.entity != adaption::Entity::ENTITY_CELL) {
            continue;
        }
 
        // Skip received cells when the cache is non-volatile
        //
        // If the cache is non-volatile, data on exchanged cells has been communicated during
        // cache adaption.
        if (adaptionInfo.type == adaption::Type::TYPE_PARTITION_RECV && !locationCache->isVolatile()) {
            continue;
        }
 
        // Add current cells to the process list
        for (long cellId : adaptionInfo.current) {
            LevelSetCellLocation cellLocation = getCellReferenceObject(cellId)->getCellLocation(cellId);
            locationCache->insertEntry(cellId, static_cast<char>(cellLocation));
        }
    }
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
bool LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::isCellInNarrowBand(long id)const
{
    return getCellReferenceObject(id)->isCellInNarrowBand(id);
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
bool LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::isInNarrowBand(const std::array<double,3> &point)const
{
    return getReferenceObject(point)->isInNarrowBand(point);
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
const SourceLevelSetObject * LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getCellReferencePrimaryObject(long id) const{
 
    const SourceLevelSetObject *referenceObject = getCellReferenceObject(id);
    if (!referenceObject) {
        return nullptr;
    }
 
    if (referenceObject->isPrimary()) {
        return referenceObject;
    }
 
    if( const LevelSetProxyObject *referenceProxyObject = dynamic_cast<const LevelSetProxyObject*>(referenceObject) ){
        return referenceProxyObject->getCellReferencePrimaryObject(id);
    }
 
    return nullptr;
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
const SourceLevelSetObject * LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getReferencePrimaryObject(const std::array<double, 3> &point) const{
 
    const SourceLevelSetObject *referenceObject = getReferenceObject(point);
    if (!referenceObject) {
        return nullptr;
    }
 
    if (referenceObject->isPrimary()) {
        return referenceObject;
    }
 
    if( const LevelSetProxyObject *referenceProxyObject = dynamic_cast<const LevelSetProxyObject*>(referenceObject) ){
        return referenceProxyObject->getReferencePrimaryObject(point);
    }
 
    return nullptr;
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
int LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getCellReferenceObjectId(long id) const{
 
    const SourceLevelSetObject *referenceObject = getCellReferenceObject(id);
    if (!referenceObject) {
        return levelSetDefaults::OBJECT;
    }
 
    return referenceObject->getId();
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
int LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getCellReferencePrimaryObjectId(long id) const{
 
    const SourceLevelSetObject *referenceObject = getCellReferenceObject(id);
    if (!referenceObject) {
        return levelSetDefaults::OBJECT;
    }
 
    if (referenceObject->isPrimary()) {
        return referenceObject->getId();
    }
 
    const LevelSetProxyBaseObject *referenceProxyObject = dynamic_cast<const LevelSetProxyBaseObject *>(referenceObject);
    if (!referenceProxyObject) {
        return levelSetDefaults::OBJECT;
    }
 
    return referenceProxyObject->getCellReferencePrimaryObjectId(id);
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
int LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getReferenceObjectId(const std::array<double, 3> &point) const{
 
    const SourceLevelSetObject *referenceObject = getReferenceObject(point);
    if (!referenceObject) {
        return levelSetDefaults::OBJECT;
    }
 
    return referenceObject->getId();
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
int LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getReferencePrimaryObjectId(const std::array<double, 3> &point) const{
 
    const SourceLevelSetObject *referenceObject = getReferenceObject(point);
    if (!referenceObject) {
        return levelSetDefaults::OBJECT;
    }
 
    if (referenceObject->isPrimary()) {
        return referenceObject->getId();
    }
 
 
    const LevelSetProxyBaseObject *referenceProxyObject = dynamic_cast<const LevelSetProxyBaseObject *>(referenceObject);
    if (!referenceProxyObject) {
        return levelSetDefaults::OBJECT;
    }
 
    return referenceProxyObject->getReferencePrimaryObjectId(point);
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
std::vector<const SourceLevelSetObject *> LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getPrimarySourceObjects() const{
 
    std::vector<const SourceLevelSetObject *> objects;
    for( const SourceLevelSetObject *sourceObject : getSourceObjects()){
        if( const LevelSetProxyObject *proxySourceObject = dynamic_cast<const LevelSetProxyObject*>(sourceObject) ){
            std::vector<const SourceLevelSetObject *> sourcePrimarySourceObjects = proxySourceObject->getPrimarySourceObjects();
            objects.insert(objects.end(), sourcePrimarySourceObjects.begin(), sourcePrimarySourceObjects.end());
        } else {
            objects.push_back(sourceObject);
        }
    }
 
    return objects;
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
std::vector<int> LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getSourceObjectIds() const{
 
    std::vector<const SourceLevelSetObject*> sourceObjects = getSourceObjects();
    std::size_t nSourceObjects = sourceObjects.size();
 
    std::vector<int> sourceIds(nSourceObjects);
    for (std::size_t i = 0; i < nSourceObjects; ++i) {
        const SourceLevelSetObject *sourceObject = sourceObjects[i];
        sourceIds[i] = sourceObject->getId();
    }
 
    return sourceIds;
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
std::vector<int> LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getPrimarySourceObjectIds() const{
 
    std::vector<const SourceLevelSetObject*> primaryObjects = getPrimarySourceObjects();
    std::size_t nPrimaryObjects = primaryObjects.size();
 
    std::vector<int> primaryIds(nPrimaryObjects);
    for (std::size_t i = 0; i < nPrimaryObjects; ++i) {
        const SourceLevelSetObject *primaryObject = primaryObjects[i];
        primaryIds[i] = primaryObject->getId();
    }
 
    return primaryIds;
 
}

template<typename SourceLevelSetObject, typename BaseLevelSetObject>
int LevelSetProxyObject<SourceLevelSetObject, BaseLevelSetObject>::getPrimaryObjectId(long id) const{
 
    return getReferencePrimaryObjectId(id);
 
}
 
}
 
#endif