Attempt p7 in fortran (15)

LivInTheLookingGlass · Oct 3, 2024 · 56a1a96 · 56a1a96
1 parent 03e9d29
commit 56a1a96
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/fortran/src/include/primes.f90 b/fortran/src/include/primes.f90
@@ -11,7 +11,7 @@ subroutine sieve_up_to(n)
         integer, intent(in) :: n
         integer :: p, i, limit
 
-        limit = int(sqrt(real(n)))
+        limit = ceiling(sqrt(real(n)))
         do p = 2, limit
             if (get_prime_bit(p)) then
                 do i = p * p, n, p