CODEWITHSUNDEEP
c#lapackintel-mklmathnet-numerics
Erik
Erik

Reputation: 954

Solve general eigenproblem with mathnet

I found a very promissing package (MatNet) to work with in my structural engineering program. However, while MathNet relies heavily on MKL from Intel, I cannot see how I can solve a general eigenproblem to find the structure's eigenfrequencies. Lapack has this routine,so MKL should have it too. Why not MathNet? Or: the MKL/Lapack routine seems to be not exposed to C#. Anyone can point me in the right direction?

Upvotes: 0

Views: 526

Answers (1)

Erik
Erik

Reputation: 954

Once you know how to do it, it is easy:

Convert generalized eigenproblem to standard form:

K⋅Φ  = Λ⋅M⋅Φ  →   K_⋅Φ  = Λ⋅Φ

    
Matrix C = ((Matrix)M.Cholesky().Factor);                                                        
Matrix K_ = (Matrix)C.Inverse().Multiply(K).Multiply(C.Transpose().Inverse());

From there, ust apply other MathNet methods to solve the standard problem. No more LAPACK, Intel MKL, just MathNet

Upvotes: 1

Related Questions