Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/andreasnoack/eiscor.jl

Julia wrapper of the Fortran library for efficiently solving structured matrix eigenvalue problems using unitary core transformations
https://github.com/andreasnoack/eiscor.jl

Last synced: 19 days ago
JSON representation

Julia wrapper of the Fortran library for efficiently solving structured matrix eigenvalue problems using unitary core transformations

Awesome Lists containing this project

README

        

# EiSCor
A Julia wrapper of the Fortran library [eiscor](https://github.com/jaurentz/eiscor)

## Installation:

```julia
julia> Pkg.clone("https://github.com/andreasnoack/EiSCor.jl")
julia> Pkg.build("EiSCor")
```

## Example:

```julia
julia> using EiSCor

julia> H = eye(5)[[5,1:4],:]
5x5 Array{Float64,2}:
0.0 0.0 0.0 0.0 1.0
1.0 0.0 0.0 0.0 0.0
0.0 1.0 0.0 0.0 0.0
0.0 0.0 1.0 0.0 0.0
0.0 0.0 0.0 1.0 0.0

julia> eigfact(EiSCor.UnitaryHessenbergMatrix(H))[:values]
5-element Array{Complex{Float64},1}:
-0.809017-0.587785im
-0.809017+0.587785im
0.309017-0.951057im
0.309017+0.951057im
1.0+0.0im

julia> eigfact(EiSCor.UnitaryHessenbergMatrix(H))[:vectors]
5x5 Array{Float64,2}:
-0.632456 1.33749e-16 -0.195469 0.601492 0.447214
0.511667 -0.371748 -0.632456 -3.05458e-5 0.447214
-0.19544 0.601501 -0.19541 -0.60151 0.447214
-0.19544 -0.601501 0.511685 -0.371723 0.447214
0.511667 0.371748 0.511649 0.371773 0.447214
```