Attempt p3 in fortran (5)

LivInTheLookingGlass · Oct 4, 2024 · c38d1ea · c38d1ea
1 parent 4ca1aa5
commit c38d1ea
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/fortran/src/include/primes.f90 b/fortran/src/include/primes.f90
@@ -87,6 +87,7 @@ subroutine prime_factor(num, factor)
             factor = num
         else
             factor = 2
+            call expand_sieve(int(sqrt(num), i18t) + 1)
             do while (factor <= num)
                 if (mod(num, factor) == 0) then
                     num = num / factor