use 0 for fill function

include/micm/solver/rosenbrock.inl

@@ -56,12 +56,12 @@ namespace micm
       H = std::min(H, std::abs(time_step - present_time));
 
       // compute the initial forcing at the beginning of the current time
-      initial_forcing.Fill(0.0);
+      initial_forcing.Fill(0);
       rates_.AddForcingTerms(state.rate_constants_, Y, initial_forcing);
       stats.function_calls_ += 1;
 
       // compute the negative jacobian at the beginning of the current time
-      state.jacobian_.Fill(0.0);
+      state.jacobian_.Fill(0);
       rates_.SubtractJacobianTerms(state.rate_constants_, Y, state.jacobian_);
       stats.jacobian_updates_ += 1;
 
@@ -95,7 +95,7 @@ namespace micm
               {
                 Ynew.Axpy(parameters_.a_[stage_combinations + j], K[j]);
               }
-              K[stage].Fill(0.0);
+              K[stage].Fill(0);
               rates_.AddForcingTerms(state.rate_constants_, Ynew, K[stage]);
               stats.function_calls_ += 1;
             }
@@ -256,7 +256,7 @@ namespace micm
         break;
       H /= 2;
       // Reconstruct the Jacobian matrix if a substepping is performed here
-      state.jacobian_.Fill(0.0);
+      state.jacobian_.Fill(0);
       rates_.SubtractJacobianTerms(state.rate_constants_, number_densities, state.jacobian_);
       stats.jacobian_updates_ += 1;
     }