operator_division.hpp

Go to the documentation of this file.

#ifndef STAN__AGRAD__FWD__OPERATORS__OPERATOR_DIVISION_HPP
#define STAN__AGRAD__FWD__OPERATORS__OPERATOR_DIVISION_HPP

#include <stan/agrad/fwd/fvar.hpp>
#include <stan/meta/traits.hpp>

namespace stan {

  namespace agrad {

    template <typename T>
    inline
    fvar<T>
    operator/(const fvar<T>& x1, const fvar<T>& x2) {
      return fvar<T>(x1.val_ / x2.val_, 
                     (x1.d_ * x2.val_ - x1.val_ * x2.d_) / (x2.val_ * x2.val_));
    }

    template <typename T>
    inline
    fvar<T>
    operator/(const fvar<T>& x1, const double x2) {
      return fvar<T>(x1.val_ / x2,
                     x1.d_ / x2);
    }

    template <typename T>
    inline
    fvar<T>
    operator/(const double x1, const fvar<T>& x2) {
      return fvar<T>(x1 / x2.val_, 
                     - x1 * x2.d_ / (x2.val_ * x2.val_));
                     
    }
  }
}
#endif