I have looked into this, but haven't done anything. A good candidate for a library is the cephes package from netlib, because it has many special functions for different numerical types. But to be really useful one would build something like ufunc. This would be fast and independent of dimensions. __Janko