|
Bill Allombert on Mon, 05 Jun 2023 18:23:18 +0200
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
|
Re: Abnormous memory use for gaussian gcd()?
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: Re: Abnormous memory use for gaussian gcd()?
- From: Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>
- Date: Mon, 5 Jun 2023 18:18:27 +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=1685981897; c=relaxed/relaxed; bh=KU+sdAwplNnxvfHr4Twisy6+mOSKhyc3uat/cHwzrTU=; h=DKIM-Signature:Date:From:To:Subject:Message-ID:Mail-Followup-To: References:MIME-Version:Content-Type:Content-Disposition: In-Reply-To; b=OAbNk8nEAYTsGavSLTKTII8p2BcMxRWByd+eq8PF19JKP3nchArNSeT1lyIO3EHI+gOQj7iCsU4AhG4IoN1IQnYl5OF4HLcJ2t0yUTs7CFuUTFwPEwvurwlAHmCgCijV3WNCVeqWjQs5gZX4gaKWph3FGDDuwyv5QSAJWf5L1pmt5QkDlVy2dueVJIFJJpaPM7ccNXCqQ9whR1t+6hgQMK5vXF6hq6rbB7IHa6UAuHXuGW/usMVnKE9UUch20G7Mon4vkbZqme8WPFxpdGOrXO7xYCl2eXa6lBNZjKhT1BmFl0Y+eh3ICOINkcgDWNwSTx5tnPR8niRgt2FSz7lI36LQ7TjXmaPMgomcS4ZxrLU1XkdBN3TzT0FhL2hGG6WQJ1pyA3NoDeTc6ErqktoJ3piJA6Cp0l8WhiYD9GNNea/mJFY8L1Y3DdssRXNE+tDUhpPKF15Hh93oiM6RDmqVkgZqc976dvijkmiJSI1WReGDpzFafDvxjrN08xxBY7Y2O6buAQgVqaoiJ0qa+/AORTC9iFXNL6fiiIDBVkr4JrEQ3iocpxdZyIrb6FKQQEW0HGWHfBpjml8Dn3NEL1raQSFgo9G4sCNSnv3UAL9jZqqpBZwyjV9NpGrMjivypsQ6Z6873eTyQb7LQ/aYT9HPQ/uW8F/cTWA0xgcpM/ukLY8=
- Arc-seal: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1685981897; cv=none; b=TIyYFwPhRGASlHIBjtz53D2oknMJBfLbMXPs5ZguQBhhhYvtDXP3pmYGeDBWIiM74YQCw9of9v17t55wkmomIYF8dOfpcyhU3KC9JPhxvet0pU1vUJrf4ybuRLWmSF7gCp1jip29xwcj+aMwmEHok+R29oc3yDe6diNryX+R4dsAka4SGlNWSO/FejYY2e1+Ltr3wBrcf0j323y8tM1TdwgVGeN2Cvs7UfWIpMR7VZkxpir1ue/fn6ZqMwN6GKS9JZX2I4oUX64wjtaZpdeytp3V0+TLq+y9w0B+iDP8p1lcJUb9kwX24S3uy4SVndFr9Qotq352brfYeC4iB5V4tCQWzLgEB2AMZZfFB2LPcSTKUTP5cYMQ7f4LqsK9fjH7Cu3umMu0M/izVXMa39oPtXk3NrNfOhY8hRNILyKExrHcmcXCP9Ii4bUIh3E9WUObS9HUplCHsQOe69xstN6e99G4q3VvcI7/opVldZwayg+aSvjHTW+D6KFcSIHXfF4Jec2D1WvfsJCX3aJb2CcC5ORzZp3Vgs9hpEjuTT2VTSKfJ+Rru8GoKjYuJnWNp7sJItsp1yZc+GNX5yAQWTVOOctqBpBHymJfFSkAl1jW/gaTsfrBN2LHqGqkLLfeOmEhau9QFzABMB3L04oqaJQsmMbpcuCndWb7OTo/K0/2TU8=
- Authentication-results: smail; dmarc=none header.from=math.u-bordeaux.fr
- Authentication-results: smail; arc=none
- Delivery-date: Mon, 05 Jun 2023 18:23:18 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=math.u-bordeaux.fr; s=2022; t=1685981897; bh=KU+sdAwplNnxvfHr4Twisy6+mOSKhyc3uat/cHwzrTU=; h=Date:From:To:Subject:References:In-Reply-To:From; b=g4/OUiggRRSbRkZkmpUp2xDopbMgHBaUVbcMYS3Lv8iNlH2p8sQVa58rysUaLyh+a zCyHf3tNDPWHr4ckVI+khGTdiW4vnpvlil1qsj8OMW1Tly0TOKKHcMWRvrzwNKt6AX ADkBqL6YLPcLDt5nrhwdVw3dgh/nZt2wWrPm7v6urIhGa27yiTmG7Uy4If1Kt55Hu1 JLCYXJhvYLixe5z7iu3wsumxtCOtQqfrYit7ZMU8ss1TMo3P97QpGOQc51GWdYkyEG pkv3wJaCC95zVcTZPOmo5/qqNFQQMyGizqkMk0eTYZBnmVcuc3ho9VGJnQM55nXME3 vakTShr5R7FWq+cn7DKN0UU1+YlABEWigQz48RdkIYfXr9Mmapp9bYlGuSR07AnqYt G6n3tWjl0145qfJ75n6CRX3D9GXe/LEAzDyVmx2a2RZlC6MvwH7t/fo2DxQMUAQ52r BrHhXPyU08R4WknYAQ48wHW5Xep+caurPM7lR3/k8LBiG+V9RiSYaH+EZDT9v6REpB vzdc7c01sEQqfSPcvywk8Hve0rVMsp7aN7dXBySq2WCGcarfLHeQ+bRGsLw7OgiVg5 id7G+OLziZSgaGiRLieoczw+/s775b0tzfNT2KfzoBXXik4vOm5CsVOwRZB4SZd/h8 Z1rqcfZlDl5edQSz6bFFSDII=
- In-reply-to: <85b829f5974500a9215307f172e05724@stamm-wilbrandt.de>
- Mail-followup-to: pari-users@pari.math.u-bordeaux.fr
- References: <85b829f5974500a9215307f172e05724@stamm-wilbrandt.de>
On Mon, Jun 05, 2023 at 05:51:37PM +0200, hermann@stamm-wilbrandt.de wrote:
> Transpiled GP code for 10000-digit number works fine:
>
> $ gp 10000.gp
> ...
> 39178247477405595462999004573259238447193580865522753006504996585341802644743745733315578184178833806642718982619192126965967031226198644279480298699929676[+++]
> *** last result computed in 4,709 ms.
> ?
>
> But even though I give 24G of my laptop's 32G to Pari, computation does not
> complete.
> What am I doing wrong?
> Or is this a bug in Pari/GP?
> Python 3.9, Pari/GP 2.15.3.
>
> $ grep ^parisize /etc/gprc
> parisizemax = 24G
> parisize = 24G
> $ gp 36401.gp
> Reading GPRC: /etc/gprc
> GPRC Done.
> ...
> parisizemax = 24000000000, primelimit = 500000
> foobar
> *** The result history is empty.
When doing 'gp 36401.gp' there is no "result history" (that is "%1 = blah") so referencing it fails.
Probably 36401.gp is using ## (or %).
I suggest you try
./gp < 36401.gp
I do not see anything that suggest there is a memory error.
Cheers,
Bill.