I already asked this question on math.stackexchange but got no answer.

For a non-integer number of degrees of freedom $\nu > p-1$, one can simulate the central Wishart distribution $W_p(\nu, \Sigma)$ with the help of the Bartlett decomposition.

How to simulate the noncentral Wishart distribution $W_p(\nu, \Sigma, \Theta)$ for such a $\nu$ (and $\Theta$ a $p\times p$ semi-definite positive matrix)? According to Letac & Massam's tutorial on non central Wishart distributions, this distribution exists for $\Sigma=I_p$, in the sense that there is a distribution whose characteristic function is the same as the one of $W_{p}(``\text{an integer} \geq p``, I_p, \Theta)$ with $\nu$ instead of the integer. For $\Sigma \neq I_p$ I don't know actually.

For $\nu>2p-1$ it suffices to apply the equality $W_p(\nu, \Sigma, \Theta) = W_{p}(\nu-p, \Sigma) \ast W_{p}(p, \Sigma, \Theta)$. How to do in the case when $p-1 < \nu \leq 2p-1$?