Add a calculate_constraints function that uses IPF

Estimates the total number of each household or person type required to match marginal totals according to observed sample proportions.