test3.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> pablo3;

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

        pablo3.computeConnectivity();
        pablo3.write("Pablo3_iter"+to_string(static_cast<unsigned long long>(iter)));

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

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

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

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

        pablo3.adaptGlobalCoarse();
        pablo3.updateConnectivity();
        pablo3.write("Pablo3_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 = pablo3.getNumOctants();
            for (int i=0; i<nocts; i++){
                vector<vector<double> > nodes = pablo3.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))){
                        pablo3.setMarker(i, -1);
                    }
                }
            }
            pablo3.adapt();
            pablo3.updateConnectivity();
            pablo3.write("Pablo3_iter"+to_string(static_cast<unsigned long long>(iter)));
        }
#if NOMPI==0
    }
    MPI::Finalize();
#endif
}