diff --git a/libs/math/include/psemek/math/dual.hpp b/libs/math/include/psemek/math/dual.hpp index dc198927..6de9473e 100644 --- a/libs/math/include/psemek/math/dual.hpp +++ b/libs/math/include/psemek/math/dual.hpp @@ -175,6 +175,13 @@ namespace psemek::math return {std::log(d.scalar), d.delta / d.scalar}; } + template + dual sqrt(dual const & d) + { + auto s = std::sqrt(d.scalar); + return {s, d.delta / (T{2} * s)}; + } + template dual sin(dual const & d) {