Blackbody radiation curves, with support for Unitful.jl.
TODO (help welcome, not expected)
- Documentation
- Need at least one example showing easy calculations/plotting
- need another example showing fits/gradients
- Unit tests
- Gradients
- Write out some chain rules (ChainRulesCore.jl)
Using Optim.jl we can numerically maximize the blackbody function to see how closely it agrees with Wien's law.
using Optim
# function to minimize
Teff = 5796
f(λ) = -blackbody(λ, Teff)
# bounded between 1 picometer and 10 m
res = Optim.optimize(f, 1e-12, 10)
λ = Optim.minimizer(res)
isapprox(λ, 2.898e-3 / Teff; rtol=1e-4)If you would like to contribute, feel free to open a pull request. If you want to discuss something before contributing, head over to discussions and join or open a new topic. If you're having problems with something, please open an issue.