Subject: Discussion related to cado-nfs
List archive
- From: Markus Tervooren <mail@markustervooren.de>
- To: Paul Zimmermann <Paul.Zimmermann@inria.fr>
- Cc: cado-nfs@inria.fr
- Subject: Re: [cado-nfs] Issue with "purge" on C210 GNFS
- Date: Fri, 4 Oct 2024 15:23:09 +0200
- Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None smtp.pra=mail@markustervooren.de; spf=Pass smtp.mailfrom=mail@markustervooren.de; spf=Pass smtp.helo=postmaster@sxb1plsmtpa01-03.prod.sxb1.secureserver.net
- Ironport-data: A9a23:4O9KO67JuGpCppKdbSmtHwxRtNbMchMFZxGqfqrLsTDasI4TYg02e lBvJinTZJ/9Ag22Io4oOcnztx8Y6NWG/mJQOFQ5qiFnRnlG85PMXI3Jdxb+NXLIcceeEhxt4 51DZoTNIsk9HyOB+hqnPOm+oSd2j/rQGOKsArPKZysqTAU9RCtJZX6P4wIcqtcAbY+RU1/V0 T+Ln/DiBbOF59JVGmga4v7Y+E407a37421D4wJlaaxG4FaPmSceAc5Geq3oBnapGYM88sySF 7+Slunhlo/6105wUY39yu6TnmkiGOO60d2m0yIOM0SaqkEe4HZ0iM7XDdJEAW9PkTKFgttt/ 9tEsJ20WG8BM7bF8Agne0Aw/xpWY+sfodcrHVDl6ZbPlhWeLyO1qxlTJBhe0bMwqr4f7V5mr qRwxAAlNnirm++wybSnfehg7uxLwB7DYevzElk5pd3oJa5OraLrG80m1vcEtNsEvf2iKN6FD yYvhZWDWzybC/FHEg9/5JvTB45EjFGnG9FTgAr9SabafwE/ZeG+uVTgGIO9RzCEeSlatnSH+ kTtxG6oOBgfFsSAyzOlzSmAuMaayEsXWKpKfFG53uV6igSLwncWFAUfDQL++aL/klS5V5REO 0ZNvDEjprUo+VemVMfwRVu+rWLsUhw0AooAVbRgsUfQkeyOvG51BUBcJtJFQNkvtcseQDUw1 UKSg97zGztoqrCTD36QnluRhWLiY3RNfDNaDcMCZQUo/4Pyj748sgDGEIpRSrKYtJ6oBz6ll lhmqwBk2+VK0abnzZ6T9lnLhHeqqIPhVR8w/gyRX2S/7wo/apTNWmCzwUPL569YKJ2cUkGI5 CFBw5PY9PoJAdeWhCbUBv4AF6yz6u2IKiPRmxhpEoVJGymRF2CLT58MxhpFLhZTO8MmYSDra W/s4i9vzcoGVJe1VpNfb4W0AsUs6KHvE9X5S/zZBuaihLAuKGdrGwkxOCatM3DRraQ6rU0oE bG/GftA4F4BGa4+1D6tQvsF2Odyg3hlg3vJQpy91AStiPyFaXWIU7oePUGQY/p/6q6YyOk0z zq9H5XSo/m8eLSvCsUyzWL1BQtURZTcLcuow/G7jsbZfmJb9JgJUpc9O48JdY1/hLhynezV5 Hy7UUIw4AOg3yKWeV3WMys4Mu6HsXNDQZQTYnFE0bGAhCZLXGpTxP1ALsNfkUQPqbA8naIsJ xX7U5nbXa8eItg4x9jtRcOh/NU+HPharQKJJDbtbSUleIZtDwPM8N+MQ+cc3HZmM8ZDjuNn+ +fI/lqCGfIrHl0+ZO6IM6jH5w3q4hAgdBdaAhGgzi97Ihu8qeCH6kXZ05cKHi37AUmcn2rDj VfHUH/1Z4Dl+ucIzTUAvojcx6/BLge0NhAy87DztO7vbXvp7SC4zJVeUe2FWznYWSmmsO+hf OhZhbW0evEOgF8A4cI2HqdJ3JAOwYLlh4Zb6QB4Q1TNTVCgUY17Lle8gMJgi6xqx51igzWQZ H6hwNdhFIusBNLEC38UfQosUfSC36obmx7U9vUEH3/57y5WopuBdxljPTeTggx7MLdVGoMA7 co9neIz9gefqAQYGYuUgw9150CKCHgxCJs2hKE8BITtuA4gywxfU6zlCwvdwpKGW/NTOGYEf x6Wg6vjgexH50zgKnAcK1nE7dB/t78v5i9Y7QYlCQySu9zngvQX4kVgwQ4vRF4I8iQdgvNBB Gd7EmZUe4OMxm5Mr+peVTmOHwphOkWoynbpwQFUqFyDHliabU2THmgTIu3XwVs49VhbdT1l/ L217mbpfDLpXcPp1BsJRk9XhK3/fON16zH9tpirL+adE7k+RAjVsKukSG4LihngWOcapknMo 8t09+dRN4z/EwMtoJMANorL7oRIFSi4J1FDT8pxo4IPP2XXIw+p1Ra0dkueR8JqJt7xy3GeN fBAHMx0ei6F5H68lQxDXa8oCJ1orcEt/+sHK+/KJ3ZZkr6xrQhJkZP39writV8JWdlet54PG qGLfWjTA16v2WZmtFLMiOJmOWOIR8YOSyOh/eKy8cQPT4kitsM1e24M87KEhVenGyo5wEvMp yLFRavd79I6+LRWh4G2T5lyXVSlG+38RMGj0V6Vse0XSfjtLM2XlQceime/DjRsJbFLBuhGz +WcguXWgnHAkq09CV3Cup+7EKJM28W+ccxXPu/zL1hYhSGyY9Dt0TRS51GHLYF1r/0F6vmFX weYbO6CRewRUfpZx1xXbHF6ODQZAKLVcKzhhH2cq9KhNxsj6jHEfeiXrSLRUWJmdyEzK8LfD C3wsK2Q/dx2ltlHKyIFIPBEOKVGBmHfd5EoTeCsignAPFKU2guDnpDAiSsf7Sr6DyjYMcTiv rPAaBvMVDWzn6Dq0OBpkZZ+jjtPPkZf2+NqL1MvxIJruR+HDUotD+cUAbMZAL56zw3w05DZY mnWTW0AUC/SYxVNQS/e0v/CADjGXvcvPPX9LRwXp3KkUT+8Xt68MeEw5xVe7GdTURq96uOec PU12GD6Zzq1ybFXHdcj3OSx27pb96mL10Az2B7PluLpCEwjGpQM7ntqGTRNWQHhE83gkEbqJ 3A/dVtbQXOUGFLALsJ9R0F7QB0pnivj7zENXxe9xNzyv4a6zuoZ7NbdP+r18KMIbeVUBbooa E70eVCw4DGt6iROgZcqhtMnv/YlQ7bDVM23N7TqSgAujrm9oDZvdd8LmS0UCtou4khDGlfaj SOh+GU6GF/DEk1KxbmK0k8czvqdiJ7X4+3h12YTZAMqkCDVC/DCaxLz1AbgLI3oovC6+h0DB i8Ia0rXuE2c7X3/qDJjrPkAq0aZDNtXE3TYOsztZo2niQ+iEQexC5g4u3zWFfoImJGH+mmQX KOYydYg456qRD2I2BrpmpEVb9CLap8fTx8bheaBEA1DT7ONOUCSTzp9ArOzVAUPmg17Q0oG1 dfebRrTBGo6dOWXsiNN46XsdnC31lh0O3/UJrny4QC2d4JwCqQipqM=
- Ironport-hdrordr: A9a23:gFW4Cait/7Q64zw+YEehlqPBYHBQXvoji2hC6mlwRA09T+Wxi9 2plu8KzAKcskd3ZJjvo6HiBEDyewK4yXcV2/hoAV7MZmbbUQeTXeVfBM7ZskTd8k7FmNK1t5 0NT0EWMrSZZmSSz/yKgjVQeOxQpuVvnprJuQ6T9RZQpZgBUctdBiRCe32m+pMafmR7OaY=
- Ironport-phdr: A9a23:9fsVXxweoqJKwE7XCzJDx1BlVkEcU1XcAAcZ59Idhq5Udez7ptK+Z xSZv682xwGQFazgqNt6yMPu8JrcEVQa5piAtH1QOLdtbDQizfssogo7HcSeAlf6JvO5JwYzH cBFSUM3tyrjaRsdF8nxfUDdrWOv5jAOBBr/KRB1JuPoEYLOksi7ze+/9pPPbwlSmjaxf71/I Bq3oAnMq8IbnZZsJqEtxxTGpXdFZ/5YyWR0K1yNgh3y/N2w/Jlt8yRRv/Iu6ctNWrjkcqo7U LJVEi0oP3g668P3uxbDSxCP5mYHXWUNjhVIGQnF4wrkUZr3ryD3q/By2CiePc3xULA0RTGv5 LplRRP0lCsKMSMy/XrJgcJskq1UvBOhpwR+w4HKZoGVKOF+db7Zcd8DWGZNQtpdWylHD4ihb YUAEvABMP5boYfzqFQBowWwCBeiBO3v1zFHhGP50LY10+g9CwzLxhAsE8wMvXnSsd77NL0SU eewzKTQyjvMdelW2TTg6IjTbxsvv++DXahuccrWzkkvDR3KjlONooHiOjOazOUNs2yA4uplT +2vkW0nqxxqojiz2MgskJfJhoUSylDe6yp5xpw5KsCmR0NiZ9OvDZRfuT2AOYRsXsMiX39nu Dw8yrAep5O2YjUHxIg5yxLCb/GJfZWF7xzhWeuVPTt1inJodbaiihuv8UatzunxW8u73VpUr ydIj8XAu2wR2hHP5cWJRP1w9Vqi1zaXzw3f9+JJLEMumabGNZIt36Q8m5odvEjZACP7ml36g ayVe0k+5+Sl6/rrbq/4qpOAOIJ4kB/yPrggl8G5HO81LAwOX3WA+eSgyrLj51D3QLRLj/Ioi qTVrJbXLtkBqKGjGQ9ayIMj5g6/Dzi41NQYmmEKLElKeR2ckoTpI1XOIOv9Dfejm1Sjiitny +3fMrH9DJjBNGDDkK3mfblg8UJc0wUywstC655OEbEOPu78WlLrtNPGFR82LRa0w/j5B9Vnz I8eX3yAArOBPKPTt1+H+P4vLvGRaIIWuzvxMeUp6+DzgXMjmVIQc7Ol0YUVZXygG/RpOUSZY X7igtcbFmcKuxIzTOnwh12DVT5efHOyULgm5jE4FYKmFpvDSZ63gLyf2ye7GYdWaX5AClyWD HjnaZ2IW/EXZy2KOM9ujiQEVaS9S48mzRyiqRf2y6B7IerM5i0YqZXj2cBp6OLJiB4y8iV4D 8KT02GWVGx0hWMJRzou3K9lu0B9y1GD0bJ5g/NCD9BT6elJAU8GMsvkxvF+Q+LzXwfMZNKAT lDuFu2vHDV3fNM0xtIUZk97Hf2jiArC1mykGelGuaaMAckR+7jfwmO5G8tn1XHBzqQnxw0jS 8pLHWivnK1i7BTeHZDOlF+UkeCmePJPj2b26G6fwD/W7wljWwlqXPCABChHDqO3hdHw50eZC qSrFaxiKQxZj8iLNqpNbNTty1RAXvbqftrEMCqqg2nlIxGOy/uXaZbyPX0H1XDYAUkAuwUa5 3WbLhQzHDmmpHjTCnpiGAGneFvipNF3s2jzVUoo10ePZkxl2aCy/0sQjP2YY/YaxLEZozsss StzG0yw2ZTaBonIvBJvKZ1Vesh1+1JbzSTZugh6a4SnNLxnj0UCfh5fpFPkjg1zEYBcicFy9 TZyl0xqM6WUlkhRfWrexp37J6faMWTo4h21Lajbxjkyyf6w/aECoLQ9olTn50SyE1Y6tm5gy 59T2med4ZPDCEwTV4jwWwA57UoyobaSeSQ7647OsB8keaCprj/P3c4oD+o5216hedlYKqaNC A71FYUTGcGvLOUgn1XhYAgDOahe86s9PsXucPXjuubjOe9hmBqjhH5E/Zth31iS8C1mR+iO0 5tEi/CU0w2bVivt2U+7u5OS+8gMbjUTE2yjjCn8UdcBNusrJ99NVDzoeZLko7c2z4TgUHNZ6 lO5UlYP2cvyPAGXc0S4xgpIk0IevX2gnyK8iT1yiTAg6KSFj0msi6zvcgQKPmlTSSxsl1Dpd MKxhtEedE2hdQs0jweo/1j9wLJWouJzIiOAJCUANzizNGxkXqaq4/CHas9LwJYhqyZKS/ixe 06dTaTwpV0W3mmwegkWjCB+fDastJLjmhV8g2/IN3d/ok3SfsRozAve7tjRLRJI9gIPXzIwy TzeB1znesKs4c3RjJDb9Oa3S2OmUJRXNyjt14KJ8iWhtyVmBhi2nvb7ndOCc0By1Cb60/FvX D/LsQ3gb5P10KOgPOQhckQgCFLn6sV8E51zicNp1MFWgCJDwMjOuyBY2W7oec1Ww6f/cGYAS Vtpi5bO7Q7p1VciZnOFyoTlV2mMl85oZt20eGQTiWo268FHDrvR7aQRxHsl5APi9ESMOr4nw WR4q7Nm8nMRjuAXtRB4yyycButXBkxEJWn3kAzO6dmiraJRbWLpcL6q1UM4k8ryadPK6gxaR nv9fY8vWCFq6cAqelfF0Xnb74j+e8TLcd8JqBqTiRrGyeRYYsFU9LJClW98NGTxsGdwgesyj BhG2Jigv5WcN2h34qW+GBNecDH4LZB2mHmlneNVmcCY2JqqF5NqF2AQXZfmevmvFSobqfXtM wvdWC15sHqQHqDTWBOO8Eoz5WyaCIilbjvEQRtRhcUnXhSWI1ZTxRwZTClv1IBsDRiknYTga Bsrt29XuQK+80ARjLovbUW3U3+D9ljyLG5qF97Fd0oQtkYYuyK3eYSf9r4hQHsep8P+6lbWd irDIF4URWARBh7bXRa6YObov4CQtbDfXLb2LuOSM+ve77UBB7HSmMjpisw8o1Pufo2OJiUwU a19gxUTGyogQIKB3G9TAy0Py3CXNpDd+kb6o3Ym6Jj4qqijWRqzt9HXW/0LaYQpq0jw2eDZa 4vyzG54MWoKh8lKnyeUjuhGmgZV0X0mdiHxQ+1R5WiSFeSKwPMRV1lCOmtyLJcatvJgmFATa YiB0I+zjvki0ptXQx9ETQCzw5H1I5VUcyfnbAuBXAHRaPyHPWGZmpGrJ/riGfsJ17wS71rq5 VP5WwfiJmjRzmm1EUnxYacV1nHdZFsE5Mm8ako/WDKlFoq6LEPmdocw121+wKVo1CqabyhBa X4lKBkL9vrKsUY6yr1+AzATtyU+a7PcwmDAvrOectFM6rNqGngmyLgcuSpljeILvGccA6Yp0 CrK8ow/+Rf8yrXJlWMhCF0U9lMpzMqKpRkwYP6DsMAcCTCdpk5LtSLKV1wLv4c3VYep4v4Pj IGQzuSrb20eu9PMoZlFV5iScZjedid7d0KuQmKxbkNNTCb3ZziP1woEyqjUqybT9t9g9tDtg MZcF+UdBA1zSK5cUx8/VN0afMUqAnV+yOLd1pRSoyXgyXuZDMRC4sKaCrTLUampeGzf0+QMZ gNUk+6ia95NcNSlhwo7NxF7hNiYRRqWBI8d5HQxN0ls5x8RlRo2BmwrhxC4MkXzuzlDU6bzx ER+ixMiM751pHG1uwZxfwuW4nN3yhVU+52tgDaadCP9IfWrRY8QEzDzq0U6LpL8RUBycBG2m ktncjzDQvoUhr9kPwiHkSfksIBUUb5ZRKxAO1oLwO2PIu4vyRJaozmmwklO4a3ED4FjnU0ka 8zko3VF0gNlJNk7QM6YbLJO1URVj7mSszWAzfA7mRcaPEMR6GrKIWtQ4woQLL4qYTG49LUq8 gGJhiFOY2UWSfcw5PRn6iZfc6yMwjnh3LhKNk2qf7XPf+XA4TGGyJXOGAx431hAj0Rf+Llqz covO1GZUUwi1vrZFhgEM9bDNRAAb8dW8yu2H27GuuHMzJRpeoSlQ7m4HanX7v9S3xniRll6e uZEpt4MFZStzkzCeMLuLbpejA4o+ByuP1KdSvJAZBOMlj4D5cC51p5+m4dHdVR/SS1wNzu64 rHPq0okmv2GCZ0za3Yedo4NLHErQ9G3gTBeuG5NCn+736hKrWrKpy+5vSnWADTmOpB7Y+yIY Bp3FNys0S4v+fGujUXS7o3Xfzi/bY8kosXI5aUBt5zeTehZRqRnvlzRgZdVXTqsXnKFQrvXb 9DgLoIra9LzEHOzVFey3ik0Q8nGN9GoNqGUgAvsSNUcoMyB0TskL8P4CiAGFkI6ubQY/KwlL 15mAdJzcVvyug85LaD6PAqIzoDkXTO2MTUPB/hHkbfrOPoIk2x1NbL8kDx6Etk717XlrhRLH c9VyEuFlbD7INMCNEq7Un1FJ1eQ/2xgyTcnbL1uhL94mk+AsEFAYW/RMrU3Nyoe5oF6XRTIf T13EjZqGQ/a1NWTpFTxhfZKuHEB+rQcmexd7Cqg7tmGP2LqAvXw782S6XZobMB68fcob8q9c 43f78KYxnuFENHRqlHXCnbqUaMCy54AeGQAB6MU/ANtcc0e5dgQswxrDpt4d+UJUO52+virc WY2VH9Ci3BAEdnagHpaxb7tk7rCykXJIMhkaUJU9s8Y3p1FFHQqB0FW7K67C9ePxzPCEzZbZ l1LtUIUo1hc3o5oIrK8sNaOFcEcjWcN+rQsFXGXc/sgv1rjFjPP0QS+Eavx1bXzjUQKi6iwm tgDBEwuUw4Hn7oQyBZucPYte+ERptKY62/SMx6i+ju2kq3+YwMNrK+cP1zgUNic6jG6CHVFv yRFH8kXjyuEXZUKz1gmMuBy+goKed72PB65vWBsxpw1TejiE5rzmxB69DBeAH32drgJQ+B+7 ACOA2cjOsz77sy/Z9MKBTUMsJyF9wUGygM0bXH/kMUGbZkKuGJEXSAT82zN7J3rEpQFhJUwV 9hWfrId8z/8AP8WYcjB5SRs4/qykTmBvGpg+FaimGfuQenoHqQHoitHQV9hJnzC+BRzSbF1r iGIqAqL7Q4Rna8TB6DT3xgt5m8tQtYeXHAXiDikNwghESMA6bsccfqJNZUBC/gqOU3xa1pnS q5ghhzQuxkzxyucAWQ6txMGqX2FGVdvE3lT2vC3w3Ufsp31Y2dKDcsQPHN4PnuDc0XBxUU19 F5ebUpuRpwUUM1d9fcDx4xI88HeSEGqbyYYQBhlMQF+2v1a3QtKtEHSEcgyJQ+hb/XSqQF6Z 9mWq9SkKrL18VUe4msGmPso/fwbQ2ajhRWkGYmb99S6rseNug2VZK2hdfaxamLbQSPBlw69n /EiAoWYpkA71SJeIoN51GE4bIL/DmnRMBIAK69JfyJm
- Ironport-sdr: 66ffec40_VyYkNfQtxjfgpSO5EpQvGetyp54knTh9Al+Eh4vlP0Qv1fx L8yY7LGn/kZdmiVpv92tIE7IHVNHZc7lsVlkxjQ==
Hi Paul,
thanks you for your insights on this topic.
With a target-density of 200 the matrix looks like this:
Final matrix has N=64337429 nc=64337178 (251) W=12878433733
bwc is expected to take quite a bit longer:
krylov: N=4000 ; ETA (N=2011000): Sat Nov 9 18:15:39 2024 [1.559 s/iter]
I'll continue with the density 141 matrix this time.
Markus
Am 04.10.24 um 09:35 schrieb Paul Zimmermann:
Hi Markus,
it's been a long time I did not do such a large factorization.
Date: Thu, 3 Oct 2024 21:31:24 +0200usually we try to set lim1 larger than lim0, since the norms on the algebraic
From: Markus Tervooren <mail@markustervooren.de>
Dear Paul,
I successfully obtained a suitable matrix and have now begun the BWC step.
Sieving parameters were:
-q0 1000000 -q1 920000000 -A 31 -lim0 300000000 -lim1 300000000 -lpb0 35
-lpb1 35 -mfb0 86 -mfb1 96 -ncurves0 20 -ncurves1 20
side are larger. You could try for example lim0=200000000 and lim1=400000000.
For a total of 3B raw relations. After dup2, 1.9B relations remained.an excess of 17% is quite nice!
Including all relations up to 35 bits resulted in a significantly better
matrix compared to using only 34 bits - these were removed after sieving.
The purge process reduced the total to 249 million relations:
Step 0: only singleton removal
Sing. rem.: begin with: nrows=1917941307 ncols=1948512258
excess=-30570951 at 2883.09
Sing. rem.: iter 001: nrows=1155353090 ncols=1012751300
excess=142601790 at 4037.84
Sing. rem.: iter 002: nrows=976499437 ncols=821895684 excess=154603753
at 4521.78
Sing. rem.: iter 003: nrows=929047533 ncols=773453347 excess=155594186
at 4857.31
Sing. rem.: iter 002: nrows=249083636 ncols=249083416 excess=220 atthis corresponds to an average density of 141 non-zero elements per row.
56595.46
Sing. rem.: iter 003: No more singletons, finished at 56675.36
Matrix: total 71851914 rows 71851662 cols 10096684980 coeffs
BWC is quite efficient for larger densities. You could try target-density 170
or even 200. It would be interesting to build the corresponding matrices
(merge is quite fast), and start BWC to see if it would take less than 45
days.
krylov: N=60000 ; ETA (N=2248000): Wed Oct 30 03:07:17 2024 [1.100 s/iter]yes, feel free to propose a patch for params.c210, which dates from 7 years
Expecting factors in about 45 days. The composite is part of aliquot
sequence 162126, index 4412.
Suggestions for better parameters are always welcome.
ago.
Paul
Attachment:
OpenPGP_signature.asc
Description: OpenPGP digital signature
- Re: [cado-nfs] Issue with "purge" on C210 GNFS, Markus Tervooren, 10/03/2024
- Re: [cado-nfs] Issue with "purge" on C210 GNFS, Paul Zimmermann, 10/04/2024
- Re: [cado-nfs] Issue with "purge" on C210 GNFS, Markus Tervooren, 10/04/2024
- Re: [cado-nfs] Issue with "purge" on C210 GNFS, Paul Zimmermann, 10/04/2024
Archive powered by MHonArc 2.6.19+.