| 1 |
greg |
2.1 |
/* RCSid $Id$ */
|
| 2 |
|
|
/*
|
| 3 |
|
|
** Author: Christian Reetz ([email protected])
|
| 4 |
|
|
*/
|
| 5 |
|
|
#ifndef __MUC_RANDVAR_H
|
| 6 |
|
|
#define __MUC_RANDVAR_H
|
| 7 |
|
|
|
| 8 |
|
|
/** \file muc_randvar.h
|
| 9 |
|
|
* \author Ch. Reetz ([email protected])
|
| 10 |
|
|
* \brief representation of a random variables of arbitrary dimension
|
| 11 |
|
|
*/
|
| 12 |
|
|
|
| 13 |
|
|
|
| 14 |
|
|
#ifdef __cplusplus
|
| 15 |
|
|
extern "C" {
|
| 16 |
|
|
#endif
|
| 17 |
|
|
|
| 18 |
|
|
#include "g3flist.h"
|
| 19 |
|
|
|
| 20 |
|
|
/** random variable structure*/
|
| 21 |
|
|
struct muc_rvar {
|
| 22 |
|
|
int n; /* number of samples */
|
| 23 |
|
|
double w; /* sum of weights*/
|
| 24 |
|
|
double* sum; /* sample values sum*/
|
| 25 |
|
|
double* sum_sqr; /* squared sample values sum*/
|
| 26 |
|
|
double* min; /* minimum component value*/
|
| 27 |
|
|
double* max; /* maximum component value*/
|
| 28 |
|
|
g3FList* samples; /* measured samples (if saved)*/
|
| 29 |
|
|
int save_samp; /* if true (default) save samples*/
|
| 30 |
|
|
};
|
| 31 |
|
|
|
| 32 |
|
|
/* returns dimension of random variable*/
|
| 33 |
|
|
#define muc_rvar_get_dim(rv) (g3fl_get_comp_size(rv->samples))
|
| 34 |
|
|
/* returns the number of samples*/
|
| 35 |
|
|
#define muc_rvar_get_sample_size(rv) (rv->n)
|
| 36 |
|
|
|
| 37 |
|
|
/* returns a reference to a sample list element*/
|
| 38 |
|
|
/* Take care to have save_samp flag set*/
|
| 39 |
|
|
#define muc_rvar_get_sample(rv,id) (g3fl_get(rv->samples, id))
|
| 40 |
|
|
|
| 41 |
|
|
|
| 42 |
|
|
int muc_rvar_init(struct muc_rvar* rv);
|
| 43 |
|
|
struct muc_rvar* muc_rvar_create();
|
| 44 |
|
|
void muc_rvar_free(struct muc_rvar* rv);
|
| 45 |
|
|
|
| 46 |
|
|
/* set dimension of random variable (default 1). */
|
| 47 |
|
|
/* random variable needs to be empty (\see muc_rvar_clear)*/
|
| 48 |
|
|
int muc_rvar_set_dim(struct muc_rvar* rv,int dim);
|
| 49 |
|
|
|
| 50 |
|
|
/* set storing samples on (default) or off*/
|
| 51 |
|
|
/* random varianle must be cleared to succeed*/
|
| 52 |
|
|
int muc_rvar_store_sample(struct muc_rvar* rv, int on);
|
| 53 |
|
|
|
| 54 |
|
|
/* clear random variables samples. Don't change dimension.*/
|
| 55 |
|
|
void muc_rvar_clear(struct muc_rvar* rv);
|
| 56 |
|
|
|
| 57 |
|
|
/* Don't store measured sample values, just summed values.*/
|
| 58 |
|
|
/* (Deletes samples if there are already stored some)*/
|
| 59 |
|
|
void muc_rvar_no_samples(struct muc_rvar* rv);
|
| 60 |
|
|
|
| 61 |
|
|
/* adds measured sample to sums (and stores them if \ref save_samp is true)*/
|
| 62 |
|
|
/* The weight of the sample is assumed to be 1.0*/
|
| 63 |
|
|
int muc_rvar_add_sample(struct muc_rvar* rv,double* sample);
|
| 64 |
|
|
|
| 65 |
|
|
/* adds \param w weighted sample to sums */
|
| 66 |
|
|
/* (and stores them if \ref save_samp is true)*/
|
| 67 |
|
|
int muc_rvar_add_weighted_sample(struct muc_rvar* rv,double w,double* samp);
|
| 68 |
|
|
|
| 69 |
|
|
/* returns current expected value in \param ex*/
|
| 70 |
|
|
int muc_rvar_get_ex(struct muc_rvar* rv,double* ex);
|
| 71 |
|
|
|
| 72 |
|
|
/* returns current variance in \param vx*/
|
| 73 |
|
|
int muc_rvar_get_vx(struct muc_rvar* rv,double* vx);
|
| 74 |
|
|
|
| 75 |
|
|
/* returns current bounds*/
|
| 76 |
|
|
/* \param bounds: array with len 2*component_size */
|
| 77 |
|
|
/* \return {min1, max1, min2, max2....}*/
|
| 78 |
|
|
int muc_rvar_get_bounding_box(struct muc_rvar* rv, double* bounds);
|
| 79 |
|
|
|
| 80 |
|
|
/* return sample median (using averaging for even sample numbers)*/
|
| 81 |
|
|
/* \return false if samples are not stored*/
|
| 82 |
|
|
int muc_rvar_get_median(struct muc_rvar* rv, double* median);
|
| 83 |
|
|
|
| 84 |
|
|
int muc_rvar_get_percentile(struct muc_rvar* rv, double* median, double percentile);
|
| 85 |
|
|
|
| 86 |
|
|
#ifdef __cplusplus
|
| 87 |
|
|
}
|
| 88 |
|
|
#endif
|
| 89 |
|
|
#endif
|
