J:/Tecgraf/prj/femoop/trunk/eig/eigjacob.h

// --------------------------------------------------------------------------
//
// EigJacob.h - This file contains routines for solving the eigenvalue
//              problem in a general form, i.e.: K x = l M x, by using
//              JACOBI transformations and Inverse iterations.
//
// --------------------------------------------------------------------------
//
// $Author: joaquim $
// $Revision: 1.2 $
// $Date: 1998/02/16 20:21:11 $
// $State: Exp $
//
// --------------------------------------------------------------------------

#ifndef _EIGJAC_H
#define _EIGJAC_H

#include "gbldefs.h"
#include "eig.h"

// -------------------------------------------------------------------------
// Jacobi class:

class cEigenJacobi : public cEigenProblem
{
 public:
        cEigenJacobi( void );
       ~cEigenJacobi( void );
  void  EigenValue  ( int, int *, double **, double **, double *, double ** );
  
 private:
  void QR              ( int, double **, double * );
  void InverseIteration( int, int *, double **, double *, double ** );
  void GramSchmidt     ( int, int, double *, double **);
  int  InverStop       ( int, double *, double * );

};


#endif 

---