|
Bill Allombert on Wed, 03 May 2023 00:07:52 +0200
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
|
Re: gcdext() for polynomials
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: Re: gcdext() for polynomials
- From: Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>
- Date: Wed, 3 May 2023 00:06:30 +0200
- Arc-authentication-results: i=1; smail; arc=none
- Arc-message-signature: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1683065185; c=relaxed/relaxed; bh=9fePTfBuVxYlzWAH2i/TXsD+a7lKT4CmCmxgi7xD4to=; h=DKIM-Signature:Date:From:To:Subject:Message-ID:Mail-Followup-To: References:MIME-Version:Content-Type:Content-Disposition: In-Reply-To; b=cwy30W376wCt/GjVEEEhPrF0nv1ZDc0TREIuR8iJlSjzlEQKxMBE95yXDVUv+ClIugAy1DXtAHn+PCjbbemkYFUX27wE4bXia+Ij3UNibNDTrJ/4rwDkkFnQQK4XElRLiIiDLPN5un1VKOP9/E6B+LjZKArNM7qJ7CIo4KrprR/ixzqHSX2ugvypHNiNAG4Ju66uJXIBAwaIxgbkQC1t4Sh6Yiozyzf6i/oLW92la3TfLItSogJzErdP5k29bhUBS41nA1lDvHZSdV6Hr3rIZtqBy7PMbiyzqOU7dQDBBZVjXNt/LsgMSSM6zivlRhKflfGV4mGd1/S6wDBc8mIiKphb4vjJnUkWF9/o2rAD9SmDLdMDCL7jAmvWQLd09r5zJEzfPDrFdp/VGzs1oizdQMDnzdLoRcv35J0RZ6hU5KPIq+ZjLXtwZ5OHlJ3TCiOCn62nmkTA0IBDX/g0Vv9pxZ3qj9xpU6iHJwi89rzvUp+AdSfmQ4V57CuM2CLsqVPINHn4KsoH90Wl72ATy+ubdlCvKAuxo5bbLr083SUgSNAw+E9XyaLkEX5KSijvhISlvzpTBlH1NEWm3V2RxepN2KeFWuBvyhgAwpq1Buca+8ZgVY9UhmznvBkto4UQCMlKNGrNb71nv2NFXHSdpVVrwoobscEiZcT0D10JwXtBedI=
- Arc-seal: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1683065185; cv=none; b=HwyuWhklgwYAHXP9kEyDMAD9zLWP/QvXDX2Db7HNIRR+RfcY8/obP8t+AgsGxTY+akYRofVosD4BYL3vhcLXgdmt0B7wwpnq8l+hOxIFPLQ88LI+8qq1+TehoI1swZgW3IfJ3kfwpxSb1QoQpMMDefse/gR4QqT+h9phMxpV95ZupKXPGU1gLoHL5ED8zZ9zyXUrhOiY0B9JKdELf/+LhLZGfqh0TdyIL2Wfv7UbKM373MWugMCzO/jUB2fnbqPuBwjCRXuiXYccIGlsrf2WOy+YFGHrY1OxBxH7JougTeoq/r7MA8yZwjzFD3J4F2BaLncLTnSwi1rDpZbWWBXzHLsOEMgK1RugRXW9HhKiKioYHlhbWxfkP5FPlcTtpNdOua2+JL8Ubr+KHUJuL/SlCo6WxhDvl1OQrApxzrVnBhZsj8OtrmrxuXY3NJZO96lN2ww7e3RSXOqlwvnu3b5JZV5152zmPvxDrP+yPb9FK2Srg4XyHrHsa0gfFFRTzpdiuj5qtLnvCB69vjgrtRb5Yjavs64tVtHssHfX8eBuwPkw6Kbjf41EEr0hGAs3zXc7k7i59x1hTtVYpDkWwtOg27Vr0Ah4uzkE9CTAVe/aaC25pzmQKNLf0LsnuNrTYuyi6BlA+ywl4Ug6tE4uFC+5OygzPvCAFGAiklbUIBdo6S4=
- Authentication-results: smail; arc=none
- Delivery-date: Wed, 03 May 2023 00:07:52 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=math.u-bordeaux.fr; s=2022; t=1683065185; bh=9fePTfBuVxYlzWAH2i/TXsD+a7lKT4CmCmxgi7xD4to=; h=Date:From:To:Subject:References:In-Reply-To:From; b=EMc7nQfgbnaF0bO1obtr6Xqk00BphzmG7xC2FrYdKeFE9aPlJUjf2RkImLHFOW+Hf F2h0HRzE5PR05ZpKg3QwDStmKVxqwHDobY4Ro/3HUz/KhB4T+3Wn5VRWamTFVTZ17a CnUE/uQ9oJTySBdU/56iMdY3OBv+226pWhACOPjgHIMMMnueA58XOCe46rsWYGT0XX EadfInJdL/WwtkaqXsGTzJqFWxNCNYS938T0Et6gijnjtk8wTnKQwAOHu06GA/RKkG fi3NwykQTz5PgC7yIpBu2jDR/fdjTI3clKbxdnkZj8htwlmabpSGYLLTytyyBft7B4 iY4tDGdlCvM2mOCWAw2JH8nhsJENo7UVRSrZDYEm1Q+o4aW0sbuB+L8pQxFYpjnqJG jipLFQqQbzBRfncho+/3fGxTey9mMPfnw0kRQE9YsGNuywxT6hnO/9Nkc7uvIkamUJ ysrGTRH/dRMeETayAniF3Pd2uMQBjos0KmuGZvvzyl7YylMpiA5stcQiSZ7FuigYCy SRH62nIb50tKc9zVVX1RSn7rkzMA7aAMc8A57j7hMAdKlpWJ6vdzz/CVIocsXkGXTW rxWzjishgQ7YiOw8jjVRIbWXA41ed+EiTja73kFl+oPSC6T3TEGQoZK96kwWzPBvRV wmlmW3+YlffQ+BlGRpWQD2xY=
- In-reply-to: <CAJkPp5P0ToGSDLYpHW0YCYFp5Zm-G+d7Q6YFkHJsaLYNq6vwZQ@mail.gmail.com>
- Mail-followup-to: pari-users@pari.math.u-bordeaux.fr
- References: <CAJkPp5P942=W+L32dtXi4o7UXN+1hH8dzuExjmyoET5N2Ccfog@mail.gmail.com> <ZE0uKigxKUxKmjyE@seventeen> <CAJkPp5P0ToGSDLYpHW0YCYFp5Zm-G+d7Q6YFkHJsaLYNq6vwZQ@mail.gmail.com>
On Sat, Apr 29, 2023 at 11:26:53AM -0400, Max Alekseyev wrote:
> I understand, but it should not be a trouble to add an optional parameter
> v, and handle it via renaming - like this:
>
{
mygcdext(f,g,v) = my(x,f_,g_,temp='temp);
x = variable(f);
if( v==x, return( gcdext(f,g) ) );
f_= subst(subst(f, x, temp), v, x );
g_= subst(subst(g, x, temp), v, x);
apply(t->subst(subst(t, x, v), temp, x), gcdext(f_,g_));
}
The number of bugs added to PARI when trying to work around variable priority
is staggering.
Exercise: find a bug in your code :)
Cheers,
Bill