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

// --------------------------------------------------------------------------
//
// Eig.h - This header file contains the public data structure definitions
//         for the eigenproblems class algorithms mechanism within FEMOOP.
//
// --------------------------------------------------------------------------
//
// $Author: joaquim $
// $Revision: 1.2 $
// $Date: 1998/02/16 20:21:11 $
// $State: Exp $
//
// --------------------------------------------------------------------------

#ifndef _EIG_H
#define _EIG_H

//#include "gbldefs.h"

#define TOL_EIG 1.0e-06


// -------------------------------------------------------------------------
//  Eigen types:

typedef enum _eigentype
{
 NONE_EIG,
 JACOBI,
 HQR,
 SUBSPACE
} eEigenType;


// -------------------------------------------------------------------------
// Base EigenProblem class:

class cEigenProblem
{
 protected:
  static int        _iFirstMode;        // First mode shape to be evaluated   
  static int        _iLastMode;         // Last  mode shape to be evaluated 
         eEigenType _eType;             // Type of algorithm
 
 public:
                         cEigenProblem( void );
  virtual               ~cEigenProblem( void );
  virtual void           EigenValue   ( int, int *, double **, double **,
                                        double *, double ** ) { };
  static  cEigenProblem *Read         ( void );
  static  void           Free         ( void ) { };
          int            GetEigenType ( void ) { return( (int)_eType ); }

 protected:
  virtual void QR               ( int, double **, double * ) { };
          void ComputeStandBase ( int, int *, double **, double **,
                                  double **, int * );
          void MountHhMatrix    ( int, double, double *, double *, double ** );
          void MountJacobiMatrix( int, int, int, double **, double [2][2] );
          void MountTotalMatrix ( int, int *, double **, double ** );
          void NormalizeSort    ( int, int, int, double *, double ** );
          int  QRStop           ( int, double ** );

 private:
  static cEigenProblem *ProcessEigenProblemAlgorithm( void );
  static void           ProcessCurves               ( void );
  static void           ReadModes                   ( void );

};

#endif 

---