Index: god.cc =================================================================== RCS file: /project/e3101f1b/newcvs/ns/god.cc,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -r1.1.1.1 -r1.2 471,473c471,473 < vector a(mb_node[i]->X(), mb_node[i]->Y(), mb_node[i]->Z()); < vector b(mb_node[j]->X(), mb_node[j]->Y(), mb_node[j]->Z()); < vector d = a - b; --- > Vector a(mb_node[i]->X(), mb_node[i]->Y(), mb_node[i]->Z()); > Vector b(mb_node[j]->X(), mb_node[j]->Y(), mb_node[j]->Z()); > Vector d = a - b; 806,807d805 < < printf("num_nodes is set %d\n", num_nodes); Index: god.h =================================================================== RCS file: /project/e3101f1b/newcvs/ns/god.h,v retrieving revision 1.1.1.1 retrieving revision 1.2 diff -r1.1.1.1 -r1.2 79c79 < class vector { --- > class Vector { 81c81 < vector(double x = 0.0, double y = 0.0, double z = 0.0) { --- > Vector(double x = 0.0, double y = 0.0, double z = 0.0) { 88c88 < inline void vector::operator=(const vector a) { --- > inline void Vector::operator=(const Vector a) { 93c93 < inline void vector::operator+=(const vector a) { --- > inline void Vector::operator+=(const Vector a) { 98c98 < inline int vector::operator==(const vector a) { --- > inline int Vector::operator==(const Vector a) { 101c101 < inline int vector::operator!=(const vector a) { --- > inline int Vector::operator!=(const Vector a) { 104,105c104,105 < inline vector operator-(const vector a) { < return vector(X-a.X, Y-a.Y, Z-a.Z); --- > inline Vector operator-(const Vector a) { > return Vector(X-a.X, Y-a.Y, Z-a.Z); 107,108c107,108 < friend inline vector operator*(const double a, const vector b) { < return vector(a*b.X, a*b.Y, a*b.Z); --- > friend inline Vector operator*(const double a, const Vector b) { > return Vector(a*b.X, a*b.Y, a*b.Z); 110,111c110,111 < friend inline vector operator/(const vector a, const double b) { < return vector(a.X/b, a.Y/b, a.Z/b); --- > friend inline Vector operator/(const Vector a, const double b) { > return Vector(a.X/b, a.Y/b, a.Z/b);