diff --git a/GPy/inference/latent_function_inference/svgp.py b/GPy/inference/latent_function_inference/svgp.py
index a7ee9d14..8d99e750 100644
--- a/GPy/inference/latent_function_inference/svgp.py
+++ b/GPy/inference/latent_function_inference/svgp.py
@@ -15,7 +15,6 @@ class SVGP(LatentFunctionInference):
         L = choleskies.flat_to_triang(q_u_chol)
 
 
-        S_ein = np.einsum('ijk,ljk->ilk', L, L) #L.dot(L.T)
         S = np.empty((num_outputs, num_inducing, num_inducing))
         [np.dot(L[:,:,i], L[:,:,i].T, S[i,:,:]) for i in range(num_outputs)]
         S = S.swapaxes(0,2)