[Ggnfs-commits] SF.net SVN: ggnfs:[403] trunk/src/lasieve4/gnfs-lasieve4e.c

[<brg...@us...>]

Revision: 403
          http://ggnfs.svn.sourceforge.net/ggnfs/?rev=403&view=rev
Author:   brgladman
Date:     2011-01-15 09:02:27 +0000 (Sat, 15 Jan 2011)

Log Message:
-----------
correct bug in lattice siever

Modified Paths:
--------------
    trunk/src/lasieve4/gnfs-lasieve4e.c

Modified: trunk/src/lasieve4/gnfs-lasieve4e.c
===================================================================
--- trunk/src/lasieve4/gnfs-lasieve4e.c	2011-01-12 12:22:41 UTC (rev 402)
+++ trunk/src/lasieve4/gnfs-lasieve4e.c	2011-01-15 09:02:27 UTC (rev 403)
@@ -2530,16 +2530,22 @@
                         { unsigned long small_factors[10];
 
                           nf = rho_factor(small_factors, large_factors[s1]);
-                          for (i = 0; i < nf; i++) {
-                            mpz_set_ui(large_primes[s1][i], small_factors[i]);
-                            if (mpz_sizeinbase(large_primes[s1][i],2) > max_primebits[s1]) {
-                              n_mpqsvain[s1]++;
-                              break;
-                            }
+                          if(nf >= 1)
+                          {
+                              for (i = 0; i < nf; i++) {
+                                mpz_set_ui(large_primes[s1][i], small_factors[i]);
+                                if (mpz_sizeinbase(large_primes[s1][i],2) > max_primebits[s1]) {
+                                  n_mpqsvain[s1]++;
+                                  break;
+                                }
+                              }
+                              if (i >= nf)
+                                  nlp[s1] = nf;
                           }
-                          if ((i >= nf) && (nf >= 1))
-                            nlp[s1] = nf;
-                          else { nlp[s1]=0; }
+                          else
+                          { 
+                              nlp[s1] = 0; 
+                          }
                         }
 #else
                           nlp[s1] = 0;


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.