COMPUTE_FORCE
real4 force1 = (real4) (-dEdX, -dEdY, -dEdZ, 0);