Return to gausspage.c CVS log

Up to [Net/2] / Net2 / tests / benchmarks / bench-4.1

Annotation of Net2/tests/benchmarks/bench-4.1/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: 
        !            10: float  sqrt(), log(), rnd(), cos(), gauss();
        !            11: char   *valloc();
        !            12: int    rand();
        !            13: 
        !            14: main(argc, argv)
        !            15:        char *argv[];
        !            16: {
        !            17:        register int pn, i, niter, delta;
        !            18:        register char *pages;
        !            19:        float sd = 10.0;
        !            20:        int npages = 4096, pagesize, debug = 0;
        !            21:        char *name;
        !            22: 
        !            23:        name = argv[0];
        !            24:        argc--, argv++;
        !            25: again:
        !            26:        if (argc < 1) {
        !            27: usage:
        !            28:                printf(
        !            29: "usage: %s [ -d ] [ -p #pages ] [ -s standard-deviation ] iterations\n", name);
        !            30:                exit(1);
        !            31:        }
        !            32:        if (strcmp(*argv, "-s") == 0) {
        !            33:                argc--, argv++;
        !            34:                if (argc < 1)
        !            35:                        goto usage;
        !            36:                sscanf(*argv, "%f", &sd);
        !            37:                if (sd <= 0) {
        !            38:                        printf("%s: Bad standard deviation.\n", *argv);
        !            39:                        exit(2);
        !            40:                }
        !            41:                argc--, argv++;
        !            42:                goto again;
        !            43:        }
        !            44:        if (strcmp(*argv, "-p") == 0) {
        !            45:                argc--, argv++;
        !            46:                if (argc < 1)
        !            47:                        goto usage;
        !            48:                npages = atoi(*argv);
        !            49:                if (npages <= 0) {
        !            50:                        printf("%s: Bad page count.\n", *argv);
        !            51:                        exit(2);
        !            52:                }
        !            53:                argc--, argv++;
        !            54:                goto again;
        !            55:        }
        !            56:        if (strcmp(*argv, "-d") == 0) {
        !            57:                argc--, argv++;
        !            58:                debug++;
        !            59:                goto again;
        !            60:        }
        !            61:        niter = atoi(*argv);
        !            62:        pagesize = getpagesize();
        !            63:        pages = valloc(npages*pagesize);
        !            64:        if (pages == (char *)0) {
        !            65:                printf("Can't allocate %d pages (%2.1f megabytes).\n",
        !            66:                    npages, (npages*pagesize) / (1024. * 1024.));
        !            67:                exit(3);
        !            68:        }
        !            69:        pn = 0;
        !            70:        for (i = 0; i < niter; i++) {
        !            71:                delta = gauss(sd, 0.0);
        !            72:                while (pn + delta < 0 || pn + delta > npages)
        !            73:                        delta = gauss(sd, 0.0);
        !            74:                pn += delta;
        !            75:                if (debug)
        !            76:                        printf("touch page %d\n", pn);
        !            77:                else
        !            78:                        pages[pn * pagesize] = 1;
        !            79:        }
        !            80: }
        !            81: 
        !            82: float
        !            83: gauss(sd, mean)
        !            84:        float sd, mean;
        !            85: {
        !            86:        register float qa, qb;
        !            87: 
        !            88:        qa = sqrt(log(rnd()) * -2.0);
        !            89:        qb = 3.14159 * rnd();
        !            90:        return (qa * cos(qb) * sd + mean);
        !            91: }
        !            92: 
        !            93: float
        !            94: rnd()
        !            95: {
        !            96:        static int seed = 1;
        !            97:        static int biggest = 0x7fffffff;
        !            98: 
        !            99:        return ((float)rand(seed) / (float)biggest);
        !           100: }

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