- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi guys;
I'm trying to compute "dynamic mode shapes" in my fortran code, so I need to find either eigenvalues of a symmetric sparse matrix or the inverse of that matrix (the first method is prefered).
Is there any routine in BLAS or any spesific procedure that could help me?
thanks in advance
Hossein
Link Copied
5 Replies
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hossein, there is no such functionality for the sparse routines at all.
We have the similar functionality in the current version for the dense matrix only.
--Gennady
We have the similar functionality in the current version for the dense matrix only.
--Gennady
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Quoting - Gennady Fedorov (Intel)
Hossein, there is no such functionality for the sparse routines at all.
We have the similar functionality in the current version for the dense matrix only.
--Gennady
We have the similar functionality in the current version for the dense matrix only.
--Gennady
Dear Gennady; thanks for your quick reply. I found a way to generate inverse of my sparse matrix. Hopefully, it will do the job.
Regards;
Hossein
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Quoting - pourmatin85
Dear Gennady; thanks for your quick reply. I found a way to generate inverse of my sparse matrix. Hopefully, it will do the job.
Regards;
Hossein
Hossein, may I ask you, how did you solve this problem?
--Gennady
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Quoting - Gennady Fedorov (Intel)
Hossein, may I ask you, how did you solve this problem?
--Gennady
Well, suppose B is a (N,1) vector. Every components of B is set to zero but it's Ith component, which is 1. The resulting vector of the product of a matrix, say K, and B is the Ith column of K.
So, all you need to do, to find the inverse of K, is to solve N equation systems with PARDISO. It means that in the Ith iteration out of N iterations, the Ith component of B is 1, and the result is the Ith column of the inverse of K:
x*K=B -> x=inv(K)*B -> x=inv(K(1:N,I))
By the way, MKL doesn't have any routine for matrix-matrix product, in which both matrices are sparse (is that right?). Isn't it weird!! I have to do the same trick to find the answer of mymatrix-matrix products.
Hossein
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Quoting - pourmatin85
By the way, MKL doesn't have any routine for matrix-matrix product, in which both matrices are sparse (is that right?). Isn't it weird!! I have to do the same trick to find the answer of mymatrix-matrix products.
see this thread.
A.
Reply
Topic Options
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page