diff --git a/acb_hypgeom/legendre_p_uiui_rec.c b/acb_hypgeom/legendre_p_uiui_rec.c
index cf9ba3d0..79fcf7ae 100644
--- a/acb_hypgeom/legendre_p_uiui_rec.c
+++ b/acb_hypgeom/legendre_p_uiui_rec.c
@@ -57,7 +57,7 @@ acb_hypgeom_legendre_p_uiui_rec(acb_t res, ulong n, ulong m, const acb_t z, slon
     if (m == 0)
         arb_one(acb_realref(t));
     else
-        arb_fac2_ui(acb_realref(t), 2 * m - 1, prec);
+        arb_doublefac_ui(acb_realref(t), 2 * m - 1, prec);
 
     if (m % 2)
         arb_neg(acb_realref(t), acb_realref(t));
diff --git a/arb.h b/arb.h
index fe3c07eb..a8f73d8d 100644
--- a/arb.h
+++ b/arb.h
@@ -617,7 +617,7 @@ void arb_atanh(arb_t z, const arb_t x, slong prec);
 void arb_asinh(arb_t z, const arb_t x, slong prec);
 void arb_acosh(arb_t z, const arb_t x, slong prec);
 void arb_fac_ui(arb_t z, ulong n, slong prec);
-void arb_fac2_ui(arb_t z, ulong n, slong prec);
+void arb_doublefac_ui(arb_t z, ulong n, slong prec);
 void arb_bin_ui(arb_t z, const arb_t n, ulong k, slong prec);
 void arb_bin_uiui(arb_t z, ulong n, ulong k, slong prec);
 void arb_fib_fmpz(arb_t z, const fmpz_t n, slong prec);
diff --git a/arb/fac2_ui.c b/arb/doublefac_ui.c
similarity index 93%
rename from arb/fac2_ui.c
rename to arb/doublefac_ui.c
index 1d6db717..f899bf51 100644
--- a/arb/fac2_ui.c
+++ b/arb/doublefac_ui.c
@@ -26,7 +26,7 @@
 #include "arb.h"
 
 void
-arb_fac2_ui(arb_t res, ulong n, slong prec)
+arb_doublefac_ui(arb_t res, ulong n, slong prec)
 {
     if (n % 2 == 0)
     {
@@ -37,7 +37,7 @@ arb_fac2_ui(arb_t res, ulong n, slong prec)
     {
         arb_t t;
         arb_init(t);
-        arb_fac2_ui(t, n - 1, prec + 5);
+        arb_doublefac_ui(t, n - 1, prec + 5);
         arb_fac_ui(res, n, prec + 5);
         arb_div(res, res, t, prec);
         arb_clear(t);
diff --git a/arb/test/t-fac2_ui.c b/arb/test/t-doublefac_ui.c
similarity index 94%
rename from arb/test/t-fac2_ui.c
rename to arb/test/t-doublefac_ui.c
index 1f0a9f86..cf1ea4fc 100644
--- a/arb/test/t-fac2_ui.c
+++ b/arb/test/t-doublefac_ui.c
@@ -30,7 +30,7 @@ int main()
     slong iter;
     flint_rand_t state;
 
-    flint_printf("fac2_ui....");
+    flint_printf("doublefac_ui....");
     fflush(stdout);
 
     flint_randinit(state);
@@ -52,8 +52,8 @@ int main()
         if (n + 1 == 0 || n + 2 == 0)
             n -= 2;
 
-        arb_fac2_ui(a, n, prec1);
-        arb_fac2_ui(b, n + 2, prec1);
+        arb_doublefac_ui(a, n, prec1);
+        arb_doublefac_ui(b, n + 2, prec1);
         arb_mul_ui(c, a, n + 2, prec2);
 
         if (!arb_overlaps(b, c))
diff --git a/doc/source/arb.rst b/doc/source/arb.rst
index 691b3f17..fe126bc0 100644
--- a/doc/source/arb.rst
+++ b/doc/source/arb.rst
@@ -1031,7 +1031,7 @@ Gamma function and factorials
 
     Computes the factorial `z = n!` via the gamma function.
 
-.. function:: void arb_fac2_ui(arb_t z, ulong n, slong prec)
+.. function:: void arb_doublefac_ui(arb_t z, ulong n, slong prec)
 
     Computes the double factorial `z = n!!` via the gamma function.