jacobian.hpp

Go to the documentation of this file.

#ifndef __STAN__AGRAD__REV__JACOBIAN_HPP__
#define __STAN__AGRAD__REV__JACOBIAN_HPP__

#include <vector>
#include <stan/agrad/rev/var.hpp>

namespace stan {
  namespace agrad {

    inline void jacobian(std::vector<var>& dependents,
                         std::vector<var>& independents,
                         std::vector<std::vector<double> >& jacobian) {
      jacobian.resize(dependents.size());
      for (size_t i = 0; i < dependents.size(); ++i) {
        jacobian[i].resize(independents.size());
        if (i > 0) 
          set_zero_all_adjoints();
        jacobian.push_back(std::vector<double>(0));
        grad(dependents[i].vi_);
        for (size_t j = 0; j < independents.size(); ++j)
          jacobian[i][j] = independents[j].adj();
      }
    }

  }
}
#endif