investigate why pij is not respected in case2

[samthiriot]

data(dwellings_households)
prepared <- matching.prepare(dwellings_households$sample.A, dwellings_households$sample.B, dwellings_households$pdi, dwellings_households$pdj, dwellings_households$pij)
solved <- matching.solve(prepared, nA=50000, nB=40000, nu.A=1, phi.A=1, delta.A=1, gamma=0, delta.B=1, phi.B=1, nu.B=1)
generated <- matching.generate(solved, dwellings_households$sample.A, dwellings_households$sample.B)
plot(generated, dwellings_households$sample.A$sample, dwellings_households$sample.B$sample)