If so, then TVM won’t call CBLAS functions because it doesn’t know the compute function is actually doing GEMM. You need to call the dense operator in Relay/TOPI:
from tvm import relay
a = relay.var('a', shape=(10, 5))
b = relay.var('b', shape=(20, 5))
c = relay.nn.dense(a, b) # shape of c is (10, 20)
Note that the dense op was implemented to perform A(M, K) x B(N, K) instead of normal A(M, K) x B(K, N) computation.