Question: eigenvalue problem about owl HOT 2 CLOSED

Hi @wkwkes

you are right, there wasn't any yet, but it is very easy to add actually. I added four functions a moment ago and you can try them out if you are using docker image or build from master branch yourself.

The four functions in Linalg module are

val eigen_symm : mat_d -> mat_d    (* real symmetric, returns eigenvalues *)
val eigen_symmv : mat_d -> mat_d * mat_d    (* real symmetric, returns eigenvalues and vectors *)
val eigen_nonsymm : mat_d -> mat_z    (* real asymmetric, returns eigenvalues *)
val eigen_nonsymmv : mat_d -> mat_z * mat_z    (* real symmetric, returns eigenvalues and vectors *)

Here is an example:

let x = Mat.semidef 4;;
let v, z = Linalg.eigen_symmv x;;

let z' = Linalg.inv z;;
Mat.(z' $@ x $@ z);;

I hope this suits your need for the time being. The linalg module will be redesigned soon in the future to provide better apis :)

from owl.

This is what I wanted.
Thank you for the so quick response!

from owl.