This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]

Re: Diagonalization of complex matrices

To: Gerard Jungman <jungman at lanl dot gov>
Subject: Re: Diagonalization of complex matrices
From: Brian Gough <bjg at network-theory dot co dot uk>
Date: Tue, 26 Sep 2000 22:57:44 +0100 (BST)
Cc: gsl-discuss at sources dot redhat dot com
References: <20000925214016.A21923@pccp.com.ar><14801.1312.75613.220192@localhost><39D11321.888C3029@lanl.gov>
Reply-To: gsl-discuss at sources dot redhat dot com

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

References:
- Diagonalization of complex matrices
  - From: Pablo Bianucci
- Re: Diagonalization of complex matrices
  - From: Brian Gough
- Re: Diagonalization of complex matrices
  - From: Gerard Jungman

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]