Re: Non square matrix transpose ?

[Brian Gough <bjg at network-theory dot co dot uk>]

Hi,

As you found, there has not been a function to do rectangular
transposes.

I have added a function gsl_matrix_transpose_memcpy(dest,src) to
handle the rectangular case.   It is in CVS now.

It operates out of place. (BLAS matrices use a unit stride in one
direction, so in general the in-place transpose of a submatrix view
cannot itself be represented as a submatrix view).

regards
Brian Gough

Jose Miguel Buenaposada Biencinto writes:
 > Hello !
 > 
 > First of all I have to ask you to answer to me directly because I'm not
 > in the list.
 > 
 > I'm not a mathematician, only a poor computer scientist ;-), so I'm
 > pretty sure to make a trivial question. Anyway I need some help
 > with gsl.
 > 
 > What is the best way to get the transpose of a non square matrix using
 > gsl ? I can figure out how to do it using too much gsl matrix functions
 > in
 > a loop but, maybe there is a better solution. I'm not familiar with BLAS
 > 
 > and I only want to use gsl for some machine vision algoritms so I can
 > not understand the reason to only provide an in place matrix transpose
 > function only suitable for square matrices.
 > 
 > Can anybody explain me this subject ?
 > 
 > Thanks for this great library.
 > Regards.
 > 
 > --
 > -------------------------------------------------------------------
 > Jose Miguel Buenaposada Biencinto
 > Facultad de Informática                  | Computer Science School
 > Universidad Politécnica de Madrid  (UPM) | Madrid Tech. University.
 > Campus de Montegancedo s/n
 > 28660 Boadilla del Monte, Madrid (ESPAÑA)| (SPAIN)
 > Voice    +34 91 336 69 47
 > e-mail:  jmbuena@dia.fi.upm.es
 > web   :  http://www.dia.fi.upm.es/~jmbuena
 > -------------------------------------------------------------------
 > 
 > 
 > 
 >