Merge pull request #770 from stan-dev/bugfix/issue-769-MatrixExpReference

Add passing by reference for matrix_exp function.

stan/math/prim/mat/fun/matrix_exp.hpp

@@ -18,7 +18,7 @@ namespace math {
  */
 template <typename T>
 inline Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> matrix_exp(
-    const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> A) {
+    const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic>& A) {
   check_nonzero_size("matrix_exp", "input matrix", A);
   check_square("matrix_exp", "input matrix", A);