Code clean up

Author: jchristopherson

src/allan.f90

@@ -1,5 +1,6 @@
 submodule (fstats) allan
     use fstats_errors
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 ! REF: Yadav, Shrikanth & Shastri, Saurav & Chakravarthi, Ghanashyam & Kumar, 

src/distributions_chisquared.f90

@@ -1,4 +1,5 @@
 submodule (fstats) distributions_chisquared
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 pure module elemental function cs_pdf(this, x) result(rst)

src/distributions_f.f90

@@ -1,5 +1,6 @@
 submodule (fstats) distributions_f
     use ieee_arithmetic
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 pure module elemental function fd_pdf(this, x) result(rst)

src/distributions_normal.f90

@@ -1,4 +1,5 @@
 submodule (fstats) distributions_normal
+    implicit none
     real(real64), parameter :: pi = 2.0d0 * acos(0.0d0)
 contains
 ! ------------------------------------------------------------------------------

src/distributions_t.f90

@@ -1,5 +1,6 @@
 submodule (fstats) distributions_t
     use ieee_arithmetic
+    implicit none
     real(real64), parameter :: pi = 2.0d0 * acos(0.0d0)
 contains
 ! ------------------------------------------------------------------------------

src/experimental_design_implementation.f90

@@ -1,5 +1,6 @@
 submodule (fstats) experimental_design_implementation
     use fstats_errors
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 module subroutine get_full_factorial_matrix_size_int32(vars, m, n, err)

src/levenberg_marquardt.f90

@@ -3,6 +3,7 @@
 ! 1. https://people.duke.edu/~hpgavin/ExperimentalSystems/lm.pdf
     use linalg
     use fstats_errors
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
     module subroutine regression_jacobian_1(fun, xdata, params, &
@@ -803,7 +804,7 @@ subroutine lm_solve(fun, xdata, ydata, p, weights, maxP, minP, controls, &
         allocate(character(len = 512) :: errmsg)
         write(errmsg, 100) "Memory allocation error code ", flag, "."
         call err%report_error("lm_solve", &
-            trim(errmsg), ML_OUT_OF_MEMORY_ERROR)
+            trim(errmsg), FS_MEMORY_ERROR)
         return
 
         ! Formatting

src/regression_implementation.f90

@@ -1,6 +1,7 @@
 submodule (fstats) regression_implementation
     use linalg
     use fstats_errors
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 module subroutine coefficient_matrix_real64(order, intercept, x, c, err)

src/special_functions_beta.f90

@@ -1,4 +1,5 @@
 submodule (fstats) special_functions_beta
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 pure elemental module function beta_real64(a, b) result(rst)

src/special_functions_digamma.f90

@@ -1,5 +1,6 @@
 submodule (fstats) special_functions_digamma
     use ieee_arithmetic
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 pure elemental module function digamma_real64(x) result(rst)

src/special_functions_gamma.f90

@@ -1,5 +1,6 @@
 submodule (fstats) special_functions_gamma
     use ieee_arithmetic
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 ! REF: https://people.math.sc.edu/Burkardt/f_src/special_functions/special_functions.f90

src/statistics_implementation.f90

@@ -2,6 +2,7 @@
     use linalg, only : sort
     use ieee_arithmetic
     use fstats_errors
+    implicit none
 contains
 ! ------------------------------------------------------------------------------
 pure module function mean_real64(x) result(rst)
@@ -346,7 +347,7 @@ pure module function quantile_real64(x, q) result(rst)
     real(real64), parameter :: one = 1.0d0
 
     ! Local Variables
-    real(real64) :: a, b, c
+    real(real64) :: a, b, c, tol
     integer(int32) :: n, ib
 
     ! Initialization
@@ -376,7 +377,7 @@ pure module function quantile_real32(x, q) result(rst)
     real(real32), parameter :: one = 1.0
 
     ! Local Variables
-    real(real32) :: a, b, c
+    real(real32) :: a, b, c, tol
     integer(int32) :: n, ib
 
     ! Initialization
@@ -614,7 +615,7 @@ module subroutine t_test_paired_real64(x1, x2, stat, p, dof, err)
     class(errors), pointer :: errmgr
     type(errors), target :: deferr
     real(real64) :: v1, v2, m1, m2, sd, cov, a, b, x
-    integer(int32) :: n1, n2, n
+    integer(int32) :: i, n1, n2, n
 
     ! Initialization
     if (present(err)) then
@@ -671,7 +672,7 @@ module subroutine t_test_paired_real32(x1, x2, stat, p, dof, err)
     class(errors), pointer :: errmgr
     type(errors), target :: deferr
     real(real32) :: v1, v2, m1, m2, sd, cov, a, b, x
-    integer(int32) :: n1, n2, n
+    integer(int32) :: i, n1, n2, n
 
     ! Initialization
     if (present(err)) then
@@ -881,7 +882,7 @@ module function anova_2_factor(x) result(rst)
 
     ! Local Variables
     integer(int32) :: i, j, jj, k, r, n
-    real(real64) :: factorMean
+    real(real64) :: factorMean, sum_all
     real(real64), allocatable :: xpack(:)
 
     ! Initialization
@@ -994,6 +995,7 @@ module function anova_model_fit(nmodelparams, ymeas, ymod, err) result(rst)
     ! Local Variables
     integer(int32) :: n, flag
     real(real64), allocatable :: ypack(:)
+    real(real64) :: sum_all
     class(errors), pointer :: errmgr
     type(errors), target :: deferr