Re: [Rd] transpose of complex matrices in R

2010-07-30 Thread William Dunlap
> -Original Message- > From: r-devel-boun...@r-project.org > [mailto:r-devel-boun...@r-project.org] On Behalf Of Olaf Mersmann > Sent: Friday, July 30, 2010 3:01 AM > To: Robin Hankin > Cc: R-devel@r-project.org; Martin Maechler > Subject: Re: [Rd] transpose of complex matrices in R > > H

Re: [Rd] transpose of complex matrices in R

2010-07-30 Thread Olaf Mersmann
Hi, On 30.07.2010, at 11:35, Robin Hankin wrote: > 3. Try to define a t.complex() function: > t.complex <- function(x){t(Conj(x))} > (also fails because of recursion) Try this version: t.complex <- function(x) { xx <- Conj(x) .Internal(t.default(xx)) } You get infinite recursion in

[Rd] transpose of complex matrices in R

2010-07-30 Thread Prof. John C Nash
When one is working with complex matrices, "transpose" very nearly always means *Hermitian* transpose, that is, A[i,j] <- Conj(A[j,i]). One often writes A^* for the Hermitian transpose. I believe that I have actually (many years ago) used a true complex transpose, but I agree that one more of

Re: [Rd] transpose of complex matrices in R

2010-07-30 Thread Robin Hankin
Hello Peter thanks for this. On 07/30/2010 11:01 AM, peter dalgaard wrote: What's wrong with t.complex<- function(x) t.default(Conj(x)) M<- matrix(rnorm(4)+1i*rnorm(4),2) M It's not going to help with the cross products though. As a general matter, in my book, transpose is tran

Re: [Rd] transpose of complex matrices in R

2010-07-30 Thread Duncan Murdoch
Robin Hankin wrote: Hello everybody When one is working with complex matrices, "transpose" very nearly always means *Hermitian* transpose, that is, A[i,j] <- Conj(A[j,i]). One often writes A^* for the Hermitian transpose. I have only once seen a "real-life" case where transposition does not

Re: [Rd] transpose of complex matrices in R

2010-07-30 Thread peter dalgaard
On Jul 30, 2010, at 11:35 AM, Robin Hankin wrote: > Hello everybody > > When one is working with complex matrices, "transpose" very nearly always > means > *Hermitian* transpose, that is, A[i,j] <- Conj(A[j,i]). > One often writes A^* for the Hermitian transpose. > > I have only once seen a

[Rd] transpose of complex matrices in R

2010-07-30 Thread Robin Hankin
Hello everybody When one is working with complex matrices, "transpose" very nearly always means *Hermitian* transpose, that is, A[i,j] <- Conj(A[j,i]). One often writes A^* for the Hermitian transpose. I have only once seen a "real-life" case where transposition does not occur simultaneously