The test function looks naive, it is not so. Ordinary Simpson 7½ D. Same test function in the transformed domain boosts accuracy to 14½ D.
That would greatly enhance Smath vs competition. I spare many hours trying to code Romberg, what a f...n mess. Rather deceiving, though sometimes as good as Adaptive RemToDo if I can convince myself !