Exemplo n.º 1
0
    double eval() const
    {
      switch (deriv)
	{
	case 0:
	  return (rp | nhat) - ( Sigma + wallnHatPosition());
	case 1:
	  return (vp | nhat) - velnHatWall();
	case 2:
	  return Delta * Omega * Omega * std::cos(Omega * t); 
	default:
	  M_throw() << "Invalid access";
	}
    }
Exemplo n.º 2
0
 double F_firstDeriv() const
 {    
   return (vp | nhat) - velnHatWall();
 }
Exemplo n.º 3
0
 Vector wallVelocity() const
 {
   return  nhat * velnHatWall();
 }
Exemplo n.º 4
0
    virtual bool test_root() const
    {
      return (((vp | nhat) - velnHatWall()) 
	      * ((rp | nhat) - wallnHatPosition())) > 0;
    }