00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #if !defined (octave_SCHUR_h)
00024 #define octave_SCHUR_h 1
00025
00026 #if defined (__GNUG__) && defined (USE_PRAGMA_INTERFACE_IMPLEMENTATION)
00027 #pragma interface
00028 #endif
00029
00030 #include <iostream>
00031 #include <string>
00032
00033 #include "dMatrix.h"
00034
00035 class
00036 SCHUR
00037 {
00038 public:
00039
00040 SCHUR (void)
00041 : schur_mat (), unitary_mat () { }
00042
00043 SCHUR (const Matrix& a, const std::string& ord, bool calc_unitary = true)
00044 : schur_mat (), unitary_mat () { init (a, ord, calc_unitary); }
00045
00046 SCHUR (const Matrix& a, const std::string& ord, int& info,
00047 bool calc_unitary = true)
00048 : schur_mat (), unitary_mat () { info = init (a, ord, calc_unitary); }
00049
00050 SCHUR (const SCHUR& a)
00051 : schur_mat (a.schur_mat), unitary_mat (a.unitary_mat) { }
00052
00053 SCHUR& operator = (const SCHUR& a)
00054 {
00055 if (this != &a)
00056 {
00057 schur_mat = a.schur_mat;
00058 unitary_mat = a.unitary_mat;
00059 }
00060 return *this;
00061 }
00062
00063 ~SCHUR (void) { }
00064
00065 Matrix schur_matrix (void) const { return schur_mat; }
00066
00067 Matrix unitary_matrix (void) const { return unitary_mat; }
00068
00069 friend std::ostream& operator << (std::ostream& os, const SCHUR& a);
00070
00071 typedef int (*select_function) (const double&, const double&);
00072
00073 private:
00074
00075 Matrix schur_mat;
00076 Matrix unitary_mat;
00077
00078 select_function selector;
00079
00080 int init (const Matrix& a, const std::string& ord, bool calc_unitary);
00081 };
00082
00083 #endif
00084
00085
00086
00087
00088
00089