```py q_B_ECI = Quaternion(np.array((0.0, 0.0, 0.0, 1.0)), Quaternion.Format.XYZS) ``` returns a `[0.0, 0.0, 0.0, 0.0]` quaternion. Seems to be due to a faulty `np.array -> Vector4d` conversion. Workaround is to use: ```py q_B_ECI = Quaternion.XYZS(0.0, 0.0, 0.0, 1.0) ``` instead.