serialMultiDataField2D.hh

Go to the documentation of this file.

/* This file is part of the Palabos library.
 *
 * Copyright (C) 2011 FlowKit Sarl
 * Avenue de Chailly 23
 * 1012 Lausanne, Switzerland
 * E-mail contact: contact@flowkit.com
 *
 * The most recent release of Palabos can be downloaded at 
 * <http://www.palabos.org/>
 *
 * The library Palabos is free software: you can redistribute it and/or
 * modify it under the terms of the GNU Affero General Public License as
 * published by the Free Software Foundation, either version 3 of the
 * License, or (at your option) any later version.
 *
 * The library 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 Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef SERIAL_MULTI_DATA_FIELD_2D_HH
#define SERIAL_MULTI_DATA_FIELD_2D_HH

#include "serialMultiDataField2D.h"

/* *************** Class SerialScalarAccess2D ************************ */

namespace plb {
template<typename T>
SerialScalarAccess2D<T>::SerialScalarAccess2D()
    : locatedBlock(0)
{ }

template<typename T>
T& SerialScalarAccess2D<T>::getDistributedScalar (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,ScalarField2D<T>*>& fields )
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk (
            iX,iY, locatedBlock, localX, localY );
    PLB_PRECONDITION( ok );
    return fields[locatedBlock] -> get(localX,localY);
}

template<typename T>
T const& SerialScalarAccess2D<T>::getDistributedScalar (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,ScalarField2D<T>*> const& fields ) const
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk
            (iX,iY, locatedBlock, localX, localY);
    PLB_PRECONDITION( ok );
    return fields.find(locatedBlock)->second -> get(localX,localY);
}


template<typename T>
SerialScalarAccess2D<T>* SerialScalarAccess2D<T>::clone() const {
    return new SerialScalarAccess2D(*this);
}


/* *************** Class SerialTensorAccess2D ************************ */

template<typename T, int nDim>
SerialTensorAccess2D<T,nDim>::SerialTensorAccess2D()
    : locatedBlock(0)
{ }

template<typename T, int nDim>
Array<T,nDim>& SerialTensorAccess2D<T,nDim>::getDistributedTensor (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,TensorField2D<T,nDim>*>& fields )
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk
            (iX,iY, locatedBlock, localX, localY);
    PLB_PRECONDITION( ok );
    return fields[locatedBlock] -> get(localX,localY);
}

template<typename T, int nDim>
Array<T,nDim> const& SerialTensorAccess2D<T,nDim>::getDistributedTensor (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,TensorField2D<T,nDim>*> const& fields ) const
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk
            (iX,iY, locatedBlock, localX, localY);
    PLB_PRECONDITION( ok );
    return fields.find(locatedBlock)->second -> get(localX,localY);
}


template<typename T, int nDim>
SerialTensorAccess2D<T,nDim>* SerialTensorAccess2D<T,nDim>::clone() const {
    return new SerialTensorAccess2D(*this);
}



/* *************** Class SerialNTensorAccess2D ************************ */

template<typename T>
SerialNTensorAccess2D<T>::SerialNTensorAccess2D()
    : locatedBlock(0)
{ }

template<typename T>
T* SerialNTensorAccess2D<T>::getDistributedNTensor (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,NTensorField2D<T>*>& fields )
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk
            (iX,iY, locatedBlock, localX, localY);
    PLB_PRECONDITION( ok );
    typename std::map<plint,NTensorField2D<T>*>::const_iterator it = fields.find(locatedBlock);
    PLB_ASSERT( it != fields.end() );
    return it->second -> get(localX,localY);
}

template<typename T>
T const* SerialNTensorAccess2D<T>::getDistributedNTensor (
            plint iX, plint iY,
            MultiBlockManagement2D const& multiBlockManagement,
            std::map<plint,NTensorField2D<T>*> const& fields ) const
{
    plint localX, localY;
#ifdef PLB_DEBUG
    bool ok =
#endif
        multiBlockManagement.findInLocalBulk
            (iX,iY, locatedBlock, localX, localY);
    PLB_PRECONDITION( ok );
    typename std::map<plint,NTensorField2D<T>*>::const_iterator it = fields.find(locatedBlock);
    PLB_ASSERT( it != fields.end() );
    return it->second -> get(localX,localY);
}


template<typename T>
SerialNTensorAccess2D<T>* SerialNTensorAccess2D<T>::clone() const {
    return new SerialNTensorAccess2D(*this);
}

}  // namespace plb

#endif  // SERIAL_MULTI_DATA_FIELD_2D_HH