Analytic functions with blackbody

See #3077 
Yikes! Maybe a combination with single-precision and not having `expm1`, but still...