Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Already on GitHub? Sign in to your account

Avoid divide by zero in average length, std #64

Open
wants to merge 1 commit into
from
Jump to file or symbol
Failed to load files and symbols.
+10 −2
Split
View
@@ -388,9 +388,17 @@ int main (int argc, char **argv) {
double dev = sqrt((((double)joincnt)*tlensq-pow((double)tlen,2)) / ((double)joincnt*((double)joincnt-1)) );
+ double avg = (double) tlen / (double) joincnt;
+ if(joincnt == 0) {
+ dev = 0;
+ avg = 0;
+ }
+ if(joincnt == 1) {
+ dev = 0;
+ }
printf("Total reads: %d\n", nrec);
printf("Total joined: %d\n", joincnt);
- printf("Average join len: %.2f\n", (double) tlen / (double) joincnt);
+ printf("Average join len: %.2f\n", avg);
printf("Stdev join len: %.2f\n", dev);
printf("Version: %s\n", VERSION);
@@ -1,5 +1,5 @@
Total reads: 1
Total joined: 1
Average join len: 50.00
-Stdev join len: -nan
+Stdev join len: 0.00
Version: 1.01.759