Eigs.H

Go to the documentation of this file.

// $Revision: 1.6 $  $Author: schwehr $  $Date: 2004/11/03 00:29:23 $
/*
    Copyright (C) 2004  Kurt Schwehr

    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

*/
#ifndef _EIGS_H_
#define _EIGS_H_

#include <gsl/gsl_eigen.h> 

#include <vector>



enum EigsEnum {KMAX=1, KINT=2, KMIN=3};




size_t GetIndex(const EigsEnum which, const gsl_vector *v);


bool GetEig(const EigsEnum which, const gsl_matrix *eigenvecs, const gsl_vector *eigenvals,
            float vec[3], float &val);

bool GetEigs(const gsl_matrix *eigenvec, const gsl_vector *eigenval, float newEigs[9]);



class S_Engine {
public:

  S_Engine();

  ~S_Engine();  

  S_Engine(const std::vector<float> &s);

  bool setS(const std::vector<float> &s);

  bool getXYZ(const EigsEnum which, std::vector<float> &xyz) const;

  bool isValid() const {return(valid);}

  // FIX: could implement these...
  //S_Engine(const float s[6]);
  //S_Engine(const float s[7]);
  //bool getS(std::vector<float> &s);
  // Lisa Style ???
  //bool getEigs(std::vector<float> &V, vector<float> &T); ///< min, int, max sorted
  //bool getEig(const EigsEnum which, std::vector<float> &V, float &T);

  //bool getTau(const EigsEnum which, float &T);
  //bool getXYZ(const EigsEnum which, float &x, float &y, float &z);
  //bool getPTR(std::vector<float> &ptr);
  //bool getPTR(float &phi, float &theta, float &radius);
  //bool getLDI(float &length, float &dec, float &inc);
private:
  bool init(); 
  bool valid; 
  std::vector<float> s; 

  gsl_eigen_symmv_workspace *w; 
  gsl_matrix *A; 
  gsl_vector *eigenval; 
  gsl_matrix *eigenvec; 

  float eigs[9];

  //vector<float> xyz; ///< size==9, min, int, max triples
  //vector<float> ptr; ///< size==9, min, int, max triples
  //vector<float> eigs; ///< same format as s_eigs
};

#endif // _EIGS_H_

---