From 980177dab9be3fdaab9ef9e4b93c059b3ab338cf Mon Sep 17 00:00:00 2001
From: martinjrobins <martinjrobins@gmail.com>
Date: Sat, 30 May 2026 22:52:10 +0000
Subject: [PATCH 1/4] version 0.12.0

---
 diffsl/Cargo.toml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/diffsl/Cargo.toml b/diffsl/Cargo.toml
index aed7ec6..56555f2 100644
--- a/diffsl/Cargo.toml
+++ b/diffsl/Cargo.toml
@@ -1,6 +1,6 @@
 [package]
 name = "diffsl"
-version = "0.11.11"
+version = "0.12.0"
 edition.workspace = true
 license.workspace = true
 authors.workspace = true

From a0412abaf75cde30d873e44b143c9df7f0c51abe Mon Sep 17 00:00:00 2001
From: martinjrobins <martinjrobins@gmail.com>
Date: Sat, 6 Jun 2026 12:23:18 +0000
Subject: [PATCH 2/4] fix: typo in dpow_f64

---
 diffsl/src/execution/functions.rs | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/diffsl/src/execution/functions.rs b/diffsl/src/execution/functions.rs
index 1cfc8f7..34f5af2 100644
--- a/diffsl/src/execution/functions.rs
+++ b/diffsl/src/execution/functions.rs
@@ -315,7 +315,7 @@ extern "C" fn pow_f64(x: f64, y: f64) -> f64 {
     x.powf(y)
 }
 extern "C" fn dpow_f64(x: f64, dx: f64, y: f64, dy: f64) -> f64 {
-    x.powf(y - 1.0) * (y * dx + x * dx.ln() * dy)
+    x.powf(y - 1.0) * (y * dx + x * x.ln() * dy)
 }
 
 extern "C" fn min_f64(x: f64, y: f64) -> f64 {
@@ -465,7 +465,7 @@ extern "C" fn pow_f32(x: f32, y: f32) -> f32 {
     x.powf(y)
 }
 extern "C" fn dpow_f32(x: f32, dx: f32, y: f32, dy: f32) -> f32 {
-    x.powf(y - 1.0) * (y * dx + x * dx.ln() * dy)
+    x.powf(y - 1.0) * (y * dx + x * x.ln() * dy)
 }
 
 extern "C" fn min_f32(x: f32, y: f32) -> f32 {

From e5ae0b067e6e377c6ec0b83e7c772ff636b7d1c1 Mon Sep 17 00:00:00 2001
From: martinjrobins <martinjrobins@gmail.com>
Date: Sat, 6 Jun 2026 13:30:12 +0000
Subject: [PATCH 3/4] fix: x==0 guard for pow

---
 diffsl/src/execution/functions.rs | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/diffsl/src/execution/functions.rs b/diffsl/src/execution/functions.rs
index 34f5af2..c872c2a 100644
--- a/diffsl/src/execution/functions.rs
+++ b/diffsl/src/execution/functions.rs
@@ -315,7 +315,12 @@ extern "C" fn pow_f64(x: f64, y: f64) -> f64 {
     x.powf(y)
 }
 extern "C" fn dpow_f64(x: f64, dx: f64, y: f64, dy: f64) -> f64 {
-    x.powf(y - 1.0) * (y * dx + x * x.ln() * dy)
+    let d_dy = if x == 0.0 {
+        0.0 // lim x→0+ x^y * ln(x) = 0 for y > 0, avoids NaN from 0 * (-inf)
+    } else {
+        x * x.ln() * dy
+    };
+    x.powf(y - 1.0) * (y * dx + d_dy)
 }
 
 extern "C" fn min_f64(x: f64, y: f64) -> f64 {
@@ -465,7 +470,12 @@ extern "C" fn pow_f32(x: f32, y: f32) -> f32 {
     x.powf(y)
 }
 extern "C" fn dpow_f32(x: f32, dx: f32, y: f32, dy: f32) -> f32 {
-    x.powf(y - 1.0) * (y * dx + x * x.ln() * dy)
+    let d_dy = if x == 0.0 {
+        0.0 // lim x→0+ x^y * ln(x) = 0 for y > 0, avoids NaN from 0 * (-inf)
+    } else {
+        x * x.ln() * dy
+    };
+    x.powf(y - 1.0) * (y * dx + d_dy)
 }
 
 extern "C" fn min_f32(x: f32, y: f32) -> f32 {

From 8d98ba5603f5bf1c52227b62ce8e9997f9fb8efc Mon Sep 17 00:00:00 2001
From: martinjrobins <martinjrobins@gmail.com>
Date: Mon, 8 Jun 2026 09:38:46 +0000
Subject: [PATCH 4/4] fix: pow(0, y) gradient bug in Enzyme, add test

---
 diffsl/Enzyme                    | 2 +-
 diffsl/src/execution/compiler.rs | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/diffsl/Enzyme b/diffsl/Enzyme
index db74433..6f0b4c9 160000
--- a/diffsl/Enzyme
+++ b/diffsl/Enzyme
@@ -1 +1 @@
-Subproject commit db7443336d020d3184cdc7eab00288fc34983e78
+Subproject commit 6f0b4c9945ada9461d054d7afd6455f78dacaa9f
diff --git a/diffsl/src/execution/compiler.rs b/diffsl/src/execution/compiler.rs
index 9ffe743..a9ed093 100644
--- a/diffsl/src/execution/compiler.rs
+++ b/diffsl/src/execution/compiler.rs
@@ -2875,7 +2875,7 @@ mod tests {
         state_and_const_grad2: "r_i { 3 * y, 2 * y }" expect "r" vec![3., 2.] ; vec![5.] ;  vec![0., 0.] ; vec![0.],
         state_and_const_grad3: "r_i { 2 * p, 3 }" expect "r" vec![2., 0.] ; vec![2.] ;  vec![2., 0.] ; vec![2.],
         state_and_const_grad4: "r_i { 3 * p, 2 * p }" expect "r" vec![3., 2.] ; vec![5.] ;  vec![3., 2.] ; vec![5.],
-
+        pow_zero_grad: "r { pow(y - p, p) }" expect "r" vec![0.] ; vec![0.] ; vec![-1.] ; vec![-1.],
     }
 
     macro_rules! tensor_test_big_state {