# | Line 59 | Line 59 | double Thermo::getKinetic(){ | |
---|---|---|
59 | k = (i+2)%3; | |
60 | kinetic += aJ[j]*aJ[j]/I[j][j] + aJ[k]*aJ[k]/I[k][k]; | |
61 | } else { | |
62 | < | for (j=0; j<3; j++) |
63 | < | kinetic += aJ[j]*aJ[j] / I[j][j]; |
62 | > | for (j=0; j<3; j++) |
63 | > | kinetic += aJ[j]*aJ[j] / I[j][j]; |
64 | } | |
65 | < | } |
65 | > | } |
66 | } | |
67 | #ifdef IS_MPI | |
68 | MPI_Allreduce(&kinetic,&kinetic_global,1,MPI_DOUBLE, | |
69 | MPI_SUM, MPI_COMM_WORLD); | |
70 | kinetic = kinetic_global; | |
71 | #endif //is_mpi | |
72 | < | |
72 | > | |
73 | kinetic = kinetic * 0.5 / e_convert; | |
74 | ||
75 | return kinetic; |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |