diff --git a/matlab/default_option_values.m b/matlab/default_option_values.m
index 39b77382301a0265ba1a149e6ba93230f711a863..712fd6ef185b1cd349de351c45d5b3f20e29a0c9 100644
--- a/matlab/default_option_values.m
+++ b/matlab/default_option_values.m
@@ -758,9 +758,6 @@ options_.endogenous_prior_restrictions.moment={};
 % OSR Optimal Simple Rules
 options_.osr.opt_algo=4;
 
-% use GPU
-options_.gpu = false;
-
 %Geweke convergence diagnostics
 options_.convergence.geweke.taper_steps=[4 8 15];
 options_.convergence.geweke.geweke_interval=[0.2 0.5];
diff --git a/matlab/dyn_first_order_solver.m b/matlab/dyn_first_order_solver.m
index 23fb281a701494ce33f480fbf72d640e86cf17ce..41d49cbdd543f4ea0ba539dedd37e2915d127229 100644
--- a/matlab/dyn_first_order_solver.m
+++ b/matlab/dyn_first_order_solver.m
@@ -159,11 +159,6 @@ if task ~= 1 && (DynareOptions.dr_cycle_reduction || DynareOptions.dr_logarithmi
                    'coefficient matrix for current variables isn''t invertible'])
         end
     end
-    if DynareOptions.gpu
-        gpuArray(A1);
-        gpuArray(B1);
-        gpuArray(C1);
-    end
     A1 = [aa(row_indx,index_m ) zeros(ndynamic,nfwrd)];
     B1 = [aa(row_indx,index_0m) aa(row_indx,index_0p) ];
     C1 = [zeros(ndynamic,npred) aa(row_indx,index_p)];