From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 16338 invoked by alias); 4 Jun 2009 17:40:51 -0000 Received: (qmail 16320 invoked by uid 22791); 4 Jun 2009 17:40:47 -0000 X-SWARE-Spam-Status: No, hits=-2.0 required=5.0 tests=BAYES_00,SARE_MSGID_LONG40,SPF_PASS X-Spam-Check-By: sourceware.org Received: from mail-ew0-f211.google.com (HELO mail-ew0-f211.google.com) (209.85.219.211) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Thu, 04 Jun 2009 17:40:39 +0000 Received: by ewy7 with SMTP id 7so1301994ewy.24 for ; Thu, 04 Jun 2009 10:40:35 -0700 (PDT) MIME-Version: 1.0 Received: by 10.216.18.205 with SMTP id l55mr889235wel.194.1244137235349; Thu, 04 Jun 2009 10:40:35 -0700 (PDT) In-Reply-To: <4A28028B.8020403@iki.fi> References: <4a00655d0906030731l8632674o9abb33dca522e689@mail.gmail.com> <4A28028B.8020403@iki.fi> From: Rhys Ulerich Date: Thu, 04 Jun 2009 17:40:00 -0000 Message-ID: <4a00655d0906041040y135755cdlc877ba570c4308da@mail.gmail.com> Subject: Re: Implementing hybrid explicit/implicit ODE stepper To: Tuomo Keskitalo Cc: Brian Gough , gsl-discuss@sourceware.org Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable 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: 2009-q2/txt/msg00013.txt.bz2 > Sorry, your description of the method is not clear to me. Can you send me > the reference as e.g. pdf? Sure, scanned images (sorry no PDF) of the three page appendix detailing the method are at http://users.ices.utexas.edu/~rhys/smr_appendix/page1.jpg http://users.ices.utexas.edu/~rhys/smr_appendix/page2.jpg http://users.ices.utexas.edu/~rhys/smr_appendix/page3.jpg Though the appendix talks about A. Wray's submission to J. Comput. Phys., that paper was never published. This integrator tends to get used in DNS of Navier-Stokes since the memory savings become significant for very large fields. - Rhys On Thu, Jun 4, 2009 at 12:21 PM, Tuomo Keskitalo w= rote: > On 06/04/2009 10:18 AM, Brian Gough wrote: > >> At Wed, 3 Jun 2009 09:31:40 -0500, >> Rhys Ulerich wrote: >>> >>> What's the "right way" to shoehorn this into the gsl_odeiv_system >>> struc? =A0L(u) as gsl_odeiv_system.function with N(u) passed within >>> gsl_odeiv_system.params ? =A0The reverse? =A0Neither? >> >> That's a good question, I'm not sure - maybe Tuomo has some >> suggestions? > > > -- > Tuomo.Keskitalo@iki.fi > http://iki.fi/tuomo.keskitalo >