add missing fpwrap pow

Author: fredrik-johansson

arb_fpwrap.h

@@ -51,6 +51,9 @@ int arb_fpwrap_cdouble_log(complex_double * res, complex_double x, int flags);
 int arb_fpwrap_double_log1p(double * res, double x, int flags);
 int arb_fpwrap_cdouble_log1p(complex_double * res, complex_double x, int flags);
 
+int arb_fpwrap_double_pow(double * res, double x, double y, int flags);
+int arb_fpwrap_cdouble_pow(complex_double * res, complex_double x, complex_double y, int flags);
+
 int arb_fpwrap_double_sqrt(double * res, double x, int flags);
 int arb_fpwrap_cdouble_sqrt(complex_double * res, complex_double x, int flags);
 

arb_fpwrap/fpwrap.c

@@ -869,6 +869,9 @@ DEF_CDOUBLE_FUN_1(log, acb_log)
 DEF_DOUBLE_FUN_1(log1p, arb_log1p)
 DEF_CDOUBLE_FUN_1(log1p, acb_log1p)
 
+DEF_DOUBLE_FUN_2(pow, arb_pow)
+DEF_CDOUBLE_FUN_2(pow, acb_pow)
+
 DEF_DOUBLE_FUN_1(sqrt, arb_sqrt)
 DEF_CDOUBLE_FUN_1(sqrt, acb_sqrt)
 

arb_fpwrap/test/t-fpwrap.c

@@ -182,6 +182,9 @@ int main()
         CHECK_DOUBLE(arb_fpwrap_double_log1p(&res, x, flags));
         CHECK_CDOUBLE(arb_fpwrap_cdouble_log1p(&cres, cx, flags));
 
+        CHECK_DOUBLE(arb_fpwrap_double_pow(&res, x, y, flags));
+        CHECK_CDOUBLE(arb_fpwrap_cdouble_pow(&cres, cx, cy, flags));
+
         CHECK_DOUBLE(arb_fpwrap_double_sqrt(&res, x, flags));
         CHECK_CDOUBLE(arb_fpwrap_cdouble_sqrt(&cres, cx, flags));
 

doc/source/arb_fpwrap.rst

@@ -133,6 +133,9 @@ Elementary functions
 .. function:: int arb_fpwrap_double_log1p(double * res, double x, int flags)
               int arb_fpwrap_cdouble_log1p(complex_double * res, complex_double x, int flags)
 
+.. function:: int arb_fpwrap_double_pow(double * res, double x, double y, int flags)
+              int arb_fpwrap_cdouble_pow(complex_double * res, complex_double x, complex_double y, int flags)
+
 .. function:: int arb_fpwrap_double_sqrt(double * res, double x, int flags)
               int arb_fpwrap_cdouble_sqrt(complex_double * res, complex_double x, int flags)