17 #include <deal.II/base/geometric_utilities.h> 18 #include <deal.II/base/exceptions.h> 21 DEAL_II_NAMESPACE_OPEN
31 <<
"The radius <" << arg1 <<
"> can not be negative.");
35 <<
"The azimuth angle <" << arg1 <<
"> is not in [0,2Pi).");
39 <<
"The polar angle <" << arg1 <<
"> is not in [0,Pi].");
43 std::array<double,dim>
46 std::array<double,dim> scoord;
49 scoord[0] = position.
norm();
51 scoord[1] = std::atan2(position(1),position(0));
60 if (scoord[0] > std::numeric_limits<double>::min())
61 scoord[2] = std::acos(position(2)/scoord[0]);
68 template <std::
size_t dim>
84 ccoord[0] = scoord[0] * std::cos(scoord[1]);
85 ccoord[1] = scoord[0] * std::sin(scoord[1]);
93 ccoord[0] = scoord[0] * std::sin(scoord[2]) * std::cos(scoord[1]);
94 ccoord[1] = scoord[0] * std::sin(scoord[2]) * std::sin(scoord[1]);
95 ccoord[2] = scoord[0] * std::cos(scoord[2]);
107 #include "geometric_utilities.inst" 112 DEAL_II_NAMESPACE_CLOSE
numbers::NumberTraits< Number >::real_type norm() const
static ::ExceptionBase & SphericalPolar(double arg1)
#define DeclException1(Exception1, type1, outsequence)
std::array< double, dim > to_spherical(const Point< dim > &point)
#define Assert(cond, exc)
static ::ExceptionBase & SphericalAzimuth(double arg1)
static ::ExceptionBase & NegativeRadius(double arg1)
static ::ExceptionBase & ExcNotImplemented()
Point< dim > from_spherical(const std::array< double, dim > &scoord)