---

EllipsoidalAnomaly.h

/****************************************************************************
 *                                                                          *
 *  Program: FDDlib                                                         *
 *  Version: 1.1                                                            *
 *                                                                          *
 *  Copyright (C) 2002 - 2004 by Eric Miller and Dana Brooks                *
 *  All rights reserved.                                                    *
 *                                                                          *
 *  This software is Version 1.1 of the fddlib tomography toolbox.          *
 *  It is not to be redistributed or used for any commercial purpose        *
 *  without the prior written consent of the authors and Northeastern       *
 *  University.                                                             *
 *                                                                          *
 *  This software is provided as is, and any express or implied warranties, *
 *  including but not limited to the implied warranty of merchantability    *
 *  and the implied warranty of fitness for a particular purpose, are dis-  *
 *  claimed.  In no event shall the authors or Northeastern University be   *
 *  liable for any direct, indirect, incidental, special, exemplary, or     *
 *  consequential damages (including but not limited to procurement of      * 
 *  substitute goods or services; loss of use, data, or profits; or busi-   *
 *  ness interruption) however caused and on any theory of liability,       *
 *  whether in contract, strict liability, or tort (including negligence or *
 *  otherwise) arising in any way from the use of this software, even if    *
 *  advised of the possibility of such damage.                              *
 *                                                                          *
 *  Portions of this code benefit from ideas from Kyle Guilbert, Greg       *
 *  Boverman, Derek Uluski, David Kaeli, and Jennifer Black                 *
 *                                                                          *
 ****************************************************************************/

#ifndef _ELLIPSOIDAL_ANOMALY_H_
#define _ELLIPSOIDAL_ANOMALY_H_

#include "Anomaly3D.h"
#include "CartesianNode3D.h"
#include "Property.h"
#include <string>
#include "DenseMatrix.h"

namespace FDDlib {

class EllipsoidalAnomaly: public Anomaly3D
{
protected:

  Location3D center_;
  
  double xlen_;

  double ylen_;

  double zlen_;  

  double theta1_;
  
  double theta2_;
  
  double theta3_;
  
  DenseMatrix<double> D;

  DenseMatrix<double> U;
  
  DenseMatrix<double> UD_neg2U_t;
  
  void updateData();

public:
  
  EllipsoidalAnomaly(const Location3D& center,
                     double xlen, double ylen, double zlen,
                     double theta1, double theta2, double theta3,
                     Property* prop) throw(std::string);
                     
  Location3D getCenter() const;

  void setCenter(const Location3D& center);
  
  std::vector<double> getLengths() const;
  
  void setLengths(double xlen, double ylen, double zlen) throw(std::string);
  
  std::vector<double> getEulerAngles() const;
  
  void setEulerAngles(double theta1, double theta2, double theta3);
  
  bool encloses(const CartesianNode3D& node) const;
  
};  
}
#endif

---