newarp_UpperHessenbergEigen_bones.hpp (armadillo-10.8.2.tar.xz) | : | newarp_UpperHessenbergEigen_bones.hpp (armadillo-11.0.0.tar.xz) | ||
---|---|---|---|---|
skipping to change at line 22 | skipping to change at line 22 | |||
// distributed under the License is distributed on an "AS IS" BASIS, | // distributed under the License is distributed on an "AS IS" BASIS, | |||
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. | // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. | |||
// See the License for the specific language governing permissions and | // See the License for the specific language governing permissions and | |||
// limitations under the License. | // limitations under the License. | |||
// ------------------------------------------------------------------------ | // ------------------------------------------------------------------------ | |||
namespace newarp | namespace newarp | |||
{ | { | |||
//! Calculate the eigenvalues and eigenvectors of an upper Hessenberg matrix. | //! Calculate the eigenvalues and eigenvectors of an upper Hessenberg matrix. | |||
//! This class is a wrapper of the Lapack functions `_lahqr` and `_trevc`. | //! This class is uses lapack::lahqr() and lapack::trevc() | |||
template<typename eT> | template<typename eT> | |||
class UpperHessenbergEigen | class UpperHessenbergEigen | |||
{ | { | |||
private: | private: | |||
blas_int n; | uword n_rows; | |||
Mat<eT> mat_Z; // In the first stage, H = ZTZ', Z is an ortho gonal matrix | Mat<eT> mat_Z; // In the first stage, H = ZTZ', Z is an ortho gonal matrix | |||
// In the second stage, Z will be overwritten by the eigenvectors of H | // In the second stage, Z will be overwritten by the eigenvectors of H | |||
Mat<eT> mat_T; // H = ZTZ', T is a Schur form matrix | Mat<eT> mat_T; // H = ZTZ', T is a Schur form matrix | |||
Col< std::complex<eT> > evals; // eigenvalues of H | Col< std::complex<eT> > evals; // eigenvalues of H | |||
bool computed; | bool computed; | |||
public: | public: | |||
//! Default constructor. Computation can | //! Default constructor. Computation can | |||
//! be performed later by calling the compute() method. | //! be performed later by calling the compute() method. | |||
End of changes. 2 change blocks. | ||||
2 lines changed or deleted | 2 lines changed or added |