Function uk_money_dirty_from_snapshot 

fn uk_money_dirty_from_snapshot(snap: &UkDmoSnapshot, rho: f64) -> f64

The DMO case formula at a trial real yield:

P = [d₁ + d₂(uw) + a·c·w²/(2(1−w))·(1 − w^{n−1})]·(uw)^{r/s}
    + 100·a·u^{r/s}·w^{r/s+n}

with w = 1/(1+ρ/2) and u = (1+π)^{−1/2}.