Function erfc_inv

extlibs/catch/include/catch/catch.hpp:7449–7451

            }

            double erfc_inv(double x) {
                return erf_inv(1.0 - x);
            }

            double normal_quantile(double p) {
                static const double ROOT_TWO = std::sqrt(2.0);