/**
 *  \file Ellipsoid3D.h   \brief Simple 3D ellipsoid class.
 *
 *  Copyright 2007-2010 IMP Inventors. All rights reserved.
 *
 */
#ifndef IMPALGEBRA_ELLIPSOID_3D_H
#define IMPALGEBRA_ELLIPSOID_3D_H

#include "Cylinder3D.h"
#include <cmath>
#include <IMP/constants.h>
#include "internal/cgal_predicates.h"
#include "BoundingBoxD.h"
#include "Transformation3D.h"

IMPALGEBRA_BEGIN_NAMESPACE

/** Represent an ellispoid in 3D.
    \geometry
  */
class IMPALGEBRAEXPORT Ellipsoid3D {
public:
  Ellipsoid3D(){
#if IMP_BUILD < IMP_FAST
    radii_[0]= std::numeric_limits<double>::quiet_NaN();
#endif
  }
  Ellipsoid3D(const VectorD<3>& center,double radius_x,
              double radius_y, double radius_z,
              const Rotation3D &rot);
  double get_radius(unsigned int i) const {
    IMP_USAGE_CHECK(!is_nan(radii_[0]),
              "Attempt to use uninitialized ellipsoid.");
    return radii_[i];
  }
  const VectorD<3> &get_center() const {return center_;}
  const Rotation3D &get_rotation() const {
    return rot_;
  }
  const Transformation3D get_transformation() const  {
    return Transformation3D(rot_, center_);
  }
  IMP_SHOWABLE_INLINE({
      out << "(" << spaces_io(center_) << ": " << get_radius(0)
          << ", " << get_radius(1) << ", " << get_radius(2)
          << ")";
    });
private:
  VectorD<3> center_;
  double radii_[3];
  Rotation3D rot_;
};

IMP_VOLUME_GEOMETRY_METHODS(Ellipsoid3D,
                            IMP_NOT_IMPLEMENTED,
                            return 4.0/3.0 * PI *g.get_radius(0)
                            *g.get_radius(1)*g.get_radius(2),
                            {
                              VectorD<3> v(g.get_radius(0),
                                         g.get_radius(1),
                                         g.get_radius(2));
                              BoundingBoxD<3> b(-v, v);
                              return get_transformed(b, g.get_transformation());
                            });

IMPALGEBRA_END_NAMESPACE

#endif /* IMPALGEBRA_ELLIPSOID_3D_H */