test12.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

        Class_Para_Tree<2> pablo12;

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

        pablo12.computeConnectivity();
        pablo12.write("Pablo12_iter0");

        for (int iter=1; iter<3; iter++){
            pablo12.adaptGlobalRefine();

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

#if NOMPI==0

        pablo12.loadBalance();
#endif

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

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

            pablo12.adapt();

#if NOMPI==0

            pablo12.loadBalance();
#endif

            pablo12.updateConnectivity();
            pablo12.write("Pablo12_iter"+to_string(static_cast<unsigned long long>(iter)));
        }
#if NOMPI==0
    }

    MPI::Finalize();
#endif
}