Projektarbeit: Discontinuous Galerkin (first order) with hyper.deal

adjust_examples.py

@@ -9,7 +9,7 @@
 
 dir_path = os.path.dirname(os.path.realpath(__file__))
 GLPK_DIR= str(dir_path) + "/glpk"
-examples = ["combi_example", "combi_example_faults", "gene_distributed", "gene_distributed_linear"]
+examples = ["combi_example", "combi_example_faults", "gene_distributed", "gene_distributed_linear","dealii_example"]
 for example in examples:
     pfilein = open(str(dir_path)+ "/distributedcombigrid/examples/" + example + "/Makefile.template" ,'r')
     temp = pfilein.read()

compile.sh

@@ -1,2 +1,2 @@
 export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$(pwd)/glpk/lib
-scons -j 8 ISGENE=0 VERBOSE=1 COMPILE_BOOST_TESTS=1 RUN_BOOST_TESTS=1 RUN_CPPLINT=0 BUILD_STATICLIB=0 CC=mpicc.mpich FC=mpifort.mpich CXX=mpicxx.mpich OPT=1 TIMING=0 UNIFORMDECOMPOSITION=1 ENABLEFT=0 #DEBUG_OUTPUT=1
+scons -j 20 ISGENE=0 VERBOSE=1 COMPILE_BOOST_TESTS=1 RUN_BOOST_TESTS=1 RUN_CPPLINT=0 BUILD_STATICLIB=0 CC=mpicc FC=mpifort CXX=mpicxx OPT=1 TIMING=0 UNIFORMDECOMPOSITION=1 ENABLEFT=0 DOC=0  #DEBUG_OUTPUT=1

distributedcombigrid/examples/combi_example/DataConverter.cpp

