diff --git a/src/shogun/preprocessor/HessianLocallyLinearEmbedding.cpp b/src/shogun/preprocessor/HessianLocallyLinearEmbedding.cpp
index 0f47a22f4b4..85e0cdce0f1 100755
--- a/src/shogun/preprocessor/HessianLocallyLinearEmbedding.cpp
+++ b/src/shogun/preprocessor/HessianLocallyLinearEmbedding.cpp
@@ -84,10 +84,7 @@ SGMatrix<float64_t> CHessianLocallyLinearEmbedding::apply_to_feature_matrix(CFea
 	SG_FREE(local_neighbors_idxs);
 
 	// init W (weight) matrix
-	float64_t* W_matrix = SG_MALLOC(float64_t, N*N);
-	for (i=0; i<N; i++)
-		for (j=0; j<N; j++)
-			W_matrix[i*N+j]=0.0;
+	float64_t* W_matrix = SG_CALLOC(float64_t, N*N);
 
 	// init matrices and norm factor to be used
 	float64_t* local_feature_matrix = SG_MALLOC(float64_t, m_k*dim);
diff --git a/src/shogun/preprocessor/LocallyLinearEmbedding.cpp b/src/shogun/preprocessor/LocallyLinearEmbedding.cpp
index 574beca07ef..895228ec752 100755
--- a/src/shogun/preprocessor/LocallyLinearEmbedding.cpp
+++ b/src/shogun/preprocessor/LocallyLinearEmbedding.cpp
@@ -82,10 +82,7 @@ SGMatrix<float64_t> CLocallyLinearEmbedding::apply_to_feature_matrix(CFeatures*
 	SG_FREE(local_neighbors_idxs);
 
 	// init W (weight) matrix
-	float64_t* W_matrix = SG_MALLOC(float64_t, N*N);
-	for (i=0; i<N; i++)
-		for (j=0; j<N; j++)
-			W_matrix[i*N+j]=0.0;
+	float64_t* W_matrix = SG_CALLOC(float64_t, N*N);
 
 	// init matrices and norm factor to be used
 	float64_t* z_matrix = SG_MALLOC(float64_t, m_k*dim);