Return to factor.c CVS log

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

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

1.1       root        1: #include <math.h>
                      2: 
                      3: extern double modf();
                      4: 
                      5: #ifdef pdp11
                      6: #define setfps() ldfps(0240)   /* double, chop */
                      7: #define ifdiv(d,n) if((n/d)*d==n)
                      8: 
                      9: #else
                     10: #define setfps()
                     11: 
                     12: #ifdef vax
                     13: #define ifdiv(d,n) junk = n/d;\
                     14:                   asm(" emodd _junk,$0,$0f1.0,r0,_junk");\
                     15:                   if(junk==0.)
                     16: #else
                     17: #define ifdiv(d,n) if(modf(n/d,&junk)==0.)
                     18: #endif
                     19: 
                     20: #endif
                     21: 
                     22: double wheel[] = {
                     23:         2,10, 2, 4, 2, 4, 6, 2, 6, 4,
                     24:         2, 4, 6, 6, 2, 6, 4, 2, 6, 4,
                     25:         6, 8, 4, 2, 4, 2, 4, 8, 6, 4,
                     26:         6, 2, 4, 6, 2, 6, 6, 4, 2, 4,
                     27:         6, 2, 6, 4, 2, 4, 2,10,
                     28: };
                     29: double n;
                     30: double d;
                     31: double junk;
                     32: double s;
                     33: 
                     34: main(argc,argv)
                     35: char **argv;
                     36: {
                     37:        extern double atof();
                     38:        register double *wheelp;
                     39:        if(argc>1)
                     40:                n = atof(argv[1]);
                     41:        setfps();
                     42:        do {
                     43:                if(argc<=1 && scanf("%lf",&n)<=0)
                     44:                        break;
                     45:                if(n<=0)
                     46:                        break;
                     47:                printf("%.0f\n",n);
                     48:                s = sqrt(n) + 1;
                     49:                try(2.);
                     50:                try(3.);
                     51:                try(5.);
                     52:                try(7.);
                     53:                d = 1.;
                     54:                wheelp = &wheel[1];
                     55:                for(;;) {
                     56:                        d += *wheelp;
                     57:                        ifdiv(d,n)
                     58:                                try(d);
                     59:                        if(++wheelp >= wheel + sizeof(wheel)/sizeof(*wheel)) {
                     60:                                wheelp = wheel;
                     61:                                if(d>s)
                     62:                                        break;
                     63:                        }
                     64:                }
                     65:                if(n>1)
                     66:                        printf("     %.0f\n",n);
                     67:                printf("\n");
                     68:        } while(argc<=1);
                     69:        return(0);
                     70: }
                     71: 
                     72: try(p)
                     73: double p;
                     74: {
                     75:        static double quot;
                     76:        while(modf(n/p,&quot)==0) {
                     77:                printf("     %.0f\n",p);
                     78:                n = quot;
                     79:        }
                     80:        s = sqrt(n) + 1;
                     81: }

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