Since we are short of time let us just make thin wrappers around qutip for these issues.
rand_herm_ndarray
rand_herm_ndarray(ndarray_shape)
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
test_herm_ndarray(ndarray)
ndarray
Returns true if array is Hermitian
a = rand_herm_ndarray((9, 9)) test_herm_ndarray(a)
True