Численное интегрирование в Python с адаптивной квадратурой векторизованной функции

Я ищу супер-пупер числовую квадратурную функцию. Он должен иметь следующие три свойства:

• Адаптивный - он автоматически регулирует плотность точек выборки для соответствия подынтегральному выражению. Это абсолютно необходимо, потому что мое подынтегральное выражение очень неоднородно и дорого для вычисления.
• Векторизованное - оно вызывает подынтегральное выражение в списках точек выборки, а не по одной точке за раз, для эффективности.
• Возможность обрабатывать вектор-функции - все компоненты векторной подынтегральной функции вычисляются одновременно без дополнительных затрат, поэтому нет смысла интегрировать все компоненты по отдельности.

Кроме того, это должно быть:

• 2D - интеграл I Я хочу вычислить двойной интеграл по плоской области, и я хочу иметь возможность указать общий (относительный) допуск для всего интеграла и правильно управлять бюджетом ошибок.

Кто-нибудь знает о библиотеке, в которой есть такая функция? Даже два или три из четырех свойств были бы лучше, чем ничего.

Я использую Python и SciPy, поэтому, если он уже работает с Python, это бонус. (Но я также могу написать склеивающий код, чтобы он мог вызывать мое интегрируемое выражение, если это необходимо.)