This is the same as the one in el.c, but uses native D reals instead of the soft long double ones.
See Implementation
This is the same as the one in el.c, but uses native D reals instead of the soft long double ones.