Performance regression for transitive closure on large `rel[loc,loc]`

[jurgenvinju]

Large transitive closures behave quadratically, as it seems.
We traced this down to the code level and saw that the default method of IRelation.closure is simply calling compose iteratively until a fixed point.