Lecture du code C écrit pour les paquets R. Où trouver des définitions pour RANDIN, UNIF, EPS, etc?

Je suis intéressé par la compréhension précise du code C pour les paquets R. Mais j'ai rencontré un code tel que RANDIN , RANDOUT , UNIF , EPS etc (sont-ils des macros?) pour lequel je ne sais pas où trouver les implémentations/définitions.

Où puis-je trouver le code (et les explications le cas échéant) pour ces «expressions» capitalisées?

Un exemple est la fonction VR_onlineSOM dans le paquet R de 'classe'. Cette fonction a le code suivant écrit en C:

void
VR_onlineSOM(double *data, double *codes, double *nhbrdist,
         double *alpha, double *radii,
         Sint *pn, Sint *pp, Sint *pncodes, Sint *rlen)
{
    int n = *pn, p = *pp, ncodes = *pncodes;
    int i, j, k, nearest = 0 /* -Wall */, nind;
    double dm, dist, tmp;
    unsigned int cd; /* avoid spurious warning from gcc pre-4.3.0 */

    RANDIN;
    for (k = 0; k < *rlen; k++) {
    /* pick a random data point */
    i = (int)(n * UNIF);
    /* find the nearest code 'near' */
    nind = 0; dm = DOUBLE_XMAX;
    for (cd = 0; cd < ncodes; cd++) {
        dist = 0.0;
        for (j = 0; j < p; j++) {
        tmp = data[i + j*n] - codes[cd + j*ncodes];
        dist += tmp * tmp;
        }
        if (dist <= dm * (1 + EPS)) {
        if (dist < dm * (1 - EPS)) {
            nind = 0;
            nearest = cd;
        } else {
            if(++nind * UNIF < 1.0) nearest = cd;
        }
        dm = dist;
        }
        /* update all codes within radii[k] of 'nearest' */
        for (cd = 0; cd < ncodes; cd++) {
        if(nhbrdist[cd + ncodes*nearest] > radii[k]) continue;
        for(j = 0; j < p; j++)
            codes[cd + j*ncodes] += alpha[k] *
            (data[i + j*n] - codes[cd + j*ncodes]);
        }
    }
    }
    RANDOUT;
}
1

1 Réponses

Oui, ce sont des macros et des constantes définies au début du class.c .

#define EPS 1e-4            /* relative test of equality of distances */
#define RANDIN  GetRNGstate()
#define RANDOUT PutRNGstate()
#define UNIF unif_rand()
3
ajouté