| # | 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 |