POD_example_00003.cpp

/*---------------------------------------------------------------------------*\
 *
 *  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/>.
 *
\*---------------------------------------------------------------------------*/
 
 
#include <array>
#if BITPIT_ENABLE_MPI
#include <mpi.h>
#endif
 
#include "pod.hpp"
 
using namespace bitpit;
 
void run()
{
    POD pod;
   
    for (int i=0; i<10; i++)
        pod.addSnapshot("./data", "test."+std::to_string(i));
    
    pod.setMeshType(POD::MeshType::VOLOCTREE);
    pod.setStaticMesh(true); 
    pod.setUseMean(false);
    pod.setWriteMode(POD::WriteMode::DEBUG);
    pod.setMemoryMode(POD::MemoryMode::MEMORY_NORMAL);
    pod.setEnergyLevel(99);
    
    pod.setDirectory("pod");
    pod.setName("pod.test.solver");    
   
    for (int i=0; i<5; i++)
        pod.removeLeave1outSnapshot("./data", "test."+std::to_string(2*i));
     
    pod.leave1out();
    
    std::vector<std::string> namesf {"rho","rho_e"};
    std::vector<std::array<std::string,3>> namevf {}; //{{"u_x", "u_y","u_z"}};  
    pod.setTargetErrorFields(namesf,namevf);

    pod.setErrorThreshold(0.001);
    
    pod.evalErrorBoundingBox();
}

 
int main(int argc, char *argv[]) 
{
#if BITPIT_ENABLE_MPI
    MPI_Init(&argc,&argv);
#endif    

    try {
        run();
    } catch (const std::exception &exception) {
        log::cout() << exception.what();
        exit(1);
    }
 
#if BITPIT_ENABLE_MPI
    MPI_Finalize();
#endif
 
}