Re: [Help-gsl] GSL performance compared to C++ alternatives such as Armadillo

2015-10-26 Thread Rhys Ulerich
> Armadillo would be faster. > GSL would be more reliable and wider used. Eigen may be another good alternative for C++. - Rhys

Re: [Help-gsl] GSL performance compared to C++ alternatives such as Armadillo

2015-10-26 Thread Vladimir Sakharuk
Subject: [Help-gsl] GSL performance compared to C++ alternatives such as Armadillo Hello all, First of this is not a flame post, I would like to use GSL as linear algebra library for Neural Network training, so this would involve calculating hessian matrices, normal matrix and vector operations etc

[Help-gsl] GSL performance compared to C++ alternatives such as Armadillo

2015-10-26 Thread Christopher Choi
Hello all, First of this is not a flame post, I would like to use GSL as linear algebra library for Neural Network training, so this would involve calculating hessian matrices, normal matrix and vector operations etc... I was curious as to how does GSL compares to C++ alternatives such as Armadill

Re: [Help-gsl] gsl performance

2013-10-08 Thread onefire
I know that its all "just RAM". What I wanted to say (but did not) was that I've seen code become much faster when stuff is allocated on the stack but that is because compilers can sometimes do better "magic" in such cases. I did measure running time with a loop that reuses a minimizer and the res

Re: [Help-gsl] gsl performance

2013-10-08 Thread Sam Mason
On 7 October 2013 22:02, onefire wrote: > Sam Mason wrote: > "I tried this technique with the ODE solvers in the GSL and it gave me > about 5% overall performance improvement so I dropped it from my code, > it was a fiddle to maintain and the user would barely notice the > difference. I was doing

Re: [Help-gsl] gsl performance

2013-10-07 Thread onefire
"I tried this technique with the ODE solvers in the GSL and it gave me about 5% overall performance improvement so I dropped it from my code, it was a fiddle to maintain and the user would barely notice the difference. I was doing quite a lot of other things, so maybe if your overall time is domin

Re: [Help-gsl] gsl performance

2013-10-07 Thread Sam Mason
Hi, On 7 October 2013 18:22, onefire wrote: > One, which I like, is that the routines give the user low-level control > over their progress in the sense that you can create an object, manually > iterate and observe their progress. I prefer this over having a single > function that does the work u

Re: [Help-gsl] gsl performance

2013-10-07 Thread onefire
If there is a chance of a patch getting accepted I will certainly submit it if I have one. I took a look at the code in the multimin directory last night (the last problem affects other parts of the library but my example was about the multidimensional minimizer, so I figured I would start there).

Re: [Help-gsl] gsl performance

2013-10-06 Thread Rhys Ulerich
> Rhys, I did try to use views. They do not help because the gsl routines > allocate vectors internally and there is not much that I can do about it... > except for maybe hacking gsl and changing gsl_vector_alloc myself. If from hacking around you happen to restructure the API so that a clean work

Re: [Help-gsl] gsl performance

2013-10-06 Thread onefire
Rhys, I did try to use views. They do not help because the gsl routines allocate vectors internally and there is not much that I can do about it... except for maybe hacking gsl and changing gsl_vector_alloc myself. "The main target of a library like GSL is an inherent consistency of its objects an

Re: [Help-gsl] gsl performance

2013-10-06 Thread Simon Zehnder
I write between the lines. On Oct 6, 2013, at 10:04 PM, onefire wrote: > I am aware of those. I tried the following: > 1) Call gsl_set_error_handler_off() to turn off the error handler > 2) Compile with -DHAVE_INLINE and -DGSL_RANGE_CHECK_OFF > 3) Link to a different cblas (actually I tried op

Re: [Help-gsl] gsl performance

2013-10-06 Thread onefire
I am aware of those. I tried the following: 1) Call gsl_set_error_handler_off() to turn off the error handler 2) Compile with -DHAVE_INLINE and -DGSL_RANGE_CHECK_OFF 3) Link to a different cblas (actually I tried openblas, mkl and gslcblas). However, my most interesting experiment was to modify m

Re: [Help-gsl] gsl performance

2013-10-06 Thread Simon Zehnder
Hi Gilbaerto, congrats on your performance results! A first guess of the worse performance of the gsl library would be exception throwing. The GSL is made for 'users' and this means, that a user has to be informed about the kind of exception he encounters. This can be left out if you have you

[Help-gsl] gsl performance

2013-10-05 Thread onefire
Hi all, I am creating a personal library for my C++ projects and to evaluate its performance I decided to compare it to gsl and, to my surprise, my library is much faster than gsl in many cases. For example, my Nelder-Mead implementation can be 5 or 10 faster for some functions. This has been bot