Return to primes.c CVS log

Up to [Research Unix] / researchv10no / cmd / primes

Annotation of researchv10no/cmd/primes/primes.c, revision 1.1.1.1

1.1       root        1: static double big = 7.20575940379e16;
                      2: #define ptsiz 50
                      3: #define tabsiz 1000
                      4: #define tsiz8 8000
                      5: #define whsiz 48
                      6: double nn;
                      7: char table[tabsiz];
                      8: int pt[] = {
                      9:          2,  3,  5,  7, 11, 13, 17, 19, 23, 29,
                     10:         31, 37, 41, 43, 47, 53, 59, 61, 67, 71,
                     11:         73, 79, 83, 89, 97,101,103,107,109,113,
                     12:        127,131,137,139,149,151,157,163,167,173,
                     13:        179,181,191,193,197,199,211,223,227,229,
                     14: };
                     15: double wheel[] = {
                     16:        10, 2, 4, 2, 4, 6, 2, 6, 4, 2,
                     17:         4, 6, 6, 2, 6, 4, 2, 6, 4, 6,
                     18:         8, 4, 2, 4, 2, 4, 8, 6, 4, 6,
                     19:         2, 4, 6, 2, 6, 6, 4, 2, 4, 6,
                     20:         2, 6, 4, 2, 4, 2,10, 2,
                     21: };
                     22: char bittab[] = {
                     23:        1, 2, 4, 8, 16, 32, 64, 128,
                     24: };
                     25: 
                     26: main(argc, argp)
                     27: int argc;
                     28: char *argp[];
                     29: {
                     30:        register i;
                     31:        double k, temp, v;
                     32:        double limit = big;
                     33:        double atof(), sqrt();
                     34: 
                     35:        if(argc > 1){
                     36:                nn = atof(argp[1]);
                     37:        }else{
                     38:                scanf("%lf", &nn);
                     39:        }
                     40:        if(argc > 2){
                     41:                limit = atof(argp[2]);
                     42:                if(limit < nn) exit(0);
                     43:                if(limit > big) ouch();
                     44:        }
                     45:        if((nn<0) || (nn>big)) ouch();
                     46:        if(nn == 0) exit();
                     47: 
                     48:        if(nn < 230.){
                     49:                for(i=0; i<ptsiz; i++){
                     50:                        if(pt[i] < nn) continue;
                     51:                        if(pt[i] > limit) exit(0);
                     52:                        printf("%d\n", pt[i]);
                     53:                }
                     54:        nn = 230.;
                     55:        }
                     56: 
                     57:        modf(nn/2, &temp);
                     58:        nn = 2.*temp + 1;
                     59: /*
                     60:  *     clear the sieve table.
                     61:  */
                     62:        while(1){
                     63:                for(i=0;i<tabsiz;i++)
                     64:                        table[i] = 0;
                     65: /*
                     66:  *     run the sieve.
                     67:  */
                     68:                v = sqrt(nn+tsiz8);
                     69:                mark(nn,(double)3);
                     70:                mark(nn,(double)5);
                     71:                mark(nn,(double)7);
                     72:                for(i=0,k=11;k<=v;k=k+wheel[i]){
                     73:                        mark(nn,k);
                     74:                        if(++i == whsiz) i = 0;
                     75:                }
                     76: /*
                     77:  *     now get the primes from the table
                     78:  *     and print them.
                     79:  */
                     80:                for(i=0; i<tsiz8; i += 2){
                     81:                        if(table[i>>3] & bittab[i&07]) continue;
                     82:                        temp = nn + i;
                     83:                        if(temp > limit) exit(0);
                     84:                        printf("%.0f\n", temp);
                     85:                }
                     86:                nn += tsiz8;
                     87:        }
                     88: }
                     89: 
                     90: mark(nn,k)
                     91: double nn,k;
                     92: {
                     93:        double t1;
                     94:        register j;
                     95:        register k1 = k;
                     96: 
                     97:        modf(nn/k, &t1);
                     98:        j = k*t1 - nn;
                     99:        if(j < 0)
                    100:                j += k1;
                    101:        for( ; j < tsiz8; j += k1){
                    102:                table[j>>3] |= bittab[j&07];
                    103:        }
                    104: }
                    105: 
                    106: ouch()
                    107: {
                    108:        printf("OUCH!\n");
                    109:        exit(1);
                    110: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.