@@ -912,47 +912,48 @@ TEST_F(AgradRev, matrix_compile_time_conversions) {
912912}
913913
914914TEST_F (AgradRev, assign_nan) {
915- using stan::math::var_value;
916- using var_vector = var_value<Eigen::Matrix<double ,-1 ,1 >>;
917- using stan::math::var;
918- Eigen::VectorXd x_val (10 );
919- for (int i = 0 ; i < 10 ; ++i) {
920- x_val (i) = i + 0.1 ;
921- }
922- var_vector x (x_val);
923- var_vector y = var_vector (Eigen::Matrix<double ,-1 ,1 >::Constant (10 , std::numeric_limits<double >::quiet_NaN ()));
924- y = stan::math::head (x, 10 );
925- var sigma = 1.0 ;
926- var lp = stan::math::normal_lpdf<false >(y, 0 , sigma);
927- lp.grad ();
928- Eigen::VectorXd x_ans_adj (10 );
929- for (int i = 0 ; i < 10 ; ++i) {
930- x_ans_adj (i) = -(i + 0.1 );
931- }
932- EXPECT_MATRIX_EQ (x.adj (), x_ans_adj);
933- Eigen::VectorXd y_ans_adj = Eigen::VectorXd::Zero (10 );
934- EXPECT_MATRIX_EQ (y_ans_adj, y.adj ());
915+ using stan::math::var_value;
916+ using var_vector = var_value<Eigen::Matrix<double , -1 , 1 >>;
917+ using stan::math::var;
918+ Eigen::VectorXd x_val (10 );
919+ for (int i = 0 ; i < 10 ; ++i) {
920+ x_val (i) = i + 0.1 ;
921+ }
922+ var_vector x (x_val);
923+ var_vector y = var_vector (Eigen::Matrix<double , -1 , 1 >::Constant (
924+ 10 , std::numeric_limits<double >::quiet_NaN ()));
925+ y = stan::math::head (x, 10 );
926+ var sigma = 1.0 ;
927+ var lp = stan::math::normal_lpdf<false >(y, 0 , sigma);
928+ lp.grad ();
929+ Eigen::VectorXd x_ans_adj (10 );
930+ for (int i = 0 ; i < 10 ; ++i) {
931+ x_ans_adj (i) = -(i + 0.1 );
932+ }
933+ EXPECT_MATRIX_EQ (x.adj (), x_ans_adj);
934+ Eigen::VectorXd y_ans_adj = Eigen::VectorXd::Zero (10 );
935+ EXPECT_MATRIX_EQ (y_ans_adj, y.adj ());
935936}
936937
937938TEST_F (AgradRev, assign_nullptr_vari) {
938- using stan::math::var_value;
939- using var_vector = var_value<Eigen::Matrix<double ,-1 ,1 >>;
940- using stan::math::var;
941- Eigen::VectorXd x_val (10 );
942- for (int i = 0 ; i < 10 ; ++i) {
943- x_val (i) = i + 0.1 ;
944- }
945- var_vector x (x_val);
946- var_vector y;
947- y = stan::math::head (x, 10 );
948- var sigma = 1.0 ;
949- var lp = stan::math::normal_lpdf<false >(y, 0 , sigma);
950- lp.grad ();
951- Eigen::VectorXd x_ans_adj (10 );
952- for (int i = 0 ; i < 10 ; ++i) {
953- x_ans_adj (i) = -(i + 0.1 );
954- }
955- EXPECT_MATRIX_EQ (x.adj (), x_ans_adj);
956- Eigen::VectorXd y_ans_adj = Eigen::VectorXd::Zero (10 );
957- EXPECT_MATRIX_EQ (y_ans_adj, y.adj ());
939+ using stan::math::var_value;
940+ using var_vector = var_value<Eigen::Matrix<double , -1 , 1 >>;
941+ using stan::math::var;
942+ Eigen::VectorXd x_val (10 );
943+ for (int i = 0 ; i < 10 ; ++i) {
944+ x_val (i) = i + 0.1 ;
945+ }
946+ var_vector x (x_val);
947+ var_vector y;
948+ y = stan::math::head (x, 10 );
949+ var sigma = 1.0 ;
950+ var lp = stan::math::normal_lpdf<false >(y, 0 , sigma);
951+ lp.grad ();
952+ Eigen::VectorXd x_ans_adj (10 );
953+ for (int i = 0 ; i < 10 ; ++i) {
954+ x_ans_adj (i) = -(i + 0.1 );
955+ }
956+ EXPECT_MATRIX_EQ (x.adj (), x_ans_adj);
957+ Eigen::VectorXd y_ans_adj = Eigen::VectorXd::Zero (10 );
958+ EXPECT_MATRIX_EQ (y_ans_adj, y.adj ());
958959}
0 commit comments