Changeset 4ef101 for molecuilder/src/molecule_geometry.cpp
- Timestamp:
- Nov 5, 2009, 6:14:08 PM (16 years ago)
- Children:
- fc420f
- Parents:
- 543ce4
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
molecuilder/src/molecule_geometry.cpp
r543ce4 r4ef101 33 33 ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv); 34 34 35 delete(M);36 delete(Minv);35 Free(&M); 36 Free(&Minv); 37 37 delete(Center); 38 38 return status; … … 52 52 ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv); 53 53 54 delete(M);55 delete(Minv);54 Free(&M); 55 Free(&Minv); 56 56 return status; 57 57 }; … … 233 233 ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv); 234 234 235 delete(M);236 delete(Minv);235 Free(&M); 236 Free(&Minv); 237 237 }; 238 238 … … 253 253 atom *Walker = start; 254 254 double *matrix = ReturnFullMatrixforSymmetric(cell_size); 255 double *inversematrix = InverseMatrix(cell_size); 255 256 double tmp; 256 257 bool flag; … … 266 267 #endif 267 268 Testvector.CopyVector(&Walker->x); 268 Testvector. InverseMatrixMultiplication(matrix);269 Testvector.MatrixMultiplication(inversematrix); 269 270 Translationvector.Zero(); 270 271 for (BondList::const_iterator Runner = Walker->ListOfBonds.begin(); Runner != Walker->ListOfBonds.end(); (++Runner)) { … … 293 294 if ((*Runner)->GetOtherAtom(Walker)->type->Z == 1) { 294 295 Testvector.CopyVector(&(*Runner)->GetOtherAtom(Walker)->x); 295 Testvector. InverseMatrixMultiplication(matrix);296 Testvector.MatrixMultiplication(inversematrix); 296 297 Testvector.AddVector(&Translationvector); 297 298 Testvector.MatrixMultiplication(matrix); … … 307 308 } while (!flag); 308 309 Free(&matrix); 310 Free(&inversematrix); 311 309 312 Center.Scale(1./(double)AtomCount); 310 313 };
Note:
See TracChangeset
for help on using the changeset viewer.