Error in quaternion?

Hi,

Unless I'm doing something wrong, it seems there's an error in the Euler angles returned by a jacobi.Frame created from a quaternion.
Here is a minimal code to reproduce the error where the Euler angles are compared with the ones computed by scipy (scipy also agrees with the results of our own internal library):

import numpy as np
from scipy.spatial.transform import Rotation
import jacobi
for k in range(10):
    q_xyzw = np.array([2 * np.random.rand() - 1, 2 * np.random.rand() - 1, 2 * np.random.rand() - 1, 2 * np.random.rand() - 1])
    q_xyzw = q_xyzw / np.linalg.norm(q_xyzw)
    print("")
    print("jacobi euler: ", *[f"{e:+2.5f}" for e in jacobi.Frame.from_quaternion(0, 0, 0, q_xyzw[3], q_xyzw[0], q_xyzw[1], q_xyzw[2]).euler[3:]])
    print(" scipy euler: ", *[f"{e:+2.5f}" for e in Rotation.from_quat(q_xyzw).as_euler("xyz")])