offLatticeBoundaryProcessor3D.h

00001 /* This file is part of the Palabos library.
00002  *
00003  * Copyright (C) 2011 FlowKit Sarl
00004  * Avenue de Chailly 23
00005  * 1012 Lausanne, Switzerland
00006  * E-mail contact: contact@flowkit.com
00007  *
00008  * The most recent release of Palabos can be downloaded at 
00009  * <http://www.palabos.org/>
00010  *
00011  * The library Palabos is free software: you can redistribute it and/or
00012  * modify it under the terms of the GNU Affero General Public License as
00013  * published by the Free Software Foundation, either version 3 of the
00014  * License, or (at your option) any later version.
00015  *
00016  * The library is distributed in the hope that it will be useful,
00017  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  * GNU Affero General Public License for more details.
00020  *
00021  * You should have received a copy of the GNU Affero General Public License
00022  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
00023 */
00024 
00025 #ifndef OFF_LATTICE_BOUNDARY_PROCESSOR_3D_H
00026 #define OFF_LATTICE_BOUNDARY_PROCESSOR_3D_H
00027 
00028 #include "core/globalDefs.h"
00029 #include "core/array.h"
00030 #include "atomicBlock/dataProcessingFunctional3D.h"
00031 #include "atomicBlock/atomicContainerBlock3D.h"
00032 #include "dataProcessors/metaStuffHelper.h"
00033 #include "offLattice/boundaryShapes3D.h"
00034 #include <vector>
00035 #include <utility>
00036 
00037 namespace plb {
00038 
00039 template<typename T>
00040 class CheckVoxelizationFunctional3D : public PlainReductiveBoxProcessingFunctional3D
00041 {
00042 public:
00043     CheckVoxelizationFunctional3D(BoundaryShape3D<T,Array<T,3> >* shape_, int flowType_);
00044     virtual ~CheckVoxelizationFunctional3D();
00045     CheckVoxelizationFunctional3D(CheckVoxelizationFunctional3D const& rhs);
00046     CheckVoxelizationFunctional3D& operator= (
00047             CheckVoxelizationFunctional3D const& rhs );
00048 
00052     virtual void processGenericBlocks(Box3D domain, std::vector<AtomicBlock3D*> fields);
00053     virtual CheckVoxelizationFunctional3D<T>* clone() const;
00054     virtual void getTypeOfModification(std::vector<modif::ModifT>& modified) const;
00055     virtual BlockDomain::DomainT appliesTo() const;
00056     plint getNumErrors() const;
00057 private:
00058     bool isFluid(Dot3D const& location) const;
00059     void computeCell(Dot3D const& cellLocation,
00060                      ScalarField3D<int>& flagMatrix,
00061                      Dot3D const& offset);
00062 private:
00063     plint numErrorsId;
00064     BoundaryShape3D<T,Array<T,3> >* shape;
00065     int flowType;
00066 };
00067 
00068 }  // namespace plb
00069 
00070 #endif  // OFF_LATTICE_BOUNDARY_PROCESSOR_3D_H