 #endif
 
 	double beta = 0;
-	double A[NN * NN] = {ZERO};
-	double Hm[STRIDE * STRIDE] = {ZERO};
-	double Vm[NN * STRIDE];
-	double phiHm[STRIDE * STRIDE];
-	double savedActions[NN * 5] = {ZERO};
 	while ((t < t_end) && (t + h > t)) {
 		//initial krylov subspace sizes
 		int m, m1, m2;

    		gy[i] = fy[i] - gy[i];
    	}

		double err = ZERO;
		do
		{
			//do arnoldi

		dydt (t, pr, y, fy);

		// Jacobian matrix
		double A[NN * NN] = {ZERO};
		eval_jacob (t, pr, y, A);
    	double gy[NN];
    	//gy = fy - A * y
