This PR corrects the domain calculations for compose_trans(), which appear to be incorrect.

The problem

Currently, compose_trans() pushes the domain of the first function through all transformations to calculate the domain. This procedure, when it works (i.e. if the range of each transformation in the chain is contained within the domain of the next one), yields the range of the composed set of transformations, not the domain.

This leads to incorrect domains; e.g.:

compose_trans("sqrt", "reverse")$domain
## [1] -Inf    0

The correct domain here is just the domain of sqrt(), which is c(0, Inf). We can see this by plotting it:

plot(compose_trans("sqrt", "reverse"), xlim = c(-1,4))
## Warning message:
## In trans$transform(x) : NaNs produced

This also leads to compose_trans() incorrectly erroring on valid compositions; e.g.:

compose_trans("reverse", "log")
## Error in `compose_trans()`:
## ! Sequence of transformations yields invalid domain
## Run `rlang::last_trace()` to see where the error occurred.

Yet, this should be a valid transformation, just with the domain c(-Inf, 0) instead of the original domain of c(0, Inf) (this also offers a solution to the question about reverse log in #393: the transformation is valid so long as the domain is correctly calculated).

The solution

This PR adjusts the domain calculation to push the first domain forward through the transformations, intersecting it with the domains of each transformation along the way, then pushing the resulting range back through the inverses to get the domain.

On the above two examples, with this PR, we have:

compose_trans("sqrt", "reverse")$domain
## [1]   0 Inf

and

compose_trans("reverse", "log")$domain
## [1]    -Inf -1e-100