Ejemplo n.º 1
0
TEST(AgradRev,inv_square) {
  AVAR a = 7.0;
  AVEC x = createAVEC(a);
  AVAR f = inv_square(a);
  EXPECT_FLOAT_EQ(1 / 49.0, f.val());

  VEC grad_f;
  f.grad(x,grad_f);
  EXPECT_EQ(1U,grad_f.size());
  EXPECT_FLOAT_EQ(-2.0 / 343.0, grad_f[0]);

  a = 0.0;
  x = createAVEC(a);
  f = inv_square(a);
  EXPECT_FLOAT_EQ(stan::math::positive_infinity(),f.val());

  f.grad(x,grad_f);
  EXPECT_FLOAT_EQ(stan::math::negative_infinity(),grad_f[0]);
}
Ejemplo n.º 2
0
 static inline T fun(const T& x) {
     return inv_square(x);
 }
Ejemplo n.º 3
0
TEST(AgradRev, check_varis_on_stack) {
    AVAR a = 7.0;
    test::check_varis_on_stack(inv_square(a));
}
Ejemplo n.º 4
0
 inline T0
 operator()(const T0& arg1) const {
     return inv_square(arg1);
 }