LRdf.c

Go to the documentation of this file.

/*
 * Copyright    2019    R.K. Owen, Ph.D.
 * License      see lgpl.md (Gnu Lesser General Public License)
 */
#ifdef __cplusplus
extern "C" {
#endif

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

double LRd_RAN(LR_obj *o) {
        if (o->d != LR_double) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->rnd(o);
}

double LRd_PDF(LR_obj *o, double x) {
        if (o->d != LR_double) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->pdfd(o,x);
}

double LRd_CDF(LR_obj *o, double x) {
        if (o->d != LR_double) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->cdfd(o,x);
}

float LRf_RAN(LR_obj *o) {
        if (o->d != LR_float) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->rnf(o);
}

float LRf_PDF(LR_obj *o, float x) {
        if (o->d != LR_float) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->pdff(o,x);
}

float LRf_CDF(LR_obj *o, float x) {
        if (o->d != LR_float) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->cdff(o,x);
}

int LRi_RAN(LR_obj *o) {
        if (o->d != LR_int) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->rni(o);
}

float LRi_PDF(LR_obj *o, int x) {
        if (o->d != LR_int) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->pdfi(o,x);
}

float LRi_CDF(LR_obj *o, int x) {
        if (o->d != LR_int) {
                o->errno = LRerr_BadDataType;
                return NAN;
        }
        return o->cdfi(o,x);
}

int LR_aux_new(LR_obj *o, int n) {
        if (!o->aux)
                return o->errno = LRerr_NoAuxiliaryObject;
        if (!(((LR_pcs *) o->aux)->new))
                return o->errno = LRerr_BadAuxSetup;

        return ((LR_pcs *) o->aux)->new(o,n);
}

int LR_aux_rm(LR_obj *o) {
        if (!o->aux)
                return o->errno = LRerr_NoAuxiliaryObject;
        if (!(((LR_pcs *) o->aux)->rm))
                return o->errno = LRerr_BadAuxSetup;

        return ((LR_pcs *) o->aux)->rm(o);
}

int LR_aux_set(LR_obj *o, double x, double p) {
        if (!o->aux)
                return o->errno = LRerr_NoAuxiliaryObject;
        if (!(((LR_pcs *) o->aux)->set))
                return o->errno = LRerr_BadAuxSetup;

        return ((LR_pcs *) o->aux)->set(o,x,p);
}

int LR_aux_norm(LR_obj *o) {
        int ret;
        if (!o->aux)
                return o->errno = LRerr_NoAuxiliaryObject;
        if (!(((LR_pcs *) o->aux)->normalize))
                return o->errno = LRerr_BadAuxSetup;

        return ((LR_pcs *) o->aux)->normalize(o);
}

#ifdef __cplusplus
}
#endif