Hi,

Here the code I use based on this link :

void CQuaternion::ToEulerAngles( float* X, float* Y, float* Z ) const
{
  const float SingularityTest = ( x * y ) + ( z * w );
  if( SingularityTest > 0.499f )
  {
    *Y = 2.0f * CMath::ATan2( x, w );
    *Z = CMath::HALF_PI;
    *X = 0.0f;
  }
  else if( SingularityTest < -0.499f )
  {
    *Y = -2.0f * CMath::ATan2( x, w );
    *Z = -CMath::HALF_PI;
    *X = 0.0f;
  }
  else
  {
    const float zz = z * z;
    *Y = CMath::ATan2( 2.0f * ( ( y * w ) - ( x * z ) ), 1.0f - 2.0f * ( ( y * y ) - zz ) );
    *Z = CMath::ASin( 2.0f * SingularityTest );
    *X = CMath::ATan2( 2.0f * ( ( x * w ) - ( y * z ) ), 1.0f - 2.0f * ( ( x * x ) - zz ) );
  }
}

Quaternion To Euler is called when gizmo is used to shows the new rotation on the rotation property in the editor.

It's also called when you attach one actor to another since the rotation is transformed to local space.

The problem is it can fastly gives bad result.

There is a fix to have correct result ?

Thanks

Nice catch ! That works better but it's not very accurate, after a quat to euler if the euler is modified the rotation move slightly.

In some case I have really bad euler resulting from this cacule.

Here the matrix version which ends to the same issue :

void CQuaternion::ToEulerAngles( float* X, float* Y, float* Z ) const
{
  // Convert to matrix.
  const CMatrix4 Matrix = ToMatrix();

  // Extract eulers.
  if( Matrix.m44[ 1 ][ 0 ] > 0.998f )
  {
    *Y = CMath::ATan2( Matrix.m44[ 0 ][ 2 ], Matrix.m44[ 2 ][ 2 ] );
    *Z = CMath::HALF_PI;
    *X = 0.0f;
  }
  else if( Matrix.m44[ 1 ][ 0 ] < -0.998f )
  {
    *Y = CMath::ATan2( Matrix.m44[ 0 ][ 2 ], Matrix.m44[ 2 ][ 2 ] );
    *Z = -CMath::HALF_PI;
    *X = 0.0f;
  }
  else
  {
    *Y = CMath::ATan2( -Matrix.m44[ 2 ][ 0 ], Matrix.m44[ 0 ][ 0 ] );
    *Z = CMath::ASin( Matrix.m44[ 1 ][ 0 ] );
    *X = CMath::ATan2( -Matrix.m44[ 1 ][ 2 ], Matrix.m44[ 1 ][ 1 ] );
  }
}

The only way to keep correct euler angles is to track them and compute the quaternion from them apparently.

sVec3 ToEulerZYX (const int order = 0x012) 
	{
		int a0 = (order>>8)&3;
		int a1 = (order>>4)&3;
		int a2 = (order>>0)&3;

		sVec3 euler;
		// Assuming the angles are in radians.
		if ((*this)[a0][a2] > (1.0-FP_EPSILON)) 
		{ // singularity at north pole
			euler[a0] = -atan2((*this)[a2][a1], (*this)[a1][a1]);
			euler[a1] = -3.1415926535897932384626433832795f/2.0f;
			euler[a2] = 0;
			return euler;
		}
		if ((*this)[a0][a2] < -(1.0-FP_EPSILON))
		{ // singularity at south pole
			euler[a0] = -atan2((*this)[a2][a1], (*this)[a1][a1]);
			euler[a1] = 3.1415926535897932384626433832795f/2.0f;
			euler[a2] = 0;
			return euler;
		}
		euler[a0] =	-atan2(-(*this)[a1][a2], (*this)[a2][a2]);
		euler[a1] =	-asin ( (*this)[a0][a2]);
		euler[a2] =	-atan2(-(*this)[a0][a1], (*this)[a0][a0]);
		return euler;
	}

		sQuat quat (-0.700637, 0.361543, -0.412074, 0.456716);
		quat.Normalize();
		sMat3 matrix;
		matrix.FromQuat (quat);

		matrix = matrix.Transposed(); // column / major issue ?
		
		sVec3 radians = matrix.ToEuler(0x120);
		sVec3 verify	= radians * 180.0f/PI; // {mX=133.219284 mY=31.7808933 mZ=62.0091209 ...} 

		sMat3 matrixX = sMat3::rotationX (radians[0]);
		sMat3 matrixY = sMat3::rotationY (radians[1]);
		sMat3 matrixZ = sMat3::rotationZ (radians[2]);
		
		sMat3 verifyMatrix = matrixX * matrixZ * matrixY; // == matrix

		sQuat quat (-0.700637, 0.361543, -0.412074, 0.456716);
		quat.Normalize();
		sMat3 matrix;
		matrix.FromQuat (quat);
		//matrix = matrix.Transposed();
		
		sVec3 radians = matrix.ToEuler(0x021);		
		sVec3 verify = radians * 180.0f/PI; // {mX=133.219284 mY=31.7808933 mZ=62.0091209 ...} 

		sMat3 matrixX = sMat3::rotationX (radians[0]);
		sMat3 matrixY = sMat3::rotationY (radians[1]);
		sMat3 matrixZ = sMat3::rotationZ (radians[2]);
		
		sMat3 verifyMatrix = matrixX * matrixZ * matrixY; // == matrix