#ifndef __AUTILS_HH__
#define __AUTILS_HH__

#include "AAtom.hh"
#include "AFigure.hh"

class AUtils
{
public:
  virtual double fitPoints(double rot[3][3],double trans[3],
			   AAtom** set1,AAtom** set2,const int n_points) = 0;
  virtual double eRMSD(AFigure *f1, AFigure *f2) = 0;
  virtual double reducedDet(double x1,double y1,double z1,     //| x1 y1 z1 1 |
			    double x2,double y2,double z2,     //| x2 y2 z2 1 |
			    double x3,double y3,double z3,     //| x3 y3 z3 1 |
			    double x4,double y4,double z4) = 0;//| x4 y4 z4 1 |
  virtual int    tetraType(AFigure *tet) = 0;
  virtual bool   atomInsideTetra(AAtom *atom,AFigure *tet) = 0;
  virtual bool   getCircumeSphere(AFigure *tet,
				  double &x,double &y,double &z,double &R) = 0;
  virtual bool   getCircumeSphere(double x1,double y1,double z1,
				  double x2,double y2,double z2,
				  double x3,double y3,double z3,
				  double x4,double y4,double z4,
				  double &x,double &y,double &z,double &R) = 0;
  virtual double erf(double x)  = 0;
  virtual double erfc(double x) = 0;
};

#endif