[1]} suggest a constrained HMC algorithm using the Hamiltonian
\(H(\theta , p) = \Big (\mathfrak {K}(\theta ,p)\Big ) + \Big (\mathfrak {U}(\theta )\Big ) := \left( \frac{1}{2} p^{\prime } M(\theta )^{-1}p \right) + \left( \frac{1}{2} \log | M(\theta ) | - \log \pi (\theta ) \right) + \lambda ^{\prime } g(\theta ),\) 
