ioofoam_example_00002.cpp

Example of reading,morphing and writing of a parallel OpenFOAM case mesh.Mesh is read from an OpenFOAM case. In writing, boundary mesh is written on output.

Using: IOOFOAM, MimmoGeometry, Chain

To run (MPI version only): mpirun -np 8 ioofoam_example_00002
Beware: please use np = 8 procs, since the reference sample mesh is already partitioned on 8 processors. Otherwise to use other np's but re-partition the OpenFoam sample mesh first.

visit: mimmo website

/*---------------------------------------------------------------------------*\
 *
 *  mimmo
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of mimmo.
 *
 *  mimmo 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.
 *
 *  mimmo 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 mimmo. If not, see <http://www.gnu.org/licenses/>.
 *
 \ *---------------------------------------------------------------------------*/
 
#include "mimmo_iogeneric.hpp"
#include "IOOFOAM.hpp"
#include "mimmo_manipulators.hpp"
 
// =================================================================================== //
 
void OFOAM_readwrite() {
 
    /* read a OpenFoam mesh, expose bulk volume and boundary mesh */
    mimmo::IOOFOAM * reader = new mimmo::IOOFOAM(false);
    reader->setDir("geodata/ofoam.00002");
 
    /* write surface boundary mesh to file */
    mimmo::MimmoGeometry * surfacewriter = new mimmo::MimmoGeometry(mimmo::MimmoGeometry::IOMode::WRITE);
    surfacewriter->setWriteDir(".");
    surfacewriter->setWriteFileType(FileType::SURFVTU);
    surfacewriter->setWriteFilename("ioofoam.00002");
 
    /* write bulk volume mesh to file */
    mimmo::MimmoGeometry * volumewriter = new mimmo::MimmoGeometry(mimmo::MimmoGeometry::IOMode::WRITE);
    volumewriter->setWriteDir(".");
    volumewriter->setWriteFileType(FileType::VOLVTU);
    volumewriter->setWriteFilename("ioofoam.00003");
 
    /* write openfoam mesh to file */
    mimmo::IOOFOAM * writer = new mimmo::IOOFOAM(true);
    writer->setDir("geodata/ofoam.00002");
 
    /* define block pin connections */
    mimmo::pin::addPin(reader, writer, M_GEOMOFOAM2, M_GEOMOFOAM2);
    mimmo::pin::addPin(reader, surfacewriter, M_GEOMOFOAM2, M_GEOM);
    mimmo::pin::addPin(reader, volumewriter, M_GEOMOFOAM, M_GEOM);
    mimmo::pin::addPin(reader, writer, M_GEOMOFOAM, M_GEOMOFOAM);
    mimmo::pin::addPin(reader, writer, M_UMAPIDS, M_UMAPIDS);
 
    /* define chain execution */
    mimmo::Chain c0;
    c0.addObject(reader);
    c0.addObject(surfacewriter);
    c0.addObject(volumewriter);
    c0.addObject(writer);
    c0.exec(true);
 
    /*clean up */
    delete reader;
    delete surfacewriter;
    delete volumewriter;
    delete writer;
}
 
// =================================================================================== //
 
int main( int argc, char *argv[] ) {
 
    BITPIT_UNUSED(argc);
    BITPIT_UNUSED(argv);
 
#if MIMMO_ENABLE_MPI
    MPI_Init(&argc, &argv);
#endif
 
        try{
             OFOAM_readwrite() ;
        }
        catch(std::exception & e){
            std::cout<<"ioofoam_example_00002 exit with the following errors :"<<e.what()<<std::endl;
            return 1;
        }
 
#if MIMMO_ENABLE_MPI
    MPI_Finalize();
#endif
 
    return 0;
}