# For toeplitz matrices, solve_toeplitz() should be equivalent to solve().
When dim>256, mat result is wrong. mat[j]=DeviceArray([[1., 1., 1., ..., 1., 1., 1.], [1., 1., 1., ..., 1., 1., 1.], [1., 1., 1., ..., 1., 1., 1.], ..., [1., 1., 1 ...