[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Question: gsl in C++ ??
From: |
eknecronzontas |
Subject: |
Re: [Help-gsl] Question: gsl in C++ ?? |
Date: |
Tue, 13 Mar 2007 19:25:13 -0700 (PDT) |
Hello!
In regards to Jordi's message about indexing with operator[]:
> they're indexing with
> operator[]. Not sure, but as I recall, this doesn't permit indexing
> such as a[1,2] since operator[] can only take one argument; the
> numeric community has a Fortran heritage of indexing with operator()
> instead which does allow more than one argument in C++ (here's our
> chance to start indexing from 1 à la GNU Octave!) We should also look
> into implementing gsl_vector with valarrays instead of direct memory
> manipulation, but this might break compatibility with BLAS.
I use operator[] because it allows me to use C-style arrays, which are
faster, and can be easier to optimize with for the compiler. Because of this,
some of the Hal versions of the GSL routines can be faster than their GSL
progenitors (when they are written correctly). I have tested this for a few
particular cases, but some things are missing. (This also means I need to have
separate template types for the array allocation, which is a bit confusing at
first glance.) The situation for 2-d or higher-d arrays is more complicated...
If you're not constrained to the GSL structure, then there are C++ ish
alternatives that have already been created (e.g. the template numerical
toolkit aka TNT).
Take care,
Andrew
Re: [Help-gsl] Question: gsl in C++ ??, Brian Gough, 2007/03/13
Re: [Help-gsl] Question: gsl in C++ ??, Alfredo Solano, 2007/03/14