| # | Line 199 | Line 199 | void Thermo::velocitize() { | |
|---|---|---|
| 199 | vdrift[1] += atoms[vd]->get_vy() * atoms[vd]->getMass(); | |
| 200 | vdrift[2] += atoms[vd]->get_vz() * atoms[vd]->getMass(); | |
| 201 | ||
| 202 | < | mtot = mtot + atoms[vd]->getMass(); |
| 202 | > | mtot += atoms[vd]->getMass(); |
| 203 | } | |
| 204 | ||
| 205 | for (vd = 0; vd < 3; vd++) { | |
| 206 | vdrift[vd] = vdrift[vd] / mtot; | |
| 207 | } | |
| 208 | ||
| 209 | + | |
| 210 | for(vd = 0; vd < n_atoms; vd++){ | |
| 211 | ||
| 212 | vx = atoms[vd]->get_vx(); | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |