Make pow with negative integer exponent be differentiable.

aspgomes · Oct 26, 2010 · 6e70bd3 · 6e70bd3
1 parent 985fe24
commit 6e70bd3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/taylor/taylor_math.h b/src/taylor/taylor_math.h
@@ -294,7 +294,7 @@ static taylor<T,Nvar,Ndeg> pow(const taylor<T,Nvar,Ndeg> &t, int n)
     }
   else if (n < 0)
     {
-      return pow(t,double(n));
+      return 1/pow(t,-n);
     }
   else 
     {