On Sat, Jan 30, 2010 at 2:54 PM, Bill Hart <goodwillh...@googlemail.com> wrote:
> I've placed an MPIR 1.3.1 which hopefully fixes this issue, here:
>
> http://www.mpir.org/mpir-1.3.1.tar.gz
>
> Let me know if there are further problems.
I confirmed that the problem existed in gmpy with 1.3.0 and is fixed in 1.3.1.
casevh
>
> Bill.
>
> 2010/1/30 Bill Hart <goodwillh...@googlemail.com>:
>> Hi Dan,
>>
>> Yes, there is a bug. The fix is to replace line 1130 of mpn/gcdext.c with:
>>
>> *s0size = 1 - (s0p[0] == 0);
>>
>> I will issue a service release shortly to fix this issue.
>>
>> Can you let us know if this fixes your (other) program.
>>
>> Bill.
>>
>> 2010/1/30 Bill Hart <goodwillh...@googlemail.com>:
>>> Hi Dan,
>>>
>>> That definitely looks like a bug. You are correct in that size should
>>> be set to zero.
>>>
>>> It's slightly bizarre in that we haven't altered mpz_gcdext, that I
>>> recall. But we did modify mpn_gcdext, so I assume the bug is there.
>>>
>>> I'll take a look.
>>>
>>> Bill.
>>>
>>> 2010/1/30 Dan Grayson <d...@math.uiuc.edu>:
>>>> #include <mpir.h>
>>>> #include <assert.h>
>>>>
>>>> int main () {
>>>> mpz_t a,b,d,u,v;
>>>> mpz_init(a);
>>>> mpz_init(b);
>>>> mpz_init(d);
>>>> mpz_init(u);
>>>> mpz_init(v);
>>>> mpz_set_si(a,1);
>>>> mpz_set_si(b,2);
>>>> mpz_gcdext(d,u,v,a,b);
>>>> assert(!(u->_mp_size == 1 && u->_mp_d[0] == 0));
>>>> assert(!(v->_mp_size == 1 && v->_mp_d[0] == 0));
>>>> return 0;
>>>> }
>>>>
>>>> produces an assertion failure with 1.3.0 but not with 1.2.1. I
>>>> believe that could be a bug, because the integer 0 should have the
>>>> size set to 0, according to the documentation; in any case, our
>>>> (other) program doesn't work. The OS is Mac OS X 10.6.
>>>>
>>>> --
>>>> You received this message because you are subscribed to the Google Groups
>>>> "mpir-devel" group.
>>>> To post to this group, send email to mpir-de...@googlegroups.com.
>>>> To unsubscribe from this group, send email to
>>>> mpir-devel+unsubscr...@googlegroups.com.
>>>> For more options, visit this group at
>>>> http://groups.google.com/group/mpir-devel?hl=en.
>>>>
>>>>
>>>
>>
>
> --
> You received this message because you are subscribed to the Google Groups
> "mpir-devel" group.
> To post to this group, send email to mpir-de...@googlegroups.com.
> To unsubscribe from this group, send email to
> mpir-devel+unsubscr...@googlegroups.com.
> For more options, visit this group at
> http://groups.google.com/group/mpir-devel?hl=en.
>
>
--
You received this message because you are subscribed to the Google Groups
"mpir-devel" group.
To post to this group, send email to mpir-de...@googlegroups.com.
To unsubscribe from this group, send email to
mpir-devel+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/mpir-devel?hl=en.
