+inline bool HyperSphere<T, D>::contains(const LinAl::Vector<T, D> &point) const
+{
+ return inner_product(point, point)<=radius*radius;
+}
+
+template<typename T, unsigned D>
+inline unsigned HyperSphere<T, D>::get_intersections(const Ray<T, D> &ray, SurfacePoint<T, D> *points, unsigned size) const