Return to sim.c CVS log

Up to [Research Unix] / researchv9 / cmd / cfront / libC / task

Annotation of researchv9/cmd/cfront/libC/task/sim.c, revision 1.1

1.1     ! root        1: #include "task.h"
        !             2: 
        !             3: histogram.histogram(int nb, int ll, int rr)
        !             4: {
        !             5:        register int i;
        !             6:        if (rr<=ll || nb<1) task_error(E_HISTO,0);
        !             7:        if (nb%2) nb++;
        !             8:        while ((rr-ll)%nb) rr++;
        !             9:        binsize = (rr-ll)/nb;
        !            10:        h = new int[nb];
        !            11:        while (h == 0) task_error(E_STORE,0);
        !            12:        for (i=0; i<nb; i++) h[i] = 0;
        !            13:        l = ll;
        !            14:        r = rr;
        !            15:        nbin = nb;
        !            16:        sum = 0;
        !            17:        sqsum = 0;
        !            18: }
        !            19: 
        !            20: void histogram.add(int a)
        !            21: /* add a to one of the bins, adjusting histogram, if necessary */
        !            22: {
        !            23:        register int i, j;
        !            24: 
        !            25:        /* make l <= a < r, */
        !            26:         /*      possibly expanding histogram by doubling binsize and range */
        !            27:        while (a<l) {
        !            28:                l -= r - l;
        !            29:                for (i=nbin-1, j=nbin-2; 0<=j; i--, j-=2) h[i] = h[j] + h[j+1];
        !            30:                while(i >= 0) h[i--] = 0;
        !            31:                binsize += binsize;
        !            32:        }
        !            33:        while (r<=a) {
        !            34:                r += r - l;
        !            35:                for (i=0, j=0; i<nbin/2 ; i++, j+=2) h[i] = h[j] + h[j+1];
        !            36:                while (i < nbin) h[i++] = 0;
        !            37:                binsize += binsize;
        !            38:        }
        !            39:        sum += a;
        !            40:        sqsum += a * a;
        !            41:        h[(a-l)/binsize]++;
        !            42: }
        !            43: 
        !            44: void histogram.print()
        !            45: /*
        !            46:        printout non-empty ranges
        !            47: */
        !            48: {
        !            49:        register int i;
        !            50:        register int x;
        !            51:        int d = binsize;
        !            52:        
        !            53:        for (i=0; i<nbin; i++) {
        !            54:                if (x=h[i]) {
        !            55:                        int ll = l+d*i;
        !            56:                        printf("[%d:%d) : %d\n",ll,ll+d,x);
        !            57:                }
        !            58:        }
        !            59: } 
        !            60: 
        !            61: /*
        !            62: int erand.draw()
        !            63: {
        !            64:        int k;
        !            65:        float a;
        !            66: 
        !            67:        for(k=0;;k++) {
        !            68:                register float u1, u2;
        !            69:                a = u1 = fdraw();
        !            70:                do {
        !            71:                        u2 = fdraw();
        !            72:                        if (u1 < u2) return (int) k+a;
        !            73:                        u1 = fdraw();
        !            74:                } while (u1<u2);
        !            75:        }
        !            76:        
        !            77: }
        !            78: */

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