42 |
|
grad[index++] = -dAtomGrad[0]; |
43 |
|
grad[index++] = -dAtomGrad[1]; |
44 |
|
grad[index++] = -dAtomGrad[2]; |
45 |
< |
grad[index++] = dAtomGrad[3]; |
46 |
< |
grad[index++] = dAtomGrad[4]; |
47 |
< |
grad[index++] = dAtomGrad[5]; |
45 |
> |
grad[index++] = -dAtomGrad[3]; |
46 |
> |
grad[index++] = -dAtomGrad[4]; |
47 |
> |
grad[index++] = -dAtomGrad[5]; |
48 |
|
|
49 |
|
} |
50 |
|
else{ |
128 |
|
x[index++] = eulerAngle[0]; |
129 |
|
x[index++] = eulerAngle[1]; |
130 |
|
x[index++] = eulerAngle[2]; |
131 |
< |
|
131 |
> |
|
132 |
|
} |
133 |
|
|
134 |
|
} |
158 |
|
calcForce(1, 1); |
159 |
|
dumpOut->writeDump(iteration); |
160 |
|
statOut->writeStat(iteration); |
161 |
< |
} |
161 |
> |
} |