Skip to content

Commit b24451a

Browse files
authoredJul 13, 2020
Flip sign of non-diagonal intertia tensor elements
A minus sign for non-diagonal elements was missing.
1 parent 1590ba1 commit b24451a

File tree

1 file changed

+5
-5
lines changed

1 file changed

+5
-5
lines changed
 

‎Project#01/hints/step8-solution.md

+5-5
Original file line numberDiff line numberDiff line change
@@ -94,9 +94,9 @@ int main()
9494
I(0,0) += mi * (mol.geom[i][1]*mol.geom[i][1] + mol.geom[i][2]*mol.geom[i][2]);
9595
I(1,1) += mi * (mol.geom[i][0]*mol.geom[i][0] + mol.geom[i][2]*mol.geom[i][2]);
9696
I(2,2) += mi * (mol.geom[i][0]*mol.geom[i][0] + mol.geom[i][1]*mol.geom[i][1]);
97-
I(0,1) += mi * mol.geom[i][0]*mol.geom[i][1];
98-
I(0,2) += mi * mol.geom[i][0]*mol.geom[i][2];
99-
I(1,2) += mi * mol.geom[i][1]*mol.geom[i][2];
97+
I(0,1) -= mi * mol.geom[i][0]*mol.geom[i][1];
98+
I(0,2) -= mi * mol.geom[i][0]*mol.geom[i][2];
99+
I(1,2) -= mi * mol.geom[i][1]*mol.geom[i][2];
100100
}
101101

102102
I(1,0) = I(0,1);
@@ -229,9 +229,9 @@ Torsional angles:
229229
Molecular center of mass: 0.64494926 0.00000000 2.31663792
230230
231231
Moment of inertia tensor (amu bohr^2):
232-
156.154091561645 0.000000000000 52.855584120568
232+
156.154091561645 0.000000000000 -52.855584120568
233233
0.000000000000 199.371126996236 0.000000000000
234-
52.855584120568 0.000000000000 54.459548882464
234+
-52.855584120568 0.000000000000 54.459548882464
235235
236236
Principal moments of inertia (amu * bohr^2):
237237
31.964078 178.649562 199.371127

0 commit comments

Comments
 (0)
Please sign in to comment.