testExportImport.cpp

Go to the documentation of this file.

//@HEADER
/*
*******************************************************************************

    Copyright (C) 2004, 2005, 2007 EPFL, Politecnico di Milano, INRIA
    Copyright (C) 2010 EPFL, Politecnico di Milano, Emory University

    This file is part of LifeV.

    LifeV is free software; you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    LifeV 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 LifeV.  If not, see <http://www.gnu.org/licenses/>.

*******************************************************************************
*/
//@HEADER

/*!
    @file
    @brief test for ExporterEnsight

    @author Tiziano Passerini <tiziano@mathcs.emory.edu>
    @contributor Simone Deparis
    @maintainer

    @date 13-1-2011

 */


#include <Epetra_ConfigDefs.h>
#include <Epetra_Comm.h>


#include <lifev/core/LifeV.hpp>
#include "../importExport/RossEthierSteinmanDec.hpp"
#include "../importExport/TestImportExport.hpp"

#include <lifev/core/filter/ExporterEnsight.hpp>
#include <lifev/core/filter/ExporterVTK.hpp>
#include <lifev/core/filter/ExporterHDF5.hpp>


typedef LifeV::ExporterEnsight<TestImportExport::mesh_Type> ensight_Type;
typedef LifeV::ExporterHDF5<TestImportExport::mesh_Type> hdf5_Type;
typedef LifeV::ExporterVTK<TestImportExport::mesh_Type> vtk_Type;



using namespace LifeV;

int
main ( int argc, char** argv )
{
    //MPI communicator initialization
    std::shared_ptr<Epetra_Comm> commPtr;

#ifdef HAVE_MPI
    std::cout << "MPI Initialization" << std::endl;
    MPI_Init ( &argc, &argv );
#endif

    //MPI Preprocessing
#ifdef EPETRA_MPI

    int nprocs;
    int rank;

    MPI_Comm_size ( MPI_COMM_WORLD, &nprocs );
    MPI_Comm_rank ( MPI_COMM_WORLD, &rank );

    if ( rank == 0 )
    {
        std::cout << "MPI processes: " << nprocs << std::endl;
        std::cout << "MPI Epetra Initialization ... " << std::endl;
    }
    commPtr.reset ( new Epetra_MpiComm ( MPI_COMM_WORLD ) );

    commPtr->Barrier();

#else

    std::cout << "MPI SERIAL Epetra Initialization ... " << std::endl;
    commPtr.reset ( new Epetra_SerialComm() );

#endif

    GetPot command_line (argc, argv);
    TestImportExport testImportExport ( commPtr );

    bool passed (false);

    const std::string exporterName = command_line.follow ("ensight", 2, "-e", "--exporter");

    if ( exporterName.compare ( "vtk" ) == 0 )
    {
    passed = testImportExport.run<vtk_Type, vtk_Type > ( command_line, "export" );
    }
    else if ( exporterName.compare ( "ensight" ) == 0 )
    {
    passed = testImportExport.run<ensight_Type, ensight_Type > ( command_line, "export" );
    }
#ifdef HAVE_HDF5
    else
    {
    passed = testImportExport.run<hdf5_Type, hdf5_Type > ( command_line, "export" );
    }
#endif

    // ----- End of test calls -----

#ifdef HAVE_MPI
    std::cout << "MPI Finalization" << std::endl;
    MPI_Finalize();
#endif

    if (passed)
    {
        return EXIT_SUCCESS;
    }
    else
    {
        return EXIT_FAILURE;
    }
}