diff --git a/motions/tetrahedral.cpp b/motions/tetrahedral.cpp index d6fb36d..a23ea72 100644 --- a/motions/tetrahedral.cpp +++ b/motions/tetrahedral.cpp @@ -29,8 +29,8 @@ void TetrahedralJump::initialize() { if (cos_theta != 1 && cos_theta != -1) { // std::cout << cos_theta << std::endl; rotated_position = { - 0*m_cos_beta * vec[0] + m_sin_beta * (-vec[0] * vec[2] * sin_alpha - vec[1] * cos_alpha) / norm, - 0*m_cos_beta * vec[1] + m_sin_beta * (-vec[1] * vec[2] * sin_alpha + vec[0] * cos_alpha) / norm, + m_cos_beta * vec[0] + m_sin_beta * (-vec[0] * vec[2] * sin_alpha - vec[1] * cos_alpha) / norm, + m_cos_beta * vec[1] + m_sin_beta * (-vec[1] * vec[2] * sin_alpha + vec[0] * cos_alpha) / norm, m_cos_beta * vec[2] + m_sin_beta * norm * sin_alpha }; } else {