BUG: catch when pdomains was not set by caller
FIN: reproduce LPA of Mori system for fixed total amounts and fixed domain length