test_find_bug_LB.cpp

/*
 * test_find_LB_bug.cpp
 *
 *  Created on: 7 sept. 2015
 *      Author: marco
 */




#include "preprocessor_defines.dat"
#include "Class_Global.hpp"
#include "Class_Para_Tree.hpp"
#include "User_Data_Comm.hpp"

using namespace std;

// =================================================================================== //

int main(int argc, char *argv[]) {

#if NOMPI==0
        MPI::Init(argc, argv);

        {
#endif
                int iter = 0;
                int dim = 2;
                int ga = 2;

                Class_Para_Tree<2> pablo;
                pablo.setBalance((uint32_t)0,false);

                for (iter=0; iter<ga; iter++){
                        pablo.adaptGlobalRefine();
                }

#if NOMPI==0

                pablo.loadBalance();
#endif

                double xc, yc;
                xc = yc = 0.75;
                double radius = 0.25;

                uint32_t nocts = pablo.getNumOctants();
                uint32_t nghosts = pablo.getNumGhosts();

                for (int i=0; i<nocts; i++){
                        vector<vector<double> > nodes = pablo.getNodes(i);
                        for (int j=0; j<global2D.nnodes; 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))){
                                        pablo.setMarker(i,1);
                                }
                        }
                }
                pablo.adapt();
                if(pablo.rank==3 || pablo.rank==4)
                        pablo.setMarker((uint32_t)0,1);
                pablo.adapt();
                if(pablo.rank==3 || pablo.rank==4)
                        pablo.setMarker((uint32_t)0,1);
                pablo.adapt();
                if(pablo.rank==3)
                        pablo.setMarker((uint32_t)0,1);
                pablo.adapt();
                if(pablo.rank==3)
                        pablo.setMarker((uint32_t)0,1);
                pablo.adapt();
                pablo.loadBalance();

                iter = 0;
                pablo.updateConnectivity();
                pablo.write("PabloLB_iter"+to_string(static_cast<unsigned long long>(iter)));

#if NOMPI==0
        }

        MPI::Finalize();
#endif
}