     		gy[i] = fy[i] - gy[i];
     	}
 
-		double err = ZERO;
+		double Hm[STRIDE * STRIDE] = {ZERO};
+		double Vm[NN * STRIDE] = {ZERO};
+		double phiHm[STRIDE * STRIDE] = {ZERO};
+ 		double err = ZERO;
+		double savedActions[NN * 5];
 		do
 		{
 			//do arnoldi

#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;

		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
