@fabianp : does the implementation of this function coincides with your experience in terms of numerical stability required (I am thinking of your recent experiments).