ldbl_mant_type
ldbl_mant_type mantbits, sq_mantbits;
mantbits = ((ldbl_mant_type)1 << (mantdig / 2)) - 1;
mantbits > ((ldbl_mant_type)1 << (mantdig / 2 - 1));