Since we are short of time let us just make thin wrappers around qutip for these issues.

rand_herm_ndarray[source]

rand_herm_ndarray(ndarray_shape)

Make a random hermitian ndarray of ndarray_shape

rand_herm_ndarray((2, 2))
array([[0.        +0.j        , 0.21128495-0.36823157j],
       [0.21128495+0.36823157j, 0.197371  +0.j        ]])

test_herm_ndarray[source]

test_herm_ndarray(ndarray)

Returns true if array is Hermitian

a = rand_herm_ndarray((9, 9))
test_herm_ndarray(a)
True