// EigenvalueVector.h		計(jì)算特征值特征向量頭文件
// Ver 1.0.0.0
// 版權(quán)所有(C) 何渝, 2002

// 最后修改: 2002.8.10

#ifndef _EIGENVALUEVECTOR_H		//避免多次編譯
#define _EIGENVALUEVECTOR_H

#include "Matrix.h"		//矩陣類及相關(guān)函數(shù)等的定義
#include <comm.h>		//公共頭文件
#include <math.h>		//數(shù)學(xué)頭文件

using namespace std;	//名字空間

//約化對(duì)稱陣為對(duì)稱三對(duì)角陣的豪斯荷爾德變換法
template <class _Ty>
int HouseholderTransform(matrix<_Ty>& a, matrix<_Ty>& q, 
									valarray<_Ty>& b, valarray<_Ty>& c);

//實(shí)對(duì)稱三角陣全部特征值及特征向量QR法
template <class _Ty>
int EigenvalueVectorRealTriangleQR(valarray<_Ty>& b, 
					valarray<_Ty>& c, matrix<_Ty>& q, _Ty eps, int l);

//約化一般實(shí)矩陣為赫申伯格陣的初等相似變換法
//矩陣類型應(yīng)是浮點(diǎn)型
template <class _Ty>
int HessenbergTransform(matrix<_Ty>& a);

//求赫申伯格陣全部特征值QR法
template <class _Ty>
int EigenvalueVectorHessenbergQR(matrix<_Ty>& a,  
						valarray<complex<_Ty> >& uv, _Ty eps, int jt);

//實(shí)對(duì)稱陣特征值及特征向量雅可比法
template <class _Ty>
int EigenvalueVectorRealSymmetryJacobi(matrix<_Ty>& a,  
										matrix<_Ty>& v, _Ty eps, int jt);

//實(shí)對(duì)稱陣特征值及特征向量雅可比過關(guān)法
template <class _Ty>
int EigenvalueVectorRealSymmetryJacobiB(matrix<_Ty>& a,  
												matrix<_Ty>& v, _Ty eps);

#include "EigenvalueVector.inl"		//類及相關(guān)函數(shù)的定義頭文件

#endif		// _EIGENVALUEVECTOR_H