code: def test_orthogonalize(rng):
      """Test that the orthogonalization is OK."""
      X = rng.standard_normal(size=(100, 5))
  
      X = orthogonalize(X)
  
      K = np.dot(X.T, X)
      K -= np.diag(np.diag(K))
  
      assert_almost_equal((K**2).sum(), 0, 15)
  
imports: import warnings
  import numpy
  import pytest
  from numpy.testing import assert_almost_equal, assert_array_almost_equal, assert_array_equal
  from nilearn.glm.first_level.hemodynamic_models import _calculate_tr, _hrf_kernel, _regressor_names, _resample_regressor, _sample_condition, compute_regressor, glover_dispersion_derivative, glover_hrf, glover_time_derivative, orthogonalize, spm_dispersion_derivative, spm_hrf, spm_time_derivative
label_desc: Use image node_id 6 for calling a global function with example usage: test_orthogonalize(rng) without return types