@@ -0,0 +1,252 @@
+#include <string>
+#include <vector>
+#include <math.h>
+#include <sstream>
+#include "sgpp/distributedcombigrid/combischeme/CombiMinMaxScheme.hpp"
+
+#include "sgpp/distributedcombigrid/loadmodel/LinearLoadModel.hpp"
+#include "sgpp/distributedcombigrid/manager/CombiParameters.hpp"
+#include <boost/property_tree/ptree.hpp>
+#include <boost/property_tree/json_parser.hpp>
+
+#include "sgpp/distributedcombigrid/utils/Types.hpp"
+#include "DataConverter.hpp"
+
+using namespace combigrid;
+
+
+    Converter::Converter(){};
+    void Converter::toParam(std::string _filename, std::string newname){
+        namespace pt = boost::property_tree;
+
+        // Create a root
+        pt::ptree root,ini_root,ini_ct, ini_application, ini_manager;
+
+        // Load the json file in this ptree
+        try{
+            pt::read_json(_filename, root);
+            int dimx=root.get<int>("General.DimX");
+            int dimv=root.get<int>("General.DimV");
+            int dim=dimx+dimv;
+            //lmax:
+            int refineX=root.get<int>("Case.NRefinementsX");
+            int refineV=root.get<int>("Case.NRefinementsV");
+            int subX=root.get<int>("Case.NSubdivisionsX.X");
+            int subV=root.get<int>("Case.NSubdivisionsV.X");
+            LevelVector lma;
+
+            //zur Basis 2 hoch beziehen!
+            for(int i=0; i<dimx;i++){
+                lma.push_back(pow(subX*(refineX+1),1.0/3.0)); 
+           }
+            for(int i=0; i<dimx;i++){
+                lma.push_back(pow(1.0/3.0,subV*(refineV+1)));
+            }
+
+            //lmin=lmax
+
+
+            //ncombi=10
+
+            //p for partitions
+            int partX=root.get<int>("General.PartitionX");
+            int partY=root.get<int>("General.PartitionV");
+            LevelVector p_{partX,partY};
+
+            //boundary
+
+            LevelVector bound;
+            for(int i=0;i<dimx;i++){
+                bound.push_back(true);
+            }
+            for(int i=0;i<dimv;i++){
+                bound.push_back(true);
+            }
+            std::string levelString;
+            std::ostringstream vec,p;
+            vec << lma;
+            levelString=vec.str();
+            p<<p_[0]<<" "<<p_[1]<<std::endl;
+            //leval=lmax;
+            ini_ct.put("dim",dim);
+            ini_ct.put("lmin",levelString);
+            ini_ct.put("lmax",levelString);
+            ini_ct.put("leval",levelString);
+            ini_ct.put("p",p.str());
+            ini_ct.put("ncombi",10);
+
+            ini_root.add_child("ct",ini_ct);
+            ini_application.put("dt",0.001);
+            ini_application.put("steps",100);
+            ini_root.add_child("application",ini_application);
+            ini_manager.put("ngroup",2);
+            ini_manager.put("nprocs",1);
+            ini_root.add_child("manager",ini_manager);
+            pt::write_ini(newname,ini_root);
+
+        }
+        catch(pt::json_parser::json_parser_error error){
+            std::cout<< "Not a valid json File";
+        }
+        catch(...){
+            std::cout << "Something went wrong";
+        }
+    }
+
+    //converts the data from the ini file to a json format and stores the data in newname.json
+    void Converter::toJSON(std::string _fileName, std::string newname){
+        namespace pt = boost::property_tree;
+
+        // Create a root
+        pt::ptree root;
+        pt::ptree json_root, json_general, json_case, json_SpatialDiscretization, json_TemporalDiscretization, json_postprocessing;
+        pt::ptree jsubX, jsubV,jsout, jVTK;
+
+        // Load the ini file in this ptree
+        try{
+            pt::read_ini(_fileName, root);
+            std::cout << "Name "<<_fileName;
+            int dim=root.get<int>("ct.dim",2);
+            int dimx=root.get<int>("ct.dimx",dim/2);
+            int dimv=root.get<int>("ct.dimv",dim/2);
+
+            //degree is 1 per default
+            int degx=root.get<int>("ct.degreex",1);
+            int degv=root.get<int>("ct.degreev",1);
+
+            // case is per default hyperrectangle and periodic is per default 
+            //nrefinments is set to 0 and the distribution is controlled via subdivisions
+
+            //lmax:
+            LevelVector lmax(dimx+dimv),p(dimx+dimv);
+            root.get<std::string>("ct.lmax") >> lmax;
+            root.get<std::string>("ct.p") >> p;
+
+            json_general.put("DimX", dimx);
+            json_general.put("DimV", dimv);
+            json_general.put("DegreeX", degx);
+            json_general.put("DegreeV", degx);
+            json_general.put("PartitionX", p[0]);
+            json_general.put("PartitionV", p[1]);
+            json_general.put("Case","hyperrectangle");
+            json_root.add_child("General", json_general);
+
+            json_case.put("NRefinementsX",0);
+            json_case.put("NRefinementsV",0);
+            json_case.put("PeriodicX",true);
+            json_case.put("PeriodicV",true);
+            jsubX.put("X",pow(2,lmax[0]));
+            jsubX.put("Y",pow(2,lmax[0]));
+            jsubX.put("Z",pow(2,lmax[0]));
+            jsubV.put("X",pow(2,lmax[1]));
+            jsubV.put("Y",pow(2,lmax[1]));
+            jsubV.put("Z",pow(2,lmax[1]));
+            json_case.add_child("NSubdivisionsX",jsubX);
+            json_case.add_child("NSubdivisionsV",jsubV);
+            json_root.add_child("Case", json_case);
+
+            json_SpatialDiscretization.put("TriangulationType","FullyDistributed");
+            json_SpatialDiscretization.put("MappingX",1);
+            json_SpatialDiscretization.put("MappingV",1);
+            json_root.add_child("SpatialDiscretization",json_SpatialDiscretization);
+
+            json_TemporalDiscretization.put("FinalTime",2);
+            json_TemporalDiscretization.put("CFLNumber",0.15);
+            json_root.add_child("TemporalDiscretization",json_TemporalDiscretization);
+
+            jsout.put("Tick",0.1);
+            json_postprocessing.add_child("StandardOutput",jsout);
+            jVTK.put("Enabled",false);
+            json_postprocessing.add_child("VTK",jVTK);
+
+            json_root.add_child("Postprocessing", json_postprocessing);
+            pt::write_json(newname,json_root);
+
+
+        }
+        catch(pt::ini_parser::ini_parser_error error){
+            std::cout<< "Not a valid file";
+        }
+        catch(...){
+            std::cout << "Something went wrong";
+        }
+    }
+
+    void Converter::toJSON(std::string _fileName, std::string newname,LevelVector l){
+        namespace pt = boost::property_tree;
+
+        // Create a root
+        pt::ptree root;
+        pt::ptree json_root, json_general, json_case, json_SpatialDiscretization, json_TemporalDiscretization, json_postprocessing;
+        pt::ptree jsubX, jsubV,jsout, jVTK;
+
+        // Load the ini file in this ptree
+        try{
+            pt::read_ini(_fileName, root);
+            std::cout << "Name "<<_fileName;
+            int dim=root.get<int>("ct.dim",2);
+            int dimx=root.get<int>("ct.dimx",dim/2);
+            int dimv=root.get<int>("ct.dimv",dim/2);
+
+            //degree is 1 per default
+            int degx=root.get<int>("ct.degreex",1);
+            int degv=root.get<int>("ct.degreev",1);
+
+            // case is per default hyperrectangle and periodic is per default 
+            //nrefinments is set to 0 and the distribution is controlled via subdivisions
+
+            //lmax:
+            LevelVector lmax(dimx+dimv),p(dimx+dimv);
+            lmax=l;
+            root.get<std::string>("ct.p") >> p;
+
+            json_general.put("DimX", dimx);
+            json_general.put("DimV", dimv);
+            json_general.put("DegreeX", degx);
+            json_general.put("DegreeV", degx);
+            json_general.put("PartitionX", p[0]);
+            json_general.put("PartitionV", p[1]);
+            json_general.put("Case","hyperrectangle");
+            json_root.add_child("General", json_general);
+
+            json_case.put("NRefinementsX",0);
+            json_case.put("NRefinementsV",0);
+            json_case.put("PeriodicX",true);
+            json_case.put("PeriodicV",true);
+            jsubX.put("X",pow(2,lmax[0]));
+            jsubX.put("Y",pow(2,lmax[0]));
+            jsubX.put("Z",pow(2,lmax[0]));
+            jsubV.put("X",pow(2,lmax[1]));
+            jsubV.put("Y",pow(2,lmax[1]));
+            jsubV.put("Z",pow(2,lmax[1]));
+            json_case.add_child("NSubdivisionsX",jsubX);
+            json_case.add_child("NSubdivisionsV",jsubV);
+            json_root.add_child("Case", json_case);
+
+            json_SpatialDiscretization.put("TriangulationType","FullyDistributed");
+            json_SpatialDiscretization.put("MappingX",1);
+            json_SpatialDiscretization.put("MappingV",1);
+            json_root.add_child("SpatialDiscretization",json_SpatialDiscretization);
+
+            json_TemporalDiscretization.put("FinalTime",2);
+            json_TemporalDiscretization.put("CFLNumber",0.15);
+            json_root.add_child("TemporalDiscretization",json_TemporalDiscretization);
+
+            jsout.put("Tick",0.1);
+            json_postprocessing.add_child("StandardOutput",jsout);
+            jVTK.put("Enabled",false);
+            json_postprocessing.add_child("VTK",jVTK);
+
+            json_root.add_child("Postprocessing", json_postprocessing);
+            pt::write_json(newname,json_root);
+
+
+        }
+        catch(pt::ini_parser::ini_parser_error error){
+            std::cout<< "Not a valid File";
+        }
+        catch(...){
+            std::cout << "Something went wrong";
+        }
+    }
+

