summaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
Diffstat (limited to 'src')
-rw-r--r--src/elementaryFunctions/pow/cpows.c11
-rw-r--r--src/elementaryFunctions/pow/zpows.c9
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));
}