From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 11646 invoked by alias); 6 Jun 2007 23:54:14 -0000 Received: (qmail 11638 invoked by uid 22791); 6 Jun 2007 23:54:13 -0000 X-Spam-Check-By: sourceware.org Received: from nz-out-0506.google.com (HELO nz-out-0506.google.com) (64.233.162.230) by sourceware.org (qpsmtpd/0.31) with ESMTP; Wed, 06 Jun 2007 23:54:12 +0000 Received: by nz-out-0506.google.com with SMTP id l1so291397nzf for ; Wed, 06 Jun 2007 16:54:10 -0700 (PDT) Received: by 10.115.60.1 with SMTP id n1mr958160wak.1181174049875; Wed, 06 Jun 2007 16:54:09 -0700 (PDT) Received: by 10.114.26.2 with HTTP; Wed, 6 Jun 2007 16:54:09 -0700 (PDT) Message-ID: Date: Wed, 06 Jun 2007 23:54:00 -0000 From: Michael To: help-gsl@gnu.org, gsl-discuss@sources.redhat.com Subject: Will I get speedup by providing Jacobian to the GSL ODE45 solver? MIME-Version: 1.0 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Content-Disposition: inline Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2007-q2/txt/msg00036.txt.bz2 Will I get speedup by providing Jacobian to the GSL ODE45 solver? Hi all, I am a newbie to GSL ODE solvers. I have a critical need for speed. I am hunting for ways to speed up. My questions are: 1. If after a comparison of GSL ODE solvers, I conclude that ODE45 is sufficient for my need in terms of accuracy. (And in fact I turned to GSL ODE solver from Matlab, and I have already done a comparison analysis using all the ODE solvers in Matlab. ) And I remember I don't provide Jacobian to ODE45 solver. And now will I get speed-up at the fixed accuracy level by providing the Jacobian to the GSL ODE45 solver? 2. Is it possible that after providing the Jacobian, other solvers use the Jacobian more efficiently and start to outform ODE45 in terms of speed or accuracy? I am still experimenting these things, but I would like to hear your comments too... Thanks for your help!