change

Solver/eigenSolver.cpp

@@ -184,28 +184,28 @@ bool eigenSolver::solve(int numEigenValues, std::string which)
 }
 
 
-void eigenSolver::normalize_mode(){
+void eigenSolver::normalize_mode(double scale){
 	Msg::Info("Normalize all eigenvectors");
   for (unsigned int i=0; i<_eigenVectors.size(); i++){
 		double norm = 0.;
     for (unsigned int j=0; j<_eigenVectors[i].size(); j++){
 			std::complex<double> val = _eigenVectors[i][j];
       double normval = std::abs(val);
-      if (normval >norm) 
+      if (normval >norm)
 				norm = normval;
     };
-    
-	  
+
+
     if (norm == 0) {
 			Msg::Error("zero eigenvector");
       return;
 		};
-    	
+
  		for (unsigned int j=0; j<_eigenVectors[i].size(); j++){
-    	_eigenVectors[i][j] /= norm;
+    	_eigenVectors[i][j] *= (scale/norm);
     };
 
-  };  
+  };
 };
 
 #endif

Solver/eigenSolver.h

@@ -42,7 +42,7 @@ class eigenSolver{
     return _eigenVectors[num][com];
   };
   int getNumberEigenvectors() {return _eigenVectors.size();};
-  void normalize_mode();
+  void normalize_mode(double scale=1.);
 };
 
 #else
@@ -68,7 +68,7 @@ class eigenSolver{
   void clear(){}
   std::complex<double> getEigenVectorComp(int num, int com) { return 0.; }
 	int getNumberEigenvectors() {return 0;};
-	void normalize_mode() {};
+	void normalize_mode(double scale=1.) {};
 };
 
 #endif