Fix function find_direct_relation. Use a slice based on a chunk size instead of selecting a single time-step.
Check that the constrains are fulfilled after a binary search.