From: Sarah Westcott Date: Tue, 13 Mar 2012 19:22:00 +0000 (-0400) Subject: bug fix in linear algebra that caused an infinite loop X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=commitdiff_plain;h=3114ad297d40ecd47bd3c1fe7997bd76a81e8c75 bug fix in linear algebra that caused an infinite loop --- diff --git a/linearalgebra.cpp b/linearalgebra.cpp index ab75b7e..effb1ad 100644 --- a/linearalgebra.cpp +++ b/linearalgebra.cpp @@ -222,7 +222,7 @@ double LinearAlgebra::betacf(const double a, const double b, const double x) { if (fabs(d) < FPMIN) d=FPMIN; d=1.0/d; h=d; - for (m1=1;m1<=MAXIT;m++) { + for (m1=1;m1<=MAXIT;m1++) { m2=2*m1; aa=m1*(b-m1)*x/((qam+m2)*(a+m2)); d=1.0+aa*d;