77 scalar rprime = radius_ + local.
x()*
sin(phi);
79 if ((local.
x()*
sin(phi)) > (radius_))
82 <<
"Badly defined toroidal coordinates"
127 "toroidalCS::globalToLocal(const vector&, bool) const"
142 "toroidalCS::globalToLocal(const vectorField&, bool) const"
152 os <<
"radius: " << radius() <<
endl;