eppp.c

#include <stdio.h>
#include "libran.h"

int main() {
        int     r = 0;
        LR_obj *o = LR_new(poisson, LR_int);
        int     nint = 20; 
        double  begin = 0.5, incr;
        double  delh = .0001;
        float   x, pdf, cdf, xcdf[nint];
        long    tot = 50*10007;
        LR_bin *b = LR_bin_new(nint+1);

        LR_set_all(o, "p", 15.);
        r = LR_check(o);
        printf("check r = %d\n", r);

        printf("  k     PDF     CDF\n");
if (0) {
        for (int i = 0; i <= nint; i++) {
                cdf = LRi_CDF(o,i);
                pdf = LRi_PDF(o,i);
                printf("% 3d %7.4f %7.4f\n", i,pdf,cdf);
        }
} else {
        for (int i = 0; i < nint; i++) {
                x   = begin + i;
                LR_bin_set(b,x);
                xcdf[i] = tot * LRi_PDF(o, i);
        }
        xcdf[nint] = tot * (1.0 - LRi_CDF(o, nint-1));

        for (int i = 0; i < tot; i++) {
                r = LR_bin_add(b,x = (double) LRi_RAN(o));
                if (r)
                        printf("%.3f %d\n", x, r);
        }

        printf("   i       bin     cdf     bin/cdf    diff\n");
        for (int i = 0; i <= nint; i++) {
                printf("% 3d  %8ld % 9.2f % 8.2f % 8.2f\n",
                        i, b->bins[i], xcdf[i],
                        b->bins[i]/xcdf[i], b->bins[i] - xcdf[i]);
        }
}

        LR_rm(&o);
        printf("%ld\n", sizeof(o));
}