Return to sim.c CVS log

Up to [Research Unix] / researchv10no / cmd / cfront / libC / otask

Annotation of researchv10no/cmd/cfront/libC/otask/sim.c, revision 1.1.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) ((object*)0)->task_error(E_HISTO);
                      7:        if (nb%2) nb++;
                      8:        while ((rr-ll)%nb) rr++;
                      9:        binsize = (rr-ll)/nb;
                     10:        h = new int[nb];
                     11:        for (i=0; i<nb; i++) h[i] = 0;
                     12:        l = ll;
                     13:        r = rr;
                     14:        nbin = nb;
                     15:        sum = 0;
                     16:        sqsum = 0;
                     17: }
                     18: 
                     19: void
                     20: 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
                     45: histogram::print()
                     46: /*
                     47:        printout non-empty ranges
                     48: */
                     49: {
                     50:        register int i;
                     51:        register int x;
                     52:        int d = binsize;
                     53:        
                     54:        for (i=0; i<nbin; i++) {
                     55:                if (x=h[i]) {
                     56:                        int ll = l+d*i;
                     57:                        printf("[%d:%d) : %d\n",ll,ll+d,x);
                     58:                }
                     59:        }
                     60: } 
                     61: 
                     62: /*
                     63: int
                     64: erand::draw()
                     65: {
                     66:        int k;
                     67:        float a;
                     68: 
                     69:        for(k=0;;k++) {
                     70:                register float u1, u2;
                     71:                a = u1 = fdraw();
                     72:                do {
                     73:                        u2 = fdraw();
                     74:                        if (u1 < u2) return (int) k+a;
                     75:                        u1 = fdraw();
                     76:                } while (u1<u2);
                     77:        }
                     78:        
                     79: }
                     80: */

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