--- trunk/OOPSE/libmdtools/Utility.cpp 2004/02/07 00:16:47 1039 +++ trunk/OOPSE/libmdtools/Utility.cpp 2004/05/12 21:54:14 1175 @@ -1,26 +1,24 @@ -#include "Utility.hpp" -#include -#include -#include -using namespace std; - -double dotProduct(vector& v1, vector& v2){ - double sum; - - sum = 0; - - if(v1.size() != v2.size()){ - cerr << "Utility Error: dimension of two vectors are not matched" << endl; - exit(-1); - } - - for(int i = 0; i < v1.size(); i++) - sum += v1[i]*v2[i]; - return sum; -} - -double copysign(double sign, double value){ - double result = fabs(value); - - return sign >= 0 ? result : -result; -} +#include "Utility.hpp" +#include +#include +#include +using namespace std; + +double dotProduct(vector& v1, vector& v2){ + double sum; + + sum = 0; + + if(v1.size() != v2.size()){ + cerr << "Utility Error: dimension of two vectors are not matched" << endl; + exit(-1); + } + + for(int i = 0; i < v1.size(); i++) + sum += v1[i]*v2[i]; + return sum; +} + +double norm2(vector& x){ + return sqrt(dotProduct(x, x)); +}