[seed] [mpfr] Added bessel functions

[Tim Horton <hortont src gnome org>]

commit 71e62c1cc04edff826e8dda6a146a5a5454af889
Author: Matt ARSENAULT <arsenm2 rpi edu>
Date:   Sun Jul 5 16:12:44 2009 -0400

    [mpfr] Added bessel functions

 modules/mpfr/seed-mpfr-trig.c |  199 +++++++++++++++++++++++++++++++++++++++++
 modules/mpfr/seed-mpfr.c      |    6 ++
 modules/mpfr/seed-mpfr.h      |    7 ++
 3 files changed, 212 insertions(+), 0 deletions(-)
---
diff --git a/modules/mpfr/seed-mpfr-trig.c b/modules/mpfr/seed-mpfr-trig.c
index f60ff03..fc05db9 100644
--- a/modules/mpfr/seed-mpfr-trig.c
+++ b/modules/mpfr/seed-mpfr-trig.c
@@ -882,3 +882,202 @@ SeedValue seed_mpfr_erfc (SeedContext ctx,
     return seed_value_from_int(ctx, ret, exception);
 }
 
+/* bessel functions */
+SeedValue seed_mpfr_j0 (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.j0", 2);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[1], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) )
+    {
+        op = seed_object_get_private(args[0]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.j0", "mpfr_t");
+    }
+
+    ret = mpfr_j0(rop, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_j1 (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.j1", 2);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[1], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) )
+    {
+        op = seed_object_get_private(args[0]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.j1", "mpfr_t");
+    }
+
+    ret = mpfr_j1(rop, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_jn (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret, n;
+
+    CHECK_ARG_COUNT("mpfr.jn", 3);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[2], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[1], mpfr_class) )
+    {
+        op = seed_object_get_private(args[1]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.jn", "mpfr_t");
+    }
+
+    if ( seed_value_is_number(ctx, args[0]) )
+    {
+        n = seed_value_to_int(ctx, args[0], exception);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.jn", "int");
+    }
+
+    ret = mpfr_jn(rop, n, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_y0 (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.y0", 2);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[1], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) )
+    {
+        op = seed_object_get_private(args[0]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.y0", "mpfr_t");
+    }
+
+    ret = mpfr_y0(rop, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_y1 (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.y1", 2);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[1], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) )
+    {
+        op = seed_object_get_private(args[0]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.y1", "mpfr_t");
+    }
+
+    ret = mpfr_y1(rop, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_yn (SeedContext ctx,
+                        SeedObject function,
+                        SeedObject this_object,
+                        gsize argument_count,
+                        const SeedValue args[],
+                        SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op;
+    gint ret, n;
+
+    CHECK_ARG_COUNT("mpfr.yn", 3);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[2], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[1], mpfr_class) )
+    {
+        op = seed_object_get_private(args[1]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.yn", "mpfr_t");
+    }
+
+    if ( seed_value_is_number(ctx, args[0]) )
+    {
+        n = seed_value_to_int(ctx, args[0], exception);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.yn", "int");
+    }
+
+    ret = mpfr_yn(rop, n, op, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
diff --git a/modules/mpfr/seed-mpfr.c b/modules/mpfr/seed-mpfr.c
index 66ece2e..ea6c9c5 100644
--- a/modules/mpfr/seed-mpfr.c
+++ b/modules/mpfr/seed-mpfr.c
@@ -634,6 +634,12 @@ seed_static_function mpfr_funcs[] =
     {"zeta", seed_mpfr_zeta, 0},
     {"erf", seed_mpfr_erf, 0},
     {"erfc", seed_mpfr_erfc, 0},
+    {"j0", seed_mpfr_j0, 0},
+    {"j1", seed_mpfr_j1, 0},
+    {"jn", seed_mpfr_jn, 0},
+    {"y0", seed_mpfr_y0, 0},
+    {"y1", seed_mpfr_y1, 0},
+    {"yn", seed_mpfr_yn, 0},
     {"set", seed_mpfr_set, 0},
     {"set_inf", seed_mpfr_set_inf, 0},
     {"set_nan", seed_mpfr_set_nan, 0},
diff --git a/modules/mpfr/seed-mpfr.h b/modules/mpfr/seed-mpfr.h
index c7f7ced..0d8f535 100644
--- a/modules/mpfr/seed-mpfr.h
+++ b/modules/mpfr/seed-mpfr.h
@@ -113,6 +113,13 @@ DEF_SEED_MPFR_FUNC(seed_mpfr_zeta);
 DEF_SEED_MPFR_FUNC(seed_mpfr_erf);
 DEF_SEED_MPFR_FUNC(seed_mpfr_erfc);
 
+DEF_SEED_MPFR_FUNC(seed_mpfr_j0);
+DEF_SEED_MPFR_FUNC(seed_mpfr_j1);
+DEF_SEED_MPFR_FUNC(seed_mpfr_jn);
+DEF_SEED_MPFR_FUNC(seed_mpfr_y0);
+DEF_SEED_MPFR_FUNC(seed_mpfr_y1);
+DEF_SEED_MPFR_FUNC(seed_mpfr_yn);
+
 
 #endif      /* _SEED_MFPR_H_ */