2
$\begingroup$

I'm trying to find eigenvalues of symbolic hermitian 40x40 block matrix. It is made from 20 2x2 matrices like this:

ArrayFlatten[Table[M[m,n],{m,1,20},{n,1,20}]]

My computer can't solve this problem, it takes too much time, but I have an access to a multiple core cluster in my university. In order too use it, I have to parallelize the evaluation. Mathematica doesn't it automatically, and the command Parallelize[] result in an error Parallelize::nopar1.

So, the question is, is there a way to parallelize this task?

$\endgroup$
  • $\begingroup$ I have the impression that the general process used to find eigenvalues is not parallelizable at all. It may be that special cases can be parallelizable, but you have to work out the details yourself and explicitly handle parallelization of intemediate calculations yourself. BTW you say you have 20 2x2 matrices, but your code suggests that you have 20x20 of such matrices (unless, of course, many of those 400 matrices are identical. $\endgroup$ – Sjoerd C. de Vries Jul 24 '15 at 10:20
  • 1
    $\begingroup$ Also large symbolic matrices quickly yield extremely complex results, if they can be found at all. $\endgroup$ – Sjoerd C. de Vries Jul 24 '15 at 10:28
  • 4
    $\begingroup$ Yes, if you were to obtain a symbolic result, it's likely it would be too complicated for you to do anything useful with it. $\endgroup$ – Simon Rochester Jul 24 '15 at 10:29
  • $\begingroup$ Just to get a flavour of what you may face, look here - mathematica.stackexchange.com/questions/2586/… and it is a 10x10 matrix ! $\endgroup$ – Sumit Jul 24 '15 at 12:13
  • 1
    $\begingroup$ Does the matrix have any structure? If it does, and if you know its symmetry group, you may be able exploit that to block diagonalize your matrix. Then, each block is separately diagonalizable. This might only work if your matrix is symmetric/Hermitian, but I don't know. $\endgroup$ – rcollyer Jul 24 '15 at 14:10

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.