Re: Matrix/Vector: Docu typo, feature request, & code

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

Hans Ekkehard Plesser writes:
 >  Hi!
 >  First, a tiny bug in the Matrix documentation.  In the
 > documentation of gsl_matrix_submatrix the function prototype is
 > given as
 >  gsl_matrix_view gsl_matrix_submatrix (gsl_matrix * m, size_t i,
 > size_t j, size_t n1, size_t n2)
 >  but subsequent text talks about indices k1 and k2; obviously one
 > should replace i:=k1, j:=k2.

Thanks.

 >
 >  Now a feature request: I would like to see functions to copy a
 > (contiguous) part of a matrix/vector into an arbitrary position of
 > another matrix/vector.  I append some code for double
 > matrices/vectors, implementing this functionality based on gsl
 > highlevel functions.  Maybe someone with more experiences could
 > re-code this in a more "direct" manner and generalize it to all
 > data types (I found the vector/matrix code a bit too confusing to
 > mess with it directly).

I think it's better to not put this in GSL itself, since it's a
combination of existing GSL view functions, so it's not orthogonal.
Probably best to put it in the application since the function is
small.

Brian