PABLO_example_00005.cpp

Parallel 2D adaptive mesh refinement (AMR) with data using PABLO.

Parallel 2D adaptive mesh refinement (AMR) with data using PABLOThis example is the parallel version of PABLO_example_00002.

Here one quadtree is instantiated in the same way on every process. Therefore, each process refines globally the initial quadrant 3 times. Till this moment, there's no parallel paradigm in action, actually the code is replicated on every process.

At this point a load-balance method call is performed and the grid is partitioned and distributed among the processes of the world communicator. From now on, each process owns one portion of the grid.

The test continues as in example 00002, but after every refinement/coarsening a load-balance call is introduced to keep the computational burden well distributed.

To run: ./PABLO_example_00005
To see the result visit: PABLO website

/*---------------------------------------------------------------------------*\
 *
 *  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/>.
 *
\*---------------------------------------------------------------------------*/
 
#if BITPIT_ENABLE_MPI==1
#include <mpi.h>
#endif
 
#include "bitpit_PABLO.hpp"
 
using namespace std;
using namespace bitpit;
 
// =================================================================================== //
// =================================================================================== //
 
void run()
{
    int iter = 0;
 
    PabloUniform pablo5(2);
 
    uint32_t idx=0;
    pablo5.setBalance(idx,false);
 
    pablo5.computeConnectivity();
    pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
 
    for (iter=1; iter<3; iter++){
        pablo5.adaptGlobalRefine();
        pablo5.updateConnectivity();
        pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
    }
 
#if BITPIT_ENABLE_MPI==1
    pablo5.loadBalance();
#endif
 
    double xc, yc;
    xc = yc = 0.5;
    double radius = 0.4;
 
    for (iter=3; iter<9; iter++){
        uint32_t nocts = pablo5.getNumOctants();
        for (unsigned int i=0; i<nocts; i++){
            vector<array<double,3> > nodes = pablo5.getNodes(i);
            for (int j=0; j<4; j++){
                double x = nodes[j][0];
                double y = nodes[j][1];
                if ((pow((x-xc),2.0)+pow((y-yc),2.0) <= pow(radius,2.0))){
                    pablo5.setMarker(i, 1);
                }
            }
        }
 
        pablo5.adapt();
 
#if BITPIT_ENABLE_MPI==1
        pablo5.loadBalance();
#endif
 
        pablo5.updateConnectivity();
        pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
    }
 
    pablo5.adaptGlobalCoarse();
    pablo5.updateConnectivity();
    pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
 
 
    xc = yc = 0.35;
    radius = 0.15;
 
    for (iter=10; iter<15; iter++){
        uint32_t nocts = pablo5.getNumOctants();
        for (unsigned int i=0; i<nocts; i++){
            vector<array<double,3> > nodes = pablo5.getNodes(i);
            for (int j=0; j<4; j++){
                double x = nodes[j][0];
                double y = nodes[j][1];
                if ((pow((x-xc),2.0)+pow((y-yc),2.0) <= pow(radius,2.0))){
                    pablo5.setMarker(i, -1);
                }
            }
        }
        pablo5.adapt();
 
#if BITPIT_ENABLE_MPI==1
        pablo5.loadBalance();
#endif
 
        pablo5.updateConnectivity();
        pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
    }
 
    pablo5.adaptGlobalCoarse();
    pablo5.updateConnectivity();
    pablo5.write("pablo00005_iter"+to_string(static_cast<unsigned long long>(iter)));
}
 
int main(int argc, char *argv[])
{
#if BITPIT_ENABLE_MPI==1
    MPI_Init(&argc,&argv);
#else
    BITPIT_UNUSED(argc);
    BITPIT_UNUSED(argv);
#endif
 
    int nProcs;
    int rank;
#if BITPIT_ENABLE_MPI==1
    MPI_Comm_size(MPI_COMM_WORLD, &nProcs);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
#else
    nProcs = 1;
    rank   = 0;
#endif
 
    // Initialize the logger
    log::manager().initialize(log::MODE_SEPARATE, false, nProcs, rank);
    log::cout() << log::fileVerbosity(log::LEVEL_INFO);
    log::cout() << log::disableConsole();
 
    // Run the example
    try {
        run();
    } catch (const std::exception &exception) {
        log::cout() << exception.what();
        exit(1);
    }
 
#if BITPIT_ENABLE_MPI==1
    MPI_Finalize();
#endif
}