ggnfs improvements

[bsquared]

Thank you very much!

Found and fixed, uploaded just now.

time ./gnfs-lasieve4I12e -k -v -n0 -a c103-2.job.T1 -o test.out
gnfs-lasieve4I12e (with asm64,avx-512 mmx-td,lasetup,lasched,sieve1,ecm,tds0,search0,tdsched): L1_BITS=15
FBsize 100777+0 (deg 4), 171534+0 (deg 1)
Sorted factor base on side 0: 1: 32495 2: 32397
Sorted factor base on side 1: 1: 168023
total yield: 155958, q=1327517 (0.00051 sec/rel) ETA 0h00m)
937 Special q, 1411 reduction iterations

reports: 165365787->20138278->18040597->11058350->11056552->10304931->321750 (5306059)

Total yield: 155958
milliseconds total: Sieve 24660 Sched 8950 medsched 6040
TD 25660 (Init 970, MPQS 3550) Sieve-Change 13800
TD side 0: init/small/medium/large/search: 790 2510 930 1140 6040
sieve: init/small/medium/large/search: 1470 10160 750 860 1170
TD side 1: init/small/medium/large/search: 130 2460 1010 810 4910
sieve: init/small/medium/large/search: 640 6140 990 1180 1300
aborts: 0 0
Expected yield/cost: 4.15e+04  0
p-1: 0 tests, 0 successes  ecm: 0 tests, 0 successes
MPQS-AUX 0
COF: 280255 tests, 0 ecm, 0 aux:
       0 mpqs, 0 mpqs3, 0 ecm, 0 too big
77.763u 1.444s 1:19.53 99.5%    0+0k 1648+23208io 1pf+0w

[chris2be8]

Why does c103-2.job.T1 contain:

m:

My version of factMsieve.pl contains (shortly after "sub makeJobFile"):

print OUTF "m: $M\n" if ($M);

The original version contains:

print OUTF "n: $N\nm: $M\n";

So I think that's a bug I fixed some time ago. Try removing m: from the job file and re-testing to be sure if that is the issue.

[bsquared]

Originally Posted by chris2be8

Why does c103-2.job.T1 contain:


My version of factMsieve.pl contains (shortly after "sub makeJobFile"):

print OUTF "m: $M\n" if ($M);

The original version contains:

print OUTF "n: $N\nm: $M\n";

So I think that's a bug I fixed some time ago. Try removing m: from the job file and re-testing to be sure if that is the issue.

That could be a separate issue, but there really was a problem that I fixed.

It was an inexact translation of #defines like this to AVX512 in lasieve_prepn:
#define A1MOD0(p) ((aux= absa1%p)> 0 ? p-aux : 0 )

where I neglected to account for the "else" case (where absa1%p == 0)

[wreck]

When sieve R942 using lpb34, the avx512-lasieve5 crash immediately,
while the original lasieve5 could work properly.

command is
./gnfs-lasieve4I16e -v -f 380000000 -c 1000 -o R942_16e_r_380000000_380001000.out -r R942_poly.txt -R

R942_poly.txt file attached.

[bsquared]

Originally Posted by wreck

When sieve R942 using lpb34, the avx512-lasieve5 crash immediately,
while the original lasieve5 could work properly.

command is
./gnfs-lasieve4I16e -v -f 380000000 -c 1000 -o R942_16e_r_380000000_380001000.out -r R942_poly.txt -R

R942_poly.txt file attached.

Fixed and checked in! Thanks for the report!

[wreck]

Thanks for the fix.

Another confusion is that ,
when set the rlim and alim to 500000000,
the version 551 (newest version) is about 10% slower than the 550 version (second newest version).

I'm not sure if change the _mm512_set1_epi32(ij_ub) to _mm512_set1_epu32(ij_ub)
or such modification could resolve the slower speed.