This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: gsl_multifit
Kai Trukenmueller writes:
> I benchmarked the programs on system with Athlon 1GHz, 640MB running Linux-2.4.10 with the following results:
>
> Dimensions (120,80) (300,200) (600,400) (900,600) (1200,800)
> Octave-Script 0.060s 0.410s 3.390s 12.260s 39.790s
> Compiled gsl 0.110s 2.240s 31.310s 2m7.590s seg.fault
>
> The gsl routine not just turns out to be much slower, it also seems to
> be `instable' for higher orders (the results do not convege -> inft..)
> Maybe sth. is wrong in my code. For low-orders (~<500) both results are
> equivalent.
Regarding the segfault, your program contains some fixed length arrays
which are accessed outside the array bounds -- if you allocate them
dynamically that problem should go away.