def METHOD_NAME():
    nrow = 3
    ncol = 4
    mat = h.Matrix(nrow, ncol)
    h("double a[%d][%d]" % (nrow, ncol))
    for i in range(nrow):
        for j in range(ncol):
            mat.x[i][j] = i * ncol + j
            h.a[i][j] = i * ncol + j
    rmat = mat._ref_x[0][0]
    ra = h._ref_a[0][0]
    for i in range(nrow * ncol):
        assert rmat[i] == float(i)
        assert ra[i] == float(i)