From 8fdb2354b50bedf84e0744346e974d308ded6a1b Mon Sep 17 00:00:00 2001 From: jakllsch Date: Tue, 30 Aug 2011 02:58:04 +0000 Subject: No need to cast double to ubig (aka uintmax_t) through int. This change prevents a modulo by zero in a invocation such as: primes 18446744073709551000 18446744073709551615 on a LP64 machine. --- primes/primes.c | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) (limited to 'primes') diff --git a/primes/primes.c b/primes/primes.c index c8c6ea58..5952bbf6 100644 --- a/primes/primes.c +++ b/primes/primes.c @@ -1,4 +1,4 @@ -/* $NetBSD: primes.c,v 1.18 2010/05/13 17:52:12 tnozaki Exp $ */ +/* $NetBSD: primes.c,v 1.19 2011/08/30 02:58:04 jakllsch Exp $ */ /* * Copyright (c) 1989, 1993 @@ -42,7 +42,7 @@ __COPYRIGHT("@(#) Copyright (c) 1989, 1993\ #if 0 static char sccsid[] = "@(#)primes.c 8.5 (Berkeley) 5/10/95"; #else -__RCSID("$NetBSD: primes.c,v 1.18 2010/05/13 17:52:12 tnozaki Exp $"); +__RCSID("$NetBSD: primes.c,v 1.19 2011/08/30 02:58:04 jakllsch Exp $"); #endif #endif /* not lint */ @@ -305,11 +305,10 @@ primes(ubig start, ubig stop) /* note highest useful factor and sieve spot */ if (stop-start > TABSIZE+TABSIZE) { tab_lim = &table[TABSIZE]; /* sieve it all */ - fact_lim = (int)sqrt( - (double)(start)+TABSIZE+TABSIZE+1.0); + fact_lim = sqrt((double)(start)+TABSIZE+TABSIZE+1.0); } else { tab_lim = &table[(stop-start)/2]; /* partial sieve */ - fact_lim = (int)sqrt((double)(stop)+1.0); + fact_lim = sqrt((double)(stop)+1.0); } /* sieve for factors >= 17 */ factor = 17; /* 17 is first prime to use */ -- cgit v1.2.3-56-ge451