IMPROVED: slightly cheaper constant matching operation for chr_identifier store
[chr.git] / Examples / primes.chr
blobdf9963963f44d930cf2ef8d48b2b677d57f4a5d2
1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2 %%
3 %% Sieve of eratosthenes to compute primes
4 %% thom fruehwirth 920218-20, 980311
5 %% christian holzbaur 980207 for Sicstus CHR
6 %%
7 %% ported to hProlog by Tom Schrijvers 
9 :- module(primes,[]).
10 :- use_module(library(chr)).
12 :- constraints candidate/1.
13 :- constraints prime/1.
16 candidate(1) <=> true.
17 candidate(N) <=> primes:prime(N), N1 is N - 1, primes:candidate(N1).
19 absorb @ prime(Y) \ prime(X) <=> 0 is X mod Y | true.
21 time(N):-
22         cputime(X),
23         candidate(N),
24         cputime( Now),
25         Time is Now-X,
26         write(N-Time), nl.
28 cputime( Ts) :- 
29         statistics( runtime, [Tm,_]),
30         Ts is Tm/1000.