Mercurial > repos > miller-lab > genome_diversity
view genome_diversity/src/eval2pct.c @ 0:2c498d40ecde
Uploaded
author | miller-lab |
---|---|
date | Mon, 09 Apr 2012 12:03:06 -0400 |
parents | |
children |
line wrap: on
line source
#include "lib.h" #define MAX_EVAL 1000 float E[MAX_EVAL]; int nE; int main (int argc, char **argv) { FILE *fp; char buf[500]; int i; float tot; fp = (argc== 1 ? stdin : ckopen(argv[1], "r")); while (fgets(buf, 500, fp)) { if (nE >= MAX_EVAL) fatal("Too many eigenvalues"); E[nE++] = atof(buf); } for (tot = 0.0, i = 0; i < nE; ++i) tot += E[i]; printf("Percentage explained by eigenvectors:\n"); for (i = 0 ; i < nE && E[i] > 0.0; ++i) printf("%d: %1.1f%%\n", i+1, 100.0*(float)E[i]/tot); return 0; }