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

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

#ifndef _EIGSUBSP_H
#define _EIGSUBSP_H

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

// -------------------------------------------------------------------------
// Subspace class:

class cEigenSubspace : public cEigenProblem
{
 public:
        cEigenSubspace( void );
       ~cEigenSubspace( 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 **);
  void StartingIteration( int, int , double **, double **, double ** );
  void Subspace         ( int, int *, double **, double **,
                          double *, double ** );
  int  Stop             ( int, int, double **, double ** );
  int  InverStop        ( int, double *, double * );
};

#endif
 

---