Me too.

Bit late (15 years), but the same problem. I need to generate an exponential decay and PBP has no ^/exp function.

There are several ways to do this, but none of them are particularly quick, which is also a complaint of most built-in exp functions. I think, when we get this right, the fastest answer is a lookup table/array approach.

https://nic.schraudolph.org/bib2html...audolph99.html

https://codingforspeed.com/using-fas...approximation/

This doesn't seem too difficult so will try implementing this.

double exp1(double x) {
x = 1.0 + x / 256.0;
x *= x; x *= x; x *= x; x *= x;
x *= x; x *= x; x *= x; x *= x;
return x;
}