X-Git-Url: http://git.treefish.org/~alex/phys/latlib.git/blobdiff_plain/a5f05e337d18c193ac57b4da0013d1cdc69f8faa..8d11db1cf57b103fcec701d71c9d620b1c381b73:/paraq.cpp?ds=inline

diff --git a/paraq.cpp b/paraq.cpp
index ac8fc14..0892962 100644
--- a/paraq.cpp
+++ b/paraq.cpp
@@ -41,11 +41,24 @@ string paraq::rangeString()
   for( map< string, vector<double*> >::iterator paraIt = rangeMap.begin(); paraIt != rangeMap.end(); ++paraIt )
     if( paraIt->second.size() > 0 )
       for( vector<double*>::iterator rangeIt = paraIt->second.begin(); rangeIt != paraIt->second.end(); ++rangeIt)
-	rangestring << "-" << paraIt->first << (*rangeIt)[0] << "_" << (*rangeIt)[1] << "_" << (*rangeIt)[2];
+	rangestring << "_" << paraIt->first << (*rangeIt)[0] << "-" << (*rangeIt)[1] << "-" << (*rangeIt)[2];
   
+  for (vector<paralink>::iterator linkit = linkedParas.begin(); linkit != linkedParas.end(); ++linkit)
+    rangestring << "_" << linkit->first << "--" << linkit->linktype << "--" << linkit->second;
+
   for( map<string,double>::iterator defIt = defaultPara.begin(); defIt != defaultPara.end(); ++defIt )
-    if( rangeMap.find(defIt->first) == rangeMap.end() )
-      rangestring << "-" << defIt->first << defIt->second;
+    if( rangeMap.find(defIt->first) == rangeMap.end() ) {
+      bool linkedpara = false;
+
+      for (vector<paralink>::iterator linkit = linkedParas.begin(); linkit != linkedParas.end(); ++linkit)
+	if (linkit->first == defIt->first) {
+	  linkedpara = true;
+	  break;
+	}
+
+      if (!linkedpara)
+	rangestring << "_" << defIt->first << defIt->second;
+    }
 
   return rangestring.str();
 }
@@ -81,7 +94,13 @@ void paraq::initJobList()
       while( jobList.begin()->find(paraIt->first) == jobList.begin()->end() ) {
 	for( vector<double>::iterator valIt = paraIt->second.begin(); valIt != paraIt->second.end(); ++valIt ) {
 	  jobList.push_back( *jobList.begin() );
+
 	  jobList.back()[paraIt->first] = *valIt;
+
+	  /* also set linked parameters */
+	  for (vector<paralink>::iterator linkit = linkedParas.begin(); linkit != linkedParas.end(); ++linkit)
+	    if (linkit->second == paraIt->first)
+	      jobList.back()[linkit->first] = linkTypeFunction(*valIt, linkit->linktype);
 	}	
 	jobList.erase( jobList.begin() );
       }
@@ -89,6 +108,18 @@ void paraq::initJobList()
   }
 }
 
+double paraq::linkTypeFunction (double linkval, int linktype) {
+  switch (linktype) {
+
+  case LINK_1TO1:
+    return linkval;
+
+  case LINK_NEGATIVE:
+    return -linkval;
+
+  }
+}
+
 int paraq::nextParas()
 {
   if( thisjob == -1 ) {   
@@ -114,7 +145,8 @@ double& paraq::operator[] (string paraid) {
     return defaultPara[paraid];
   }
 
-  return masterdefault;
+  cerr << "PARAQ: Parameter " << paraid << " does not exist!" << endl << flush;
+  exit(1);
 }
 
 int paraq::getTotalJobs() {