From 48e34132dfef67ab3634c680b3121b1a02fc2827 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= <sebastien@dynare.org>
Date: Fri, 9 May 2025 14:11:52 +0200
Subject: [PATCH] getPowerDeriv: for the derivative of x^0 in 0, return 0
instead of NaN
Ref. #73 and dynare!2358
---
src/DataTree.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/DataTree.cc b/src/DataTree.cc
index c17caa9e..fd2ffc09 100644
--- a/src/DataTree.cc
+++ b/src/DataTree.cc
@@ -1,5 +1,5 @@
/*
- * Copyright © 2003-2024 Dynare Team
+ * Copyright © 2003-2025 Dynare Team
*
* This file is part of Dynare.
*
@@ -912,7 +912,7 @@ DataTree::writeCHelpersDefinition(ostream& output) const
<< "getPowerDeriv(double x, double p, int k)" << endl
<< "{" << endl
<< " if (fabs(x) < " << power_deriv_near_zero
- << " && p > 0 && k > p && fabs(p-nearbyint(p)) < " << power_deriv_near_zero << ')'
+ << " && p >= 0 && k > p && fabs(p-nearbyint(p)) < " << power_deriv_near_zero << ')'
<< endl
<< " return 0.0;" << endl
<< " else" << endl
--
GitLab