* Diagonalization of complex matrices
@ 2000-09-25 17:41 Pablo Bianucci
2000-09-26 13:22 ` Brian Gough
0 siblings, 1 reply; 4+ messages in thread
From: Pablo Bianucci @ 2000-09-25 17:41 UTC (permalink / raw)
To: gsl-discuss
Hi!
The last time I checked (ver 0.6) I could not find any routines to diagonalize
complex matrices in GSL. Are there any? If there are not, what may I do to
diagonalize a complex matrix? (More specifically, I want to calculate its
logarithm).
Thanks a lot!
Pablo B.
^ permalink raw reply [flat|nested] 4+ messages in thread
* Re: Diagonalization of complex matrices
2000-09-25 17:41 Diagonalization of complex matrices Pablo Bianucci
@ 2000-09-26 13:22 ` Brian Gough
2000-09-26 14:14 ` Gerard Jungman
0 siblings, 1 reply; 4+ messages in thread
From: Brian Gough @ 2000-09-26 13:22 UTC (permalink / raw)
To: Pablo Bianucci; +Cc: gsl-discuss
At the moment we don't have any routines for complex matrices. I
would recommend using LAPACK (or CLAPACK) from netlib for this. There
are precompiled LAPACK packages available for most GNU/Linux
distributions.
Brian
Pablo Bianucci writes:
> The last time I checked (ver 0.6) I could not find any routines to
> diagonalize complex matrices in GSL. Are there any? If there are
> not, what may I do to diagonalize a complex matrix? (More
> specifically, I want to calculate its logarithm).
^ permalink raw reply [flat|nested] 4+ messages in thread
* Re: Diagonalization of complex matrices
2000-09-26 13:22 ` Brian Gough
@ 2000-09-26 14:14 ` Gerard Jungman
2000-09-26 15:01 ` Brian Gough
0 siblings, 1 reply; 4+ messages in thread
From: Gerard Jungman @ 2000-09-26 14:14 UTC (permalink / raw)
To: gsl-discuss
Brian Gough wrote:
>
> At the moment we don't have any routines for complex matrices. I
> would recommend using LAPACK (or CLAPACK) from netlib for this. There
> are precompiled LAPACK packages available for most GNU/Linux
> distributions.
Of course, how he is supposed to do this while using gsl_matrix
is unclear. What exactly does "use lapack" mean in this context?
^ permalink raw reply [flat|nested] 4+ messages in thread
* Re: Diagonalization of complex matrices
2000-09-26 14:14 ` Gerard Jungman
@ 2000-09-26 15:01 ` Brian Gough
0 siblings, 0 replies; 4+ messages in thread
From: Brian Gough @ 2000-09-26 15:01 UTC (permalink / raw)
To: Gerard Jungman; +Cc: gsl-discuss
I should have mentioned that the matrix conventions going from GSL to
LAPACK are transposed (C vs Fortran -- see the matrix chapter in the
manual for details) but taking this into account a direct call should
work fine.
e.g. to diagonalize a complex hermitian matrix,
zheev ('V', 'L', m->size1, m->data, m->tda, w, work, lwork, rwork, &info)
swapping 'L' for 'U' to take account of the transpose and keeping in
mind that the result will also be transposed.
then link with -lgsl -llapack ... etc
^ permalink raw reply [flat|nested] 4+ messages in thread
end of thread, other threads:[~2000-09-26 15:01 UTC | newest]
Thread overview: 4+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2000-09-25 17:41 Diagonalization of complex matrices Pablo Bianucci
2000-09-26 13:22 ` Brian Gough
2000-09-26 14:14 ` Gerard Jungman
2000-09-26 15:01 ` Brian Gough
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).