Custom Distribution - Multivariate SkewNormal

I believe the summation is handled my PyMC.