This project contains code to accompany the paper [1]. The code is for computing quadrature based approximants to operator convex and operator monotone functions. This is done using an implementation of the weighted differential algorithm, using adaptive barycentric representations of rational functions inspired by [2].
See the example.ipynb notebook for example usage.
[1] O. Faust, H. Fawzi, Rational approximations of operator monotone and operator convex functions, preprint (2023)
[2] S. Filip, Y. Nakatsukasa, L. N. Trefethen, B. Beckermann, Rational minimax approximation via adaptive barycentric representations, SIAM Journal on Scientific Computing (2018)