Obtain one uniform draw from d times d Hpd matrices with unit trace See Algorithm 2 in Mittelbach et al. (adjusted to complex case)
unit_trace_runif_single(d)