diff --git a/EC.jl b/EC.jl
index de38c14..efd1e12 100644
--- a/EC.jl
+++ b/EC.jl
@@ -77,7 +77,7 @@ end
 function get_reduced_matrices(EC::affine_EC, training_vecs, gram_schmidt_threshold, subsample=1:length(training_vecs); verbose=false)
     vecs = deepcopy(training_vecs[subsample])
 
-    if gram_schmidt_threshold > 0; vecs = gram_schmidt!(sampled_vecs, EC.weights, gram_schmidt_threshold; verbose=verbose); end
+    if gram_schmidt_threshold > 0; vecs = gram_schmidt!(vecs, EC.weights, gram_schmidt_threshold; verbose=verbose); end
 
     EC_basis = hcat(vecs...)
     weights_mat = spdiagm(EC.weights)