“…We emphasize that this algorithm assumes that the initial map is positive in the sense that its Jacobian is positive, that is, ∇ x h ≥ 0. To this purpose, a reordering algorithm is required (and is described in [18]). In particular, the operator (∇ x ) −1 is presented in [16].…”