Constrain template member functions in preparemaps.h

Author: contuandrea

include/maps/preparemaps.h

@@ -55,6 +55,8 @@ namespace maps{
         bool _gridyfied=false;
         double _scalespace=1.;
         double _scalefield=1.;
+        bool _has_active_flag=false;
+        std::string _active_flag="";
         //data points container
         hydra::multiarray<MP,3,hydra::host::sys_t> _extremes;
         hydra::multiarray<MP,N,hydra::host::sys_t> _data;
@@ -79,18 +81,22 @@ namespace maps{
         physmap(){
         }
 
-        physmap(std::string s, size_t skipcolumns=0, size_t skiprows=0, double scale=1., double scalespace=1., bool sparse=false){
-            this->init(s,skipcolumns,skiprows,scale,scalespace,sparse);
+        physmap(std::string s, size_t skipcolumns=0, size_t skiprows=0, double scale=1., double scalespace=1., bool sparse=false, std::string search_active=""){
+            this->init(s,skipcolumns,skiprows,scale,scalespace,sparse,search_active);
         }
 
-        void init(std::string s, size_t skipcolumns=0, size_t skiprows=0, double scale=1., double scalespace=1., bool sparse=false){
+        void init(std::string s, size_t skipcolumns=0, size_t skiprows=0, double scale=1., double scalespace=1., bool sparse=false, std::string search_active=""){
             _filename=s;
             _cskip=skipcolumns;
             _rskip=skiprows;
             _sparse=false;
             _gridyfied=false;
             _scalefield=scale;
             _scalespace=scalespace;
+            if(search_active!=""){
+                _has_active_flag = true;
+                _active_flag = search_active;
+            }
         }
 
         void prepare(){
@@ -139,70 +145,96 @@ namespace maps{
             hydra::copy(hydra::make_range(_data.begin(_3),_data.end(_3)),v);
         }
 
+        
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 2 && N==5, void >
+        fillscalar_and_active(Container &v){
+            v.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_3,_4),_data.end(_3,_4)),v);
+        }
+        
         template <typename Container>
-        void fillemobility(Container &v){
+        std::enable_if<N==11 || N==12, void >
+        fillemobility(Container &v){
             v.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_6),_data.end(_6)),v);
         }
 
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 2 && N==12, void >
+        fillemobility_and_active(Container &v){
+            v.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_6,_11),_data.end(_6,_11)),v);
+        }
+        
         template <typename Container>
-        void fillhmobility(Container &v){
+        std::enable_if<N==11 || N==12, void >
+        fillhmobility(Container &v){
             v.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_7),_data.end(_7)),v);
         }
 
-        template <typename Container>
-        void fillvector(Container &vx, Container &vy, Container &vz){
-            vx.resize(_data.size());
-            vy.resize(_data.size());
-            vz.resize(_data.size());
-            hydra::copy(hydra::make_range(_data.begin(_3),_data.end(_3)),vx);
-            hydra::copy(hydra::make_range(_data.begin(_4),_data.end(_4)),vy);
-            hydra::copy(hydra::make_range(_data.begin(_5),_data.end(_5)),vz);
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 2 && N==11, void >
+        fillhmobility_and_active(Container &v){
+            v.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_7,_11),_data.end(_7,_11)),v);
         }
 
-        template <typename Container>
-        void fillvector(Container &vect){
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 3 && (N==6 || N==7), void >
+        fillvector(Container &vect){
             vect.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_3,_4,_5),_data.end(_3,_4,_5)),vect);
         }
 
-        template <typename Container>
-        void fillvectorEfield(Container &vect){
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 4 && N==7, void >
+        fillvector_and_active(Container &vect){
+            vect.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_3,_4,_5,_6),_data.end(_3,_4,_5,_6)),vect);
+        }
+        
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 3 && (N==11 || N==12), void >
+        fillvectorEfield(Container &vect){
             vect.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_3,_4,_5),_data.end(_3,_4,_5)),vect);
         }
 
-        template <typename Container>
-        void fillvectorWfield(Container &vect){
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 4 && N==12, void >
+        fillvectorEfield_and_active(Container &vect){
+            vect.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_3,_4,_5,_11),_data.end(_3,_4,_5,_11)),vect);
+        }
+        
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 3 && (N==11 || N==12), void >
+        fillvectorWfield(Container &vect){
             vect.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_8,_9,_10),_data.end(_8,_9,_10)),vect);
         }
 
-        template <typename Container>
-        void fillallquantities(Container &fields){
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 4 && N==12, void >
+        fillvectorWfield_and_active(Container &vect){
+            vect.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_8,_9,_10,_11),_data.end(_8,_9,_10,_11)),vect);
+        }
+        
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 8 && (N==11 || N==12), void >
+        fillallquantities(Container &fields){
             fields.resize(_data.size());
             hydra::copy(hydra::make_range(_data.begin(_3,_4,_5,_6,_7,_8,_9,_10),_data.end(_3,_4,_5,_6,_7,_8,_9,_10)),fields);
         }
-        
-        template <typename Container>
-        void fillallquantities(Container &ex, Container &ey, Container &ez, Container &emob, Container &hmob, Container &wx, Container &wy, Container &wz){
-            ex.resize(_data.size());
-            ey.resize(_data.size());
-            ez.resize(_data.size());
-            wx.resize(_data.size());
-            wy.resize(_data.size());
-            wz.resize(_data.size());
-            emob.resize(_data.size());
-            hmob.resize(_data.size());
-            hydra::copy(hydra::make_range(_data.begin(_3),_data.end(_3)),ex);
-            hydra::copy(hydra::make_range(_data.begin(_4),_data.end(_4)),ey);
-            hydra::copy(hydra::make_range(_data.begin(_5),_data.end(_5)),ez);
-            hydra::copy(hydra::make_range(_data.begin(_6),_data.end(_6)),emob);
-            hydra::copy(hydra::make_range(_data.begin(_7),_data.end(_7)),hmob);
-            hydra::copy(hydra::make_range(_data.begin(_8),_data.end(_8)),wx);
-            hydra::copy(hydra::make_range(_data.begin(_9),_data.end(_9)),wy);
-            hydra::copy(hydra::make_range(_data.begin(_10),_data.end(_10)),wz);
+            
+        template <typename Container, typename ElemType = decltype(*(std::declval<Container>().begin()))>
+        std::enable_if_t< hydra::tuple_size< ElemType >::value == 9 && N==12, void >
+        fillallquantities_and_active(Container &fields){
+            fields.resize(_data.size());
+            hydra::copy(hydra::make_range(_data.begin(_3,_4,_5,_6,_7,_8,_9,_10,_11),_data.end(_3,_4,_5,_6,_7,_8,_9,_10,_11)),fields);
         }
 
     };
@@ -576,7 +608,6 @@ namespace maps{
         }
 
         void load(std::string map_path){
-            INFO_LINE("Loading single map")
             _pm.init(map_path);
             _pm.prepare();
             std::get<0>(_borders) = *(_pm.getx().begin());
@@ -592,8 +623,6 @@ namespace maps{
         }
 
         void load(std::string efield_path, std::string emob_path, std::string hmob_path, std::string wfield_path){
-            INFO_LINE("Loading separate maps")
-            
             //get efield
             _efieldm.init(efield_path); //E field is a vector
             _efieldm.prepare();