code: def nested_block_bind(
      A: T.Buffer((16, 16, 16, 16), "float32"),
      B: T.Buffer((16, 16, 16), "float32"),
  ):
      for i, j in T.grid(16, 16):
          with T.block("outer"):
              vi, vj = T.axis.remap("SS", [i, j])
              for k, l in T.grid(16, 16):
                  with T.block("inner"):
                      vk, vl = T.axis.remap("SR", [k, l])
                      with T.init():
                          B[vi, vj, vk] = 0.0
                      B[vi, vj, vk] = B[vi, vj, vk] + A[vi, vj, vk, vl]
  
imports: import sys
  import pytest
  import tvm
  import tvm.testing
  from tvm import tir
  from tvm.script import tir
  from tvm.tir.schedule.testing import verify_trace_roundtrip, assert_structural_equal_ignore_global_symbol
label_desc: Use image node_id 21 for calling a global function with example usage: nested_block_bind(A, B) without return types