Brilliant. Thanks for letting us know. On 1 February 2010 13:25, Dan Grayson <d...@math.uiuc.edu> wrote: > It seems that 1.3.1 fixes the problem. > > On Jan 30, 4: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. >> >> 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 >> >>> athttp://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.