OOPSE-1.0/libmdtools/Utility.cpp

#include "Utility.hpp"
#include <vector>
#include <iostream>
#include <math.h>
using namespace std;

double dotProduct(vector<double>& v1, vector<double>& 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<double>& x){
  return sqrt(dotProduct(x, x));
}       
Revision:	1334
Committed:	Fri Jul 16 18:58:03 2004 UTC (19 years, 11 months ago) by gezelter
File size:	468 byte(s)
Log Message:	Initial import of OOPSE-1.0 source tree
#	User	Rev	Content
1	gezelter	1334	#include "Utility.hpp"
2			#include <vector>
3			#include <iostream>
4			#include <math.h>
5			using namespace std;
6
7			double dotProduct(vector<double>& v1, vector<double>& v2){
8			double sum;
9
10			sum = 0;
11
12			if(v1.size() != v2.size()){
13			cerr << "Utility Error: dimension of two vectors are not matched" << endl;
14			exit(-1);
15			}
16
17			for(int i = 0; i < v1.size(); i++)
18			sum += v1[i]*v2[i];
19			return sum;
20			}
21
22			double norm2(vector<double>& x){
23			return sqrt(dotProduct(x, x));
24			}