rename acb_holonomic to acb_ode

Author: mezzarobba

Makefile.in

@@ -203,7 +203,7 @@ HEADER_DIRS :=                                                              \
         acb_mat         acb_poly        acb_calc        acb_hypgeom         \
         arb_fmpz_poly   arb_fpwrap                                          \
         acb_dft         acb_elliptic    acb_modular     acb_dirichlet       \
-        acb_theta       acb_holonomic                                       \
+        acb_theta       acb_ode                                             \
         dirichlet       bernoulli       hypgeom                             \
                                                                             \
         gr              gr_generic      gr_vec          gr_mat              \

dev/check_examples.sh

@@ -169,9 +169,9 @@ then
 elif test "$1" = "hilbert_matrix_ca";
 then
     echo "hilbert_matrix_ca....SKIPPED"
-elif test "$1" = "holonomic";
+elif test "$1" = "ode";
 then
-    echo "holonomic....SKIPPED"
+    echo "ode....SKIPPED"
 elif test "$1" = "huge_expr";
 then
     echo "huge_expr....SKIPPED"

examples/holonomic.c examples/ode.cexamples/holonomic.c renamed to examples/ode.c

@@ -1,15 +1,15 @@
 #include "acb_types.h"
 #include "acb_poly.h"
 #include "acb_mat.h"
-#include "acb_holonomic.h"
+#include "acb_ode.h"
 #include "gr.h"
 #include "gr_poly.h"
 
 
 void
