code: def gemm() -> None:
      A = T.alloc_buffer([16, 16], "float32")
      B = T.alloc_buffer([16, 16], "float32")
      C = T.alloc_buffer([16, 16], "float32")
      for i, j, k, ii, jj in T.grid(4, 4, 16, 4, 4):
          with T.block("update"):
              vi = T.axis.S(16, i * 4 + ii)
              vj = T.axis.S(16, j * 4 + jj)
              vk = T.axis.R(16, k)
              T.reads(A[vi, vk], B[vj, vk])
              T.writes(C[vi, vj])
              with T.init():
                  C[vi, vj] = 0
              C[vi, vj] += A[vi, vk] * B[vj, vk]
  
imports: import pytest
  import tvm
  from tvm import tir
  from tvm.script import tir
  from tvm.ir import Range
label_desc: Use image node_id 8 for calling a global function with example usage: gemm() without return types