Adjoint ode_bdf passing ode_bdf tests. integrate_ode_bdf not working yet. (design-doc #19)

Author: bbbales2

stan/math/prim/functor/ode_rk45.hpp

@@ -163,8 +163,6 @@ ode_rk45_tol(const F& f,
   bool observer_initial_recorded = false;
   size_t time_index = 0;
 
-  max_step_checker step_checker(max_num_steps);
-
   // avoid recording of the initial state which is included by the
   // conventions of odeint in the output
   auto filtered_observer
@@ -175,7 +173,6 @@ ode_rk45_tol(const F& f,
     }
     y.emplace_back(ode_store_sensitivities(f, coupled_state, y0, t0,
                                            ts[time_index], msgs, args...));
-    step_checker.reset();
     time_index++;
   };
 
@@ -189,7 +186,8 @@ ode_rk45_tol(const F& f,
           runge_kutta_dopri5<Eigen::VectorXd, double, Eigen::VectorXd, double,
                              vector_space_algebra>()),
       std::ref(coupled_system), initial_coupled_state, std::begin(ts_vec),
-      std::end(ts_vec), step_size, filtered_observer, step_checker);
+      std::end(ts_vec), step_size, filtered_observer,
+      max_step_checker(max_num_steps));
 
   return y;
 }

stan/math/rev/core/grad.hpp

@@ -7,6 +7,7 @@
 #include <stan/math/rev/core/nested_size.hpp>
 #include <stan/math/rev/core/vari.hpp>
 #include <vector>
+#include <iostream>
 
 namespace stan {
 namespace math {
@@ -27,8 +28,9 @@ namespace math {
  * @param vi Variable implementation for root of partial
  * derivative propagation.
  */
-static void grad(vari* vi) {
-  vi->init_dependent();
+static void grad(vari* vi = NULL) {
+  if(vi != NULL)
+    vi->init_dependent();
   std::vector<vari*>& var_stack = ChainableStack::instance_->var_stack_;
   size_t end = var_stack.size();
   size_t beginning = empty_nested() ? 0 : end - nested_size();