ext/mpfi/func_mpfi_extention.c in ruby-mpfi-0.0.9 vs ext/mpfi/func_mpfi_extention.c in ruby-mpfi-0.0.10

- old
+ new

@@ -2,11 +2,11 @@ #include <mpfi.h> #include <mpfi_io.h> #include "ruby_mpfi.h" #include "func_mpfi_extention.h" -void mpfi_mid_interval(mpfi_t ret, mpfi_t x){ +void mpfi_mid_interval (mpfi_t ret, mpfi_t x) { mpfr_t left, right; mpfr_init2(left, mpfi_get_prec(x)); mpfr_init2(right, mpfi_get_prec(x)); mpfr_add(left, r_mpfi_left_ptr(x), r_mpfi_right_ptr(x), GMP_RNDD); @@ -19,25 +19,24 @@ mpfr_clear(right); } /* Retrun 0 if this function puts subdivision to *ret. */ /* Otherwise, return -1. */ -int mpfi_subdivision(int num, mpfi_t *ret, mpfi_t x){ +int mpfi_subdivision (int num, mpfi_t *ret, mpfi_t x) { int i, ret_val = -1; - mpfr_t l; + mpfr_t l, x_diam; mpfr_init(l); mpfr_sub(l, r_mpfi_right_ptr(x), r_mpfi_left_ptr(x), GMP_RNDD); mpfr_div_si(l, l, num, GMP_RNDD); - mpfr_t x_diam; mpfr_init(x_diam); mpfi_diam_abs(x_diam, x); - if(mpfr_cmp(x_diam, l) > 0){ + if (mpfr_cmp(x_diam, l) > 0) { mpfr_set(r_mpfi_left_ptr(*ret), r_mpfi_left_ptr(x), GMP_RNDN); mpfr_add(r_mpfi_right_ptr(*ret), r_mpfi_left_ptr(*ret), l, GMP_RNDU); - for(i = 1; i < num - 1; i ++){ + for (i = 1; i < num - 1; i ++) { mpfr_set(r_mpfi_left_ptr(*ret + i), r_mpfi_right_ptr(*ret + i - 1), GMP_RNDN); mpfr_add(r_mpfi_right_ptr(*ret + i), r_mpfi_left_ptr(*ret + i), l, GMP_RNDU); } mpfr_set(r_mpfi_left_ptr(*ret + i), r_mpfi_right_ptr(*ret + i - 1), GMP_RNDN); @@ -47,6 +46,5 @@ mpfr_clear(x_diam); mpfr_clear(l); return ret_val; } -