/* nag_deviates_gamma_vector (g01tfc) Example Program. * * Copyright 2011, Numerical Algorithms Group. * * Mark 23, 2011. */ #include #include #include #include #define ENUMLEN 30 int main(void) { /* Integer scalar and array declarations */ Integer ltail, lp, la, lb, i, lout; Integer *ivalid = 0; Integer exit_status = 0; /* NAG structures */ NagError fail; Nag_TailProbability *tail = 0; /* Double scalar and array declarations */ double tol; double *p = 0, *a = 0, *b = 0, *g = 0; /* Character scalar and array declarations */ char ctail[ENUMLEN]; /* Initialise the error structure to print out any error messages */ INIT_FAIL(fail); printf("nag_deviates_gamma_vector (g01tfc) Example Program Results\n\n"); /* Skip heading in data file*/ scanf("%*[^\n] "); /* Read in the tolerance */ scanf("%lf%*[^\n] ", &tol); /* Read in the input vectors */ scanf("%ld%*[^\n] ", <ail); if (!(tail = NAG_ALLOC(ltail, Nag_TailProbability))) { printf("Allocation failure\n"); exit_status = -1; goto END; } for (i = 0; i < ltail; i++) { scanf("%s", ctail); tail[i] = (Nag_TailProbability) nag_enum_name_to_value(ctail); } scanf("%*[^\n] "); scanf("%ld%*[^\n] ", &lp); if (!(p = NAG_ALLOC(lp, double))) { printf("Allocation failure\n"); exit_status = -1; goto END; } for (i = 0; i < lp; i++) scanf("%lf", &p[i]); scanf("%*[^\n] "); scanf("%ld%*[^\n] ", &la); if (!(a = NAG_ALLOC(la, double))) { printf("Allocation failure\n"); exit_status = -1; goto END; } for (i = 0; i < la; i++) scanf("%lf", &a[i]); scanf("%*[^\n] "); scanf("%ld%*[^\n] ", &lb); if (!(b = NAG_ALLOC(lb, double))) { printf("Allocation failure\n"); exit_status = -1; goto END; } for (i = 0; i < lb; i++) scanf("%lf", &b[i]); scanf("%*[^\n] "); /* Allocate memory for output */ lout = MAX(ltail,MAX(lp,MAX(la,lb))); if (!(g = NAG_ALLOC(lout, double)) || !(ivalid = NAG_ALLOC(lout, Integer))) { printf("Allocation failure\n"); exit_status = -1; goto END; } /* Calculate probability */ nag_deviates_gamma_vector(ltail, tail, lp, p, la, a, lb, b, tol, g, ivalid, &fail); if (fail.code != NE_NOERROR) { printf("Error from nag_deviates_gamma_vector (g01tfc).\n%s\n", fail.message); exit_status = 1; if (fail.code != NW_IVALID) goto END; } /* Display title */ printf(" tail p a b g ivalid\n"); printf(" ---------------------------------------------------------------\n"); /* Display results */ for (i = 0; i < lout; i++) printf(" %15s %6.3f %6.2f %6.2f %7.3f %3ld\n", nag_enum_value_to_name(tail[i%ltail]),p[i%lp], a[i%la], b[i%lb], g[i], ivalid[i]); END: if(tail) NAG_FREE(tail); if(p) NAG_FREE(p); if(a) NAG_FREE(a); if(b) NAG_FREE(b); if(g) NAG_FREE(g); if(ivalid) NAG_FREE(ivalid); return(exit_status); }