[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] vector views

From: Rhys Ulerich
Subject: Re: [Help-gsl] vector views
Date: Mon, 7 May 2012 08:48:44 -0500

> Just wondering if someone can help me out with using vector views.
> I have a set of N complex numbers which I am storing in a 2*N size gsl_vector.
> ...
> So here's my question, does anyone have any tips on how to set up vector
> views so that I have pointers available that logically map to the real and 
> imaginary
> portions of the vector?

I'd create a complex-valued vector of length N using
gsl_vector_complex_alloc (say "z"), then use gsl_vector_view to create
a real-valued view of length 2N from the data (say "zd"), then use
gsl_vector_subvector_with_stride to create strided views for the real
and imaginary parts.  As you'd expect the real-valued offset will be
zero with stride two and the imaginary-valued offset is one with
stride two.  The views are documented well at
The complex-valued vector details are analogous and are declared
within <gsl/gsl_vector_complex_double.h> in a readable fashion.

You could, of course, skip my suggestion to start with a
complex-valued vector and simply start with the 2*N gsl_vector.  I'm
guessing having a complex-valued vector will be handy for other
reasons in your code.

Hope that helps,

reply via email to

[Prev in Thread] Current Thread [Next in Thread]