test13.cpp

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

using namespace std;

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

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

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

    {
#endif
        int iter = 0;

        Class_Para_Tree<2> pablo13;

        uint32_t idx=0;
        pablo13.setBalance(idx,false);

        pablo13.computeConnectivity();
        pablo13.write("Pablo13_iter"+to_string(static_cast<unsigned long long>(iter)));

        for (iter=1; iter<3; iter++){
            pablo13.adaptGlobalRefine();
            pablo13.updateConnectivity();
            pablo13.write("Pablo13_iter"+to_string(static_cast<unsigned long long>(iter)));
        }

#if NOMPI==0

        pablo13.loadBalance();
#endif

        double xc, yc;
        xc = yc = 0.5;
        double radius = 0.4;

        for (iter=3; iter<9; iter++){
            uint32_t nocts = pablo13.getNumOctants();
            for (int i=0; i<nocts; i++){
                vector<vector<double> > nodes = pablo13.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))){
                        pablo13.setMarker(i, 1);
                    }
                }
            }

            pablo13.adapt();

#if NOMPI==0

            pablo13.loadBalance();
#endif

            pablo13.updateConnectivity();
            pablo13.write("Pablo13_iter"+to_string(static_cast<unsigned long long>(iter)));
        }

        pablo13.adaptGlobalCoarse();
        pablo13.updateConnectivity();
        pablo13.write("Pablo13_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 = pablo13.getNumOctants();
            for (int i=0; i<nocts; i++){
                vector<vector<double> > nodes = pablo13.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))){
                        pablo13.setMarker(i, -1);
                    }
                }
            }
            pablo13.adapt();

#if NOMPI==0

            pablo13.loadBalance();
#endif

            pablo13.updateConnectivity();
            pablo13.write("Pablo13_iter"+to_string(static_cast<unsigned long long>(iter)));
        }
#if NOMPI==0
    }
    MPI::Finalize();
#endif
}