[mpfr] Added mpfr.fma and mpfr.fms

diff --git a/modules/mpfr/seed-mpfr-trig.c b/modules/mpfr/seed-mpfr-trig.c
index fc05db9..765af38 100644 (file)
--- a/modules/mpfr/seed-mpfr-trig.c
+++ b/modules/mpfr/seed-mpfr-trig.c
@@ -1081,3 +1081,71 @@ SeedValue seed_mpfr_yn (SeedContext ctx,
     return seed_value_from_int(ctx, ret, exception);
 }
 
+SeedValue seed_mpfr_fma (SeedContext ctx,
+                         SeedObject function,
+                         SeedObject this_object,
+                         gsize argument_count,
+                         const SeedValue args[],
+                         SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op1, op2, op3;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.fma", 4);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[3], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) &&
+         seed_value_is_object_of_class(ctx, args[1], mpfr_class) &&
+         seed_value_is_object_of_class(ctx, args[2], mpfr_class))
+    {
+        op1 = seed_object_get_private(args[0]);
+        op2 = seed_object_get_private(args[1]);
+        op3 = seed_object_get_private(args[2]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.fma", "mpfr_t");
+    }
+
+    ret = mpfr_fma(rop, op1, op2, op3, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+
+SeedValue seed_mpfr_fms (SeedContext ctx,
+                         SeedObject function,
+                         SeedObject this_object,
+                         gsize argument_count,
+                         const SeedValue args[],
+                         SeedException * exception)
+{
+    mpfr_rnd_t rnd;
+    mpfr_ptr rop, op1, op2, op3;
+    gint ret;
+
+    CHECK_ARG_COUNT("mpfr.fms", 4);
+
+    rop = seed_object_get_private(this_object);
+    rnd = seed_value_to_mpfr_rnd_t(ctx, args[3], exception);
+
+    if ( seed_value_is_object_of_class(ctx, args[0], mpfr_class) &&
+         seed_value_is_object_of_class(ctx, args[1], mpfr_class) &&
+         seed_value_is_object_of_class(ctx, args[2], mpfr_class))
+    {
+        op1 = seed_object_get_private(args[0]);
+        op2 = seed_object_get_private(args[1]);
+        op3 = seed_object_get_private(args[2]);
+    }
+    else
+    {
+        TYPE_EXCEPTION("mpfr.fms", "mpfr_t");
+    }
+
+    ret = mpfr_fms(rop, op1, op2, op3, rnd);
+
+    return seed_value_from_int(ctx, ret, exception);
+}
+

diff --git a/modules/mpfr/seed-mpfr.c b/modules/mpfr/seed-mpfr.c
index ea6c9c5..28fad31 100644 (file)
--- a/modules/mpfr/seed-mpfr.c
+++ b/modules/mpfr/seed-mpfr.c
@@ -640,6 +640,8 @@ seed_static_function mpfr_funcs[] =
     {"y0", seed_mpfr_y0, 0},
     {"y1", seed_mpfr_y1, 0},
     {"yn", seed_mpfr_yn, 0},
+    {"fma", seed_mpfr_fma, 0},
+    {"fms", seed_mpfr_fms, 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 0d8f535..8cda045 100644 (file)
--- a/modules/mpfr/seed-mpfr.h
+++ b/modules/mpfr/seed-mpfr.h
@@ -120,6 +120,9 @@ DEF_SEED_MPFR_FUNC(seed_mpfr_y0);
 DEF_SEED_MPFR_FUNC(seed_mpfr_y1);
 DEF_SEED_MPFR_FUNC(seed_mpfr_yn);
 
+DEF_SEED_MPFR_FUNC(seed_mpfr_fma);
+DEF_SEED_MPFR_FUNC(seed_mpfr_fms);
+
 
 #endif      /* _SEED_MFPR_H_ */