-_acb_holonomic_sum_swap_mat_ordinary(
+_acb_ode_sum_swap_mat_ordinary(
         acb_mat_t mat,
-        const acb_holonomic_sum_context_struct * ctx, slong s)
+        const acb_ode_sum_context_struct * ctx, slong s)
 {
     for (slong j = 0; j < ctx->nsols; j++)
     {
@@ -30,11 +30,11 @@ _acb_holonomic_sum_swap_mat_ordinary(
 void
 ordinary(void)
 {
-    acb_holonomic_sum_context_t ctx;
+    acb_ode_sum_context_t ctx;
 
     /* (x^2 + 1)*Dx^3 + 7*x
      * = (x^2 + 1)*Tx^3 + (-3*x^2 - 3)*Tx^2 + (2*x^2 + 2)*Tx + 7*x^4 */
-    acb_holonomic_sum_context_init(ctx, 4, 2, 3, 3);
+    acb_ode_sum_context_init(ctx, 4, 2, 3, 3);
     acb_poly_set_coeff_si(ctx->dop + 3, 2, 1);
     acb_poly_set_coeff_si(ctx->dop + 3, 0, 1);
     acb_poly_set_coeff_si(ctx->dop + 2, 2, -3);
@@ -43,10 +43,10 @@ ordinary(void)
     acb_poly_set_coeff_si(ctx->dop + 1, 0, 2);
     acb_poly_set_coeff_si(ctx->dop + 0, 4, 7);
 
-    /* ctx->flags |= ACB_HOLONOMIC_WANT_SERIES; */
+    /* ctx->flags |= acb_ode_WANT_SERIES; */
 
-    acb_holonomic_sum_ordinary(ctx);
-    acb_holonomic_sum_canonical_basis(ctx);
+    acb_ode_sum_ordinary(ctx);
+    acb_ode_sum_canonical_basis(ctx);
 
     acb_set_d_d(ctx->pts, 0.25, 0.25);
     acb_set_si(ctx->pts + 1, 0);
@@ -55,86 +55,86 @@ ordinary(void)
     ctx->prec = 64;
     ctx->sums_prec = 64;
 
-    acb_holonomic_sum_divconquer(ctx, 20);
+    acb_ode_sum_divconquer(ctx, 20);
 
     acb_mat_t mat;
     acb_mat_init(mat, ctx->nder, ctx->nsols);
 
     for (slong i = 0; i < ctx->npts; i++)
     {
-        _acb_holonomic_sum_swap_mat_ordinary(mat, ctx, i);
+        _acb_ode_sum_swap_mat_ordinary(mat, ctx, i);
         acb_mat_printd(mat, 8);
         /* flint_printf("%{acb_mat}\n\n", mat); */
     }
 
-    acb_holonomic_sum_context_clear(ctx);
+    acb_ode_sum_context_clear(ctx);
     acb_mat_clear(mat);
 }
 
 
 void
 series(void)
 {
-    acb_holonomic_sum_context_t ctx;
+    acb_ode_sum_context_t ctx;
 
-    acb_holonomic_sum_context_init(ctx, 2, 0, 1, 1);
+    acb_ode_sum_context_init(ctx, 2, 0, 1, 1);
     acb_poly_set_coeff_si(ctx->dop + 1, 0, 1);
     acb_poly_set_coeff_si(ctx->dop + 0, 1, -1);
 
-    acb_holonomic_sum_ordinary(ctx);
-    acb_holonomic_sum_canonical_basis(ctx);
+    acb_ode_sum_ordinary(ctx);
+    acb_ode_sum_canonical_basis(ctx);
     ctx->prec = 64;
-    ctx->flags |= ACB_HOLONOMIC_WANT_SERIES;
+    ctx->flags |= acb_ode_WANT_SERIES;
 
     slong len = 5;
-    acb_holonomic_sum_divconquer(ctx, len);
+    acb_ode_sum_divconquer(ctx, len);
     /* Clear the high part reserved for the residual. (In this special case, the
      * high part is zero because block_length = 1.) */
     acb_poly_truncate(ctx->sol[0].series, len);
 
     flint_printf("%{acb_poly}\n", ctx->sol[0].series);
 
-    acb_holonomic_sum_context_clear(ctx);
+    acb_ode_sum_context_clear(ctx);
 }
 
 
 void
 bessel_j0(void)
 {
-    acb_holonomic_sum_context_t ctx;
+    acb_ode_sum_context_t ctx;
 
     slong dop_order = 2;
-    acb_holonomic_sum_context_init(ctx, dop_order + 1, 1, dop_order, dop_order);
+    acb_ode_sum_context_init(ctx, dop_order + 1, 1, dop_order, dop_order);
     acb_poly_set_coeff_si(ctx->dop + 2, 0, 1);
     acb_poly_set_coeff_si(ctx->dop + 0, 2, 1);
 
     ctx->sing_shifts[0].n = 0;
     ctx->sing_shifts[0].mult = 2;
 
-    acb_holonomic_sum_canonical_basis(ctx);
+    acb_ode_sum_canonical_basis(ctx);
 
     acb_set_d(ctx->pts, 0.25);
 
     ctx->prec = 64;
     ctx->sums_prec = 64;
 
-    acb_holonomic_sum_divconquer(ctx, 10);
+    acb_ode_sum_divconquer(ctx, 10);
 
     for (slong m = 0; m < ctx->nsols; m++)
     {
-        acb_holonomic_sol_struct * sol = ctx->sol + m;
+        acb_ode_sol_struct * sol = ctx->sol + m;
         /* series in x */
         flint_printf("f%wd =", m);
         for (slong k = 0; k < sol->nlogs; k++)
         {
             flint_printf(" + (%{acb_poly})*log(%{acb} + x)^%wd/%wd",
-                         acb_holonomic_sol_sum_ptr(sol, 0, k),
+                         acb_ode_sol_sum_ptr(sol, 0, k),
                          ctx->pts, k, k);
         }
         flint_printf("\n");
     }
 
-    acb_holonomic_sum_context_clear(ctx);
+    acb_ode_sum_context_clear(ctx);
 }
 
 
@@ -149,7 +149,7 @@ whittaker_m(void)  /* non-integer exponent, no logs */
     acb_init(half);
     acb_poly_init(val);
 
-    acb_holonomic_sum_context_t ctx;
+    acb_ode_sum_context_t ctx;
 
     acb_set_si(kappa, 2);
     acb_set_si(mu2, 3);
@@ -159,7 +159,7 @@ whittaker_m(void)  /* non-integer exponent, no logs */
     slong dop_order = 2;
     slong len = dop_order;
 
-    acb_holonomic_sum_context_init(ctx, dop_order + 1, 1, 1, len);
+    acb_ode_sum_context_init(ctx, dop_order + 1, 1, 1, len);
 
     acb_poly_set_coeff_si(ctx->dop + 2, 0, 4);
     acb_poly_set_coeff_si(ctx->dop + 1, 0, -4);
@@ -174,26 +174,26 @@ whittaker_m(void)  /* non-integer exponent, no logs */
     ctx->sing_shifts[0].n = 0;
     ctx->sing_shifts[0].mult = 1;
 
-    acb_holonomic_sum_canonical_basis(ctx);
+    acb_ode_sum_canonical_basis(ctx);
 
     acb_set_d(ctx->pts, 1.4242);
 
     ctx->prec = prec;
     ctx->sums_prec = prec;
 
-    acb_holonomic_sum_divconquer(ctx, prec);
+    acb_ode_sum_divconquer(ctx, prec);
 
-    acb_poly_struct * f = acb_holonomic_sol_sum_ptr(ctx->sol, 0, 0);
+    acb_poly_struct * f = acb_ode_sol_sum_ptr(ctx->sol, 0, 0);
 
     flint_printf("(%{acb} + x)^(%{acb}) * (%{acb_poly})\n",
                  ctx->pts, ctx->expo, f);
 
-    _acb_holonomic_sol_value(val, ctx->expo, f, ctx->sol[0].nlogs, ctx->pts,
-                             ctx->nder, 1, ctx->prec);
+    _acb_ode_sol_value(val, ctx->expo, f, ctx->sol[0].nlogs, ctx->pts,
+                       ctx->nder, 1, ctx->prec);
 
     flint_printf("M(%{acb} + x) = %{acb_poly} + O(x^%wd)\n", ctx->pts, val, len);
 
-    acb_holonomic_sum_context_clear(ctx);
+    acb_ode_sum_context_clear(ctx);
     acb_clear(kappa);
     acb_clear(mu2);
     acb_clear(half);
@@ -203,7 +203,7 @@ whittaker_m(void)  /* non-integer exponent, no logs */
 
 void
 fundamental_matrix(const char * dop_str,
-                   const acb_holonomic_exponents_struct * expos,
+                   const acb_ode_exponents_struct * expos,
                    double pt_d)
 {
     gr_ctx_t CC, Pol, Dop;
@@ -234,7 +234,7 @@ fundamental_matrix(const char * dop_str,
     acb_init(pt);
     acb_set_d(pt, pt_d);
 
-    acb_holonomic_fundamental_matrix(mat, dop, Dop, expos, pt, 1, 0, 8, prec);
+    acb_ode_fundamental_matrix(mat, dop, Dop, expos, pt, 1, 0, 8, prec);
 
     flint_printf("%{acb_mat}\n", mat);
 
@@ -250,11 +250,11 @@ fundamental_matrix(const char * dop_str,
 void
 apery(void)
 {
-    acb_holonomic_shift_struct shift[1] = {{ .n = 0, .mult = 3 }};
-    acb_holonomic_group_struct grp[1] = {{ .nshifts = 1, .shifts = shift }};
+    acb_ode_shift_struct shift[1] = {{ .n = 0, .mult = 3 }};
+    acb_ode_group_struct grp[1] = {{ .nshifts = 1, .shifts = shift }};
     acb_init(grp->expo);
     acb_zero(grp->expo);
-    acb_holonomic_exponents_struct expos[1] = {{ .len = 1, .grps = grp }};
+    acb_ode_exponents_struct expos[1] = {{ .len = 1, .grps = grp }};
 
     fundamental_matrix(
             "(z^2 - 34*z + 1)*Tz^3 + (3*z^2 - 51*z)*Tz^2 + (3*z^2 - 27*z)*Tz + z^2 - 5*z",
@@ -271,15 +271,15 @@ multiple_shifts(void)
     /* const char * dop = "Tz^4 - 4*Tz^3 + 3*Tz^2 - z"; */
     const char * dop = "Tz^6 - 6*Tz^5 + 12*Tz^4 - 10*Tz^3 + 3*Tz^2 + z^2";
 
-    acb_holonomic_shift_struct shift[3] = {
+    acb_ode_shift_struct shift[3] = {
         { .n = 0, .mult = 2 },
         { .n = 1, .mult = 3 },
         { .n = 3, .mult = 1 },
     };
-    acb_holonomic_group_struct grp[1] = {{ .nshifts = 3, .shifts = shift }};
+    acb_ode_group_struct grp[1] = {{ .nshifts = 3, .shifts = shift }};
     acb_init(grp->expo);
     acb_zero(grp->expo);
-    acb_holonomic_exponents_struct expos[1] = {{ .len = 1, .grps = grp }};
+    acb_ode_exponents_struct expos[1] = {{ .len = 1, .grps = grp }};
 
     fundamental_matrix(dop, expos, 2.);