Karim Belabas on Sat, 26 Nov 2022 12:19:38 +0100
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
- To: pari-dev@pari.math.u-bordeaux.fr
- Subject: Re: 3^k performance
- From: Karim Belabas <Karim.Belabas@math.u-bordeaux.fr>
- Date: Sat, 26 Nov 2022 12:18:33 +0100
- Arc-authentication-results: i=1; smail; arc=none
- Arc-message-signature: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1669461512; c=relaxed/relaxed; bh=CZD1iSRfLkYwbawuTo/pepSZI9X+w0Uvz0fIZRthNI4=; h=DKIM-Signature:Date:From:To:Subject:Message-ID:Mail-Followup-To: References:MIME-Version:Content-Type:Content-Disposition: Content-Transfer-Encoding:In-Reply-To; b=0PN62gHBnqxPi0Wq8n/sG36GZ4nIc8sZVzhx+eXcYY50EPqpo+kn6FT6n0zdakBte6LhH8Eq8YenoQUbk/yWlpPrxbQKPOSRO+i+kNABmMK9oHclsja+jBStxJnwvVUTl9jQoeolV8fE51S8fB1I6TXvl2TcRhMNIre/0qx6zUGi6dq0sZNlwpQixnugs9d3x/V6ia1Gnwb2DSMv6nqYxzZ8hKHgCGnEXJm6nXp6IWZqDajsmEopDbQ557FafzLx4TrdVpxexF73cvy15LvT+ajiHcauPUMC9PjQxCXDZgEz8Bc9ejJeTayBd8QIogHfgor7CUPaNwXKQ21g3pCujt8Qyxo4X2Ew3WlDY/2j9tZxmhshCRMCA7rbSwWcNcgeJNLmcE2W7lf9ts/7eu2is8ogbjuuaoyXY1qD1FuIAyn+F7r8V3wnFWmku+8qhTgbhKDdMEwH/tA2B1BYPtZLwjRAnH/QDNdxrso2uKnV/MBCxk1Z6l8+/tN/aSP8otoJqYRdzqh3HyC1rhI4CSVGkkLrwLJ3U1xbGxU0FIulBbCqVORQbkLw68dT7k3WNM0f+SbutQbLlYI0EMFMd0l4FDgYDoBjlszLhyLYcX+q45srevLnfCMhI/yTZ0YhPgA7gokmbJYk8kJor91yYZfwJahgUGu1RCRQWWcRie2hGcU=
- Arc-seal: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1669461512; cv=none; b=OXHS/m9jumjs+frTin5mG6/5fcoWe3xN0LX3BEp4h+xzCONgx3w+yJLnNm67U4r9B7TZkgiBU8nZq/zwc8LHq4mVz+24Tqj4lVLIDhXoZo+tAbmPIBK3VTnnC/aVSwTlVA4YQJ1FljBjZK6hXnSOLqc2aQYfEa6TcfmjMM/TX0eF6QiE5BfLpZXhyP9ebf5m7v9LXOOvwQ9wDYQpM2ZC9hoLWwM/UIyZCL8ljyhchLT5LbxjZZV5s7xMx+jdQiReyRrs/fKCoPIW2vwrRWAhfYbjRcB8XxwqwUjoUTxlDLVt/hi7qpbJYkRTUa2b66NiF+RItoK7hLkMkTQYnHWiDWOecPFnz9uJ94nMBfrvtFfVZxmX5woNmTWMJMyShnYK+EJuuFbX8p93i1vaW6kvPyRCKHkD1WMHNzcQSO2d5MLS7ZylIYY5mdB5oeVgjlCUizDdttfNXiRP7CeMhRbdIXdDDTBf+Y+nXNP/fAG9F15ZfpAkQUIgp0I/0oPMWI88TUYONTgX87wY1eDmbRqViKat5qEBJNYDwijMe5r8tivlcJ8BrTjS1RSG4kEpCqAudP1hlBR6sGyyu5a0jP7DcCYUW4zTBRW3KYghB1BAb2G25pRtWwoi0oLllWxTLkx20FmLavC0yWKuuWtyDLfbMK3Z2CTIMlR2G32L+r7kBIk=
- Authentication-results: smail; dmarc=none header.from=math.u-bordeaux.fr
- Authentication-results: smail; arc=none
- Delivery-date: Sat, 26 Nov 2022 12:19:38 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=math.u-bordeaux.fr; s=2022; t=1669461512; bh=CZD1iSRfLkYwbawuTo/pepSZI9X+w0Uvz0fIZRthNI4=; h=Date:From:To:Subject:References:In-Reply-To:From; b=mVl0MGBkoEAhRUwLeIl8/pUMyOKj7adqGJU0pXCTBBCJiJnFvibb1mQDxSys9pSBr u3xO9Jw7sX6D2dGK8mBopqEyKrEB2ngQ2234VAvnr6GW17IeZaAj9ZIz77+P80NT6u HcEbNvSVXWIfsvH0TcCQMhsw1amLNYLZG/lbkxCoOKoCDhq44uWRftKQfSbal7hO+d uv5C1YXFVHnaJBKUqQ6NdR3UXv11/7n3exVCnPsGc1NKICuSNTgxsYmuJ1WTyTjiAC dy72n0kPVz9Goi5UeDxWIvGC1ckOw+tpBHXuBjmsYQDAAEFECVUTvqq34iYevGakBW 44iDA8Bxfzp9G2P34sW2t45EBVuASOnG3WgginKHrP6o5lXERluJXYoN3ZXnhfI127 wnZ/cjNa9K6gxO6LpiTIAuUSplhxX5HsmNyl8RBwKvYdZ7inXw4HO28OfddVsuUQ8i 4W9wGVy08xl3EhkG9vO+OqHX8k3PXqILuiVNVrYt+3AG5EfUMfOZvN2zL7PvGpXO1v XCIJJG23cy0qHktOjLm0n8doMppbM2D81CB2PiW3EWLPjlJUyiEziG5qT78q3SPd+m 91488iOzn94Kg+6pgb9nA3JC1FheV3JY59261g5Y2cjX4VjdQ1lqvAtKBqHW+QO/SX OQygoF6MOmjT37csql4elhvI=
- In-reply-to: <Y4HRNfRGZBQgIwoT@seventeen>
- Mail-followup-to: pari-dev@pari.math.u-bordeaux.fr
- References: <de8f200c-386e-3028-f44b-519b9ce445c1@isolution.nl> <Y4HRNfRGZBQgIwoT@seventeen>
* Bill Allombert [2022-11-26 09:47]:
[...]
> The simplest approach would be to use powers(3,10^5), with to caveat:
[...]
> - the powers function has a bug and is very slow for integers.
[...]
Here's a temporary workaround.
\\ slow
(12:16) gp > v=powers(3,10^5);
time = 2,191 ms.
\\ reasonably fast
(12:16) gp > w=powers(3,10^5,1);
time = 345 ms.
(12:16) gp > v==w
%3 = 1 \\ identical
Will investigate later.
Cheers,
K.B.
--
Karim Belabas / U. Bordeaux, vice-président en charge du Numérique
Institut de Mathématiques de Bordeaux UMR 5251 - (+33) 05 40 00 29 77
http://www.math.u-bordeaux.fr/~kbelabas/
`