code: def test_hessian_matrix_order(use_gaussian_derivatives):
      square = np.zeros((5, 5), dtype=float)
      square[2, 2] = 4
  
      Hxx, Hxy, Hyy = hessian_matrix(
          square,
          sigma=0.1,
          order="xy",
          use_gaussian_derivatives=use_gaussian_derivatives,
      )
  
      Hrr, Hrc, Hcc = hessian_matrix(
          square,
          sigma=0.1,
          order="rc",
          use_gaussian_derivatives=use_gaussian_derivatives,
      )
  
      # verify results are equivalent, just reversed in order
      assert_array_equal(Hxx, Hcc)
      assert_array_equal(Hxy, Hrc)
      assert_array_equal(Hyy, Hrr)
  
imports: import numpy
  import pytest
  from numpy.testing import assert_almost_equal, assert_array_equal, assert_equal
  from skimage import data, draw, img_as_float
  from skimage._shared._warnings import expected_warnings
  from skimage._shared.testing import run_in_parallel
  from skimage._shared.utils import _supported_float_type
  from skimage.color import rgb2gray
  from skimage.feature import corner_fast, corner_foerstner, corner_harris, corner_kitchen_rosenfeld, corner_moravec, corner_orientations, corner_peaks, corner_shi_tomasi, corner_subpix, hessian_matrix, hessian_matrix_det, hessian_matrix_eigvals, peak_local_max, shape_index, structure_tensor, structure_tensor_eigenvalues
  from skimage.morphology import cube, octagon
label_desc: Use image node_id 8 for calling a global function with example usage: test_hessian_matrix_order(use_gaussian_derivatives) without return types