Return to gausspage.c CVS log

Up to [Net/2] / Net2 / tests / benchmarks

Annotation of Net2/tests/benchmarks/gausspage.c, revision 1.1

1.1     ! root        1: /*
        !             2:  * Random page access with
        !             3:  * a gaussian distribution.
        !             4:  *
        !             5:  * Allocate a large (zero fill on demand) address
        !             6:  * space and fault the pages in a random gaussian
        !             7:  * order.
        !             8:  */
        !             9: #include <math.h>
        !            10: 
        !            11: float  rnd(), gauss();
        !            12: char   *valloc();
        !            13: int    rand();
        !            14: 
        !            15: main(argc, argv)
        !            16:        char *argv[];
        !            17: {
        !            18:        register int pn, i, niter, delta;
        !            19:        register char *pages;
        !            20:        float sd = 10.0;
        !            21:        int npages = 4096, pagesize, debug = 0;
        !            22:        char *name;
        !            23: 
        !            24:        name = argv[0];
        !            25:        argc--, argv++;
        !            26: again:
        !            27:        if (argc < 1) {
        !            28: usage:
        !            29:                printf(
        !            30: "usage: %s [ -d ] [ -k #Kb ] [ -s standard-deviation ] iterations\n", name);
        !            31:                exit(1);
        !            32:        }
        !            33:        if (strcmp(*argv, "-s") == 0) {
        !            34:                argc--, argv++;
        !            35:                if (argc < 1)
        !            36:                        goto usage;
        !            37:                sscanf(*argv, "%f", &sd);
        !            38:                if (sd <= 0) {
        !            39:                        printf("%s: Bad standard deviation.\n", *argv);
        !            40:                        exit(2);
        !            41:                }
        !            42:                argc--, argv++;
        !            43:                goto again;
        !            44:        }
        !            45:        if (strcmp(*argv, "-k") == 0) {
        !            46:                argc--, argv++;
        !            47:                if (argc < 1)
        !            48:                        goto usage;
        !            49:                npages = atoi(*argv);
        !            50:                if (npages <= 0) {
        !            51:                        printf("%s: Bad virtual memory size.\n", *argv);
        !            52:                        exit(2);
        !            53:                }
        !            54:                argc--, argv++;
        !            55:                goto again;
        !            56:        }
        !            57:        if (strcmp(*argv, "-d") == 0) {
        !            58:                argc--, argv++;
        !            59:                debug++;
        !            60:                goto again;
        !            61:        }
        !            62:        niter = atoi(*argv);
        !            63:        pagesize = getpagesize();
        !            64:        npages /= pagesize / 1024;
        !            65:        pages = valloc(npages, pagesize);
        !            66:        if (pages == (char *)0) {
        !            67:                printf("Can't allocate %d pages (%2.1f megabytes).\n",
        !            68:                    npages, (npages*pagesize) / (1024. * 1024.));
        !            69:                exit(3);
        !            70:        }
        !            71:        if (pagesize != 1024)
        !            72:                printf("Pagesize %dKb\n", pagesize / 1024);
        !            73:        pn = 0;
        !            74:        for (i = 0; i < niter; i++) {
        !            75:                delta = gauss(sd, 0.0);
        !            76:                while (pn + delta < 0 || pn + delta > npages)
        !            77:                        delta = gauss(sd, 0.0);
        !            78:                pn += delta;
        !            79:                if (debug)
        !            80:                        printf("touch page %d\n", pn);
        !            81:                else
        !            82:                        pages[pn * pagesize] = 1;
        !            83:        }
        !            84: }
        !            85: 
        !            86: float
        !            87: gauss(sd, mean)
        !            88:        float sd, mean;
        !            89: {
        !            90:        register float qa, qb;
        !            91: 
        !            92:        qa = sqrt(log(rnd()) * -2.0);
        !            93:        qb = 3.14159 * rnd();
        !            94:        return (qa * cos(qb) * sd + mean);
        !            95: }
        !            96: 
        !            97: float
        !            98: rnd()
        !            99: {
        !           100:        static int seed = 1;
        !           101:        static int biggest = 0x7fffffff;
        !           102: 
        !           103:        return ((float)rand(seed) / (float)biggest);
        !           104: }

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