init2.c:51: assertion failed: p >= 2 && p <= ((mpfr_prec_t)((mpfr_prec_t)(~(mpfr_prec_t)0)>>1))

[Steve Kargl]

On Sun, Oct 01, 2006 at 03:13:51PM +0200, Dominique Dhumieres wrote:
> With gcc version 4.2.0 20060930 (experimental), gmp 4.2.1, and mpfr 2.2.0,
> on OSX 10.3.9/G5, the following code
> 
> real x
> x = -0.375*tiny(x)
> end
> 
> fails with:
> 
> init2.c:51:  assertion failed: p >= 2 && p <= ((mpfr_prec_t)((mpfr_prec_t)(~(mpfr_prec_t)0)>>1))
> test_spacing_red.f90:0: internal compiler error: Abort trap
> Please submit a full bug report,
> with preprocessed source if appropriate.
> See <URL:http://gcc.gnu.org/bugs.html> for instructions.
> 

Program received signal SIGABRT, Aborted.
0x287f45cf in kill () from /lib/libc.so.7
(gdb) bt
#0  0x287f45cf in kill () from /lib/libc.so.7
#1  0x287f456c in raise () from /lib/libc.so.7
#2  0x287f3280 in abort () from /lib/libc.so.7
#3  0x285f4961 in mpfr_assert_fail () from /usr/local/lib/libmpfr.so.1
#4  0x285eede4 in mpfr_init2 () from /usr/local/lib/libmpfr.so.1
#5  0x285f8e41 in mpfr_subnormalize () from /usr/local/lib/libmpfr.so.1
#6  0x0804c532 in gfc_check_real_range (p=0x88403b0, kind=Variable "kind" is not available.
)
    at ../../gcc4x/gcc/fortran/arith.c:449
#7  0x0804c631 in gfc_range_check (e=0x8840380)
    at ../../gcc4x/gcc/fortran/arith.c:607
#8  0x0804d6f6 in gfc_arith_times (op1=0x88402e0, op2=0x8840330, 
    resultp=0xbfbfe398) at ../../gcc4x/gcc/fortran/arith.c:832
#9  0x0804be4c in eval_intrinsic (operator=INTRINSIC_TIMES, eval=
      {f2 = 0x804d690 <gfc_arith_times>, f3 = 0x804d690 <gfc_arith_times>}, 
    op1=0x88402e0, op2=0x8840330) at ../../gcc4x/gcc/fortran/arith.c:1479
#10 0x08060865 in gfc_simplify_expr (p=0x88403d0, type=0)
    at ../../gcc4x/gcc/fortran/expr.c:802
#11 0x00000000 in ?? ()

It in the use of the new mpfr_subnormalize() function.  I'll look
into the problem.
-- 
Steve