Get log(d) vector, see (39) in Mittelbach et al, adjusted to complex case Helping function for unit_trace_runif
unit_trace_runif
unit_trace_log_d(p, q)
Useful links