3 files changed, 6 insertions, 3 deletions
diff --git a/src/c/statisticsFunctions/correl/dcorrelfres.c b/src/c/statisticsFunctions/correl/dcorrelfres.c
index 7d4c718a..34389353 100644
--- a/src/c/statisticsFunctions/correl/dcorrelfres.c
+++ b/src/c/statisticsFunctions/correl/dcorrelfres.c
@@ -11,10 +11,11 @@
 #include "correl.h"
 #include "matrix.h"
 #include "mean.h"
+#include "matrixMultiplication.h"
 
 double dcorrelfres (double* inp1, int sr1, int sc1, double* inp2, int sr2, int sc2, double* inp3, int sr3, int sc3)
 {
-    double p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], sum = 0, prod[sr3];
+    double p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], prod[sr3];
 
     double fr[sr3*sc3];
 
diff --git a/src/c/statisticsFunctions/correl/scorrelfres.c b/src/c/statisticsFunctions/correl/scorrelfres.c
index 2d056a5b..066063cd 100644
--- a/src/c/statisticsFunctions/correl/scorrelfres.c
+++ b/src/c/statisticsFunctions/correl/scorrelfres.c
@@ -11,10 +11,11 @@
 #include "correl.h"
 #include "matrix.h"
 #include "mean.h"
+#include "matrixMultiplication.h"
 
 float scorrelfres (float* inp1, int sr1, int sc1, float* inp2, int sr2, int sc2, float* inp3, int sr3, int sc3)
 {
-    float p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], sum = 0, prod[sr3];
+    float p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], prod[sr3];
 
     float fr[sr3*sc3];
 
diff --git a/src/c/statisticsFunctions/correl/zcorrelfres.c b/src/c/statisticsFunctions/correl/zcorrelfres.c
index b026c58a..a183ec77 100644
--- a/src/c/statisticsFunctions/correl/zcorrelfres.c
+++ b/src/c/statisticsFunctions/correl/zcorrelfres.c
@@ -13,10 +13,11 @@
 #include "matrix.h"
 #include "mean.h"
 #include "doubleComplex.h"
+#include "matrixMultiplication.h"
 
 doubleComplex zcorrelfres (doubleComplex* inp1, int sr1, int sc1, doubleComplex* inp2, int sr2, int sc2, doubleComplex* inp3, int sr3, int sc3)
 {
-    doubleComplex p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], sum = 0, prod[sr3];
+    doubleComplex p1[sr3], p2[sc3], m1 = 0, m2 = 0,  s1 = 0, s2 = 0, s3 = 0, diff1[sr1*sc1], diff2[sr2*sc2], prod[sr3];
 
     doubleComplex fr[sr3*sc3];