Как Numpy / Scipy превращает функции C в векторизованные функции Python?

Насколько я понимаю, векторизованные функции numpy выполняются быстрее, чем циклы Python, потому что циклы выполняются в C или Fortran. Я хотел бы знать, где в исходном коде это происходит.

Например, биноминальная функция CDF scipy.special.bdtr принимает аргументы в виде массива k,n,p и будет возвращать символ ndarray при условии, что аргументы транслируемы. В документации сказано, что scipy.special.bdtr - это оболочка для подпрограммы в библиотеке математических функций Cephes. Копаясь в исходном коде на Github, я нашел scipy/special/cephes/bdtr.c содержащий код C для этой подпрограммы; вот что я считаю первыми тремя строками соответствующей функции C:

двойной бдтр (k, n, p)

int k, n;

двойной р;

Похоже, что основная функция C не работает с массивами, и я не могу найти исходный код, где эта функция преобразуется в функцию Python, которая работает с массивами.