code: def test_from_square_quat_matrix():
      # Ensure proper norm array broadcasting
      x = np.array(
          [[3, 0, 0, 4], [5, 0, 12, 0], [0, 0, 0, 1], [0, 0, 0, -1]]
      )
      r = Rotation.from_quat(x)
      expected_quat = x / np.array([[5], [13], [1], [1]])
      assert_array_almost_equal(r.as_quat(), expected_quat)
  
imports: from __future__ import division, print_function, absolute_import
  import pytest
  import numpy
  from numpy.testing import assert_equal, assert_array_almost_equal
  from numpy.testing import assert_allclose
  from scipy.spatial.transform import Rotation, Slerp
  from scipy.stats import special_ortho_group
  from itertools import permutations
label_desc: Use image node_id 4 for calling a global function with example usage: test_from_square_quat_matrix() without return types