distributedcombigrid/examples/combi_example/DataConverter.hpp

@@ -0,0 +1,36 @@
+/*
+this function reads in the data from a json file formated like that one from hyperdeal and returns the data 
+we need for a Combi
+
+*/
+
+#ifndef DATACONVERTER_H
+#define DATACONVERTER_H
+
+#include <string>
+#include <vector>
+#include "sgpp/distributedcombigrid/combischeme/CombiMinMaxScheme.hpp"
+
+#include "sgpp/distributedcombigrid/loadmodel/LinearLoadModel.hpp"
+#include "sgpp/distributedcombigrid/manager/CombiParameters.hpp"
+
+#include "sgpp/distributedcombigrid/utils/Types.hpp"
+
+using namespace combigrid;
+
+class Converter{
+    public:
+
+    Converter();
+    //reads the data from a json file and stores it in param file
+
+    void toParam(std::string, std::string);
+
+    //reads in data from an ini file and stores it in  newname.json file
+    void toJSON(std::string, std::string);
+
+    //
+    void toJSON(std::string, std::string, LevelVector );
+
+};
+#endif

distributedcombigrid/examples/combi_example/Make_copy

@@ -0,0 +1,26 @@
+CC=mpicxx
+CFLAGS=-std=c++14 -g -fopenmp -Wno-deprecated-declarations -Wno-unused-local-typedefs -Wno-deprecated -Wno-uninitialized -Wall -DUNIFORMDECOMPOSITION
+
+SGPP_DIR=/home/sgsscratch/goehrims/combi
+GLPK_DIR=/home/sgsscratch/goehrims/combi/glpk
+DEAL_II_DIR=/home/sgsscratch/goehrims/matrixfree-build
+DEAL_II_INC_DIR=/home/sgsscratch/goehrims/matrixfree
+LD_DEALLII=-L$(DEAL_II_DIR)/lib
+INC_DEALII=-I$(DEAL_II_INC_DIR)/include -I$(DEAL_II_DIR)/include
+LD_SGPP=-L$(SGPP_DIR)/lib/sgpp
+INC_GLPK=-I$(GLPK_DIR)/include
+LD_GLPK=-L$(GLPK_DIR)/lib
+INC_SGPP=-I$(SGPP_DIR)/distributedcombigrid/src/ 
+
+LDIR=$(LD_SGPP) $(LD_GLPK) $(LD_DEALII)
+INC=$(INC_SGPP) $(INC_GLPK) $(INC_DEALII)
+
+LIBS=-lsgppdistributedcombigrid -lboost_serialization
+
+all: combi_example 
+
+combi_example: combi_example.cpp TaskExample.hpp
+	mpicxx $(CFLAGS) $(LDIR) $(INC) -o combi_example combi_example.cpp $(LIBS)
+
+clean:
+	rm -f *.o out/* combi_example

distributedcombigrid/examples/combi_example/Makefile

@@ -0,0 +1,23 @@
+CC=mpicxx
+CFLAGS=-std=c++11 -g -fopenmp -Wno-deprecated-declarations -Wno-unused-local-typedefs -Wno-deprecated -Wno-uninitialized -Wall -DUNIFORMDECOMPOSITION
+
+SGPP_DIR=/home/sgsscratch/goehrims/combi
+GLPK_DIR=/home/sgsscratch/goehrims/combi/glpk
+
+LD_SGPP=-L$(SGPP_DIR)/lib/sgpp
+INC_GLPK=-I$(GLPK_DIR)/include
+LD_GLPK=-L$(GLPK_DIR)/lib
+INC_SGPP=-I$(SGPP_DIR)/distributedcombigrid/src/ 
+
+LDIR=$(LD_SGPP) $(LD_GLPK)
+INC=$(INC_SGPP) $(INC_GLPK)
+
+LIBS=-lsgppdistributedcombigrid -lboost_serialization
+
+all: combi_example 
+
+combi_example: combi_example.cpp TaskExample.hpp
+	mpicxx $(CFLAGS) $(LDIR) $(INC) -o combi_example combi_example.cpp $(LIBS)
+
+clean:
+	rm -f *.o out/* combi_example