Add complex pow exception : z=0

2 files changed, 19 insertions, 1 deletions

diff --git a/src/elementaryFunctions/pow/cpows.c b/src/elementaryFunctions/pow/cpows.c
index 87392706..5848c4b5 100644
--- a/src/elementaryFunctions/pow/cpows.c
+++ b/src/elementaryFunctions/pow/cpows.c
@@ -15,5 +15,14 @@
 #include "log.h"
 
 floatComplex		cpows(floatComplex z, floatComplex power) {
-    return cexps(cmuls(clogs(z), power));
+      /*Cas z=0 */
+  if ( (creals(z)==0) && (cimags(z)==0) ){
+  	/* Cas 0^0 */
+  	if ( (creals(power)==0) && (cimags(power)==0) ) return FloatComplex(1,0);
+  	/* Cas 0^x, x!=0 */
+  	return FloatComplex(0,0);
+  }
+  
+  /* Cas z!=0 */
+  return cexps(cmuls(clogs(z), power));
 }
diff --git a/src/elementaryFunctions/pow/zpows.c b/src/elementaryFunctions/pow/zpows.c
index cc532090..1a7059b2 100644
--- a/src/elementaryFunctions/pow/zpows.c
+++ b/src/elementaryFunctions/pow/zpows.c
@@ -15,5 +15,14 @@
 #include "exp.h"
 
 doubleComplex		zpows(doubleComplex z, doubleComplex power) {
+  /*Cas z=0 */
+  if ( (zreals(z)==0) && (zimags(z)==0) ){
+  	/* Cas 0^0 */
+  	if ( (zreals(power)==0) && (zimags(power)==0) ) return DoubleComplex(1,0);
+  	/* Cas 0^x, x!=0 */
+  	return DoubleComplex(0,0);
+  }
+  
+  /* Cas z!=0 */
   return zexps(zmuls(zlogs(z), power));
 }