[Help-gsl] Re: gsl ode solver - const a problem?

[Heiko Bauke]

Hi,

On Fri, 18 Sep 2009 12:19:23 +0530
"Sumesh P.T." <address@hidden> wrote:

> I use gsl ode solver in my program for integration as given in :
> http://www.gnu.org/software/gsl/manual/html_node/ODE-Example-programs.html
> Look at the definition of func
> 
> func (double t, const double y[], double f[], void *params)
> 
> y is a const. I am facing an issue because of this.

I also think, this interface has an unfortunate design. In fact, I would
prefer 

func (double t, const double y[], double f[], const void *params)

There is no (obvious) reason why func should change the parameters. The
current design forces me regularly to apply a const_cast (in C++) when
solving ODEs with GSL routines.


        Heiko

-- 
-- Willst du den Wert des Geldes kennenlernen, geh und versuche 
-- dir welches zu borgen. (Benjamin Franklin, 1706-1790)
-- Cluster Computing @ http://www.clustercomputing.de
--       Heiko Bauke @ http://www.mpi-hd.mpg.de/personalhomes/bauke