From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 21483 invoked by alias); 2 Jan 2008 13:40:11 -0000 Received: (qmail 21454 invoked by uid 22791); 2 Jan 2008 13:40:10 -0000 X-Spam-Check-By: sourceware.org Received: from smtpout3.cns.ohiou.edu (HELO smtpout3.cns.ohiou.edu) (132.235.51.106) by sourceware.org (qpsmtpd/0.31) with ESMTP; Wed, 02 Jan 2008 13:33:32 +0000 X-IronPort-AV: E=Sophos;i="4.24,234,1196658000"; d="asc'?scan'208";a="350502147" Received: from oak4a.cats.ohiou.edu (HELO oak.cats.ohiou.edu) ([132.235.8.203]) by smtpout3.cns.ohiou.edu with ESMTP; 02 Jan 2008 08:33:02 -0500 Received: from 132.235.227.198 by pm4 (PureMessage); Wed Jan 2 08:32:35 2008 Received: from [192.168.0.2] (crab-lab.zool.ohiou.edu [132.235.227.198]) (authenticated bits=0) by oak.cats.ohiou.edu (8.13.1/8.13.1) with ESMTP id m02DWY9a2523247 (version=TLSv1/SSLv3 cipher=RC4-MD5 bits=128 verify=NOT); Wed, 2 Jan 2008 08:32:34 -0500 (EST) Subject: Re: Bug report and fix for adaptive step size control in ODE solving From: "Kevin H. Hobbs" To: Frank Reininghaus Cc: gsl-discuss@sourceware.org In-Reply-To: <4772D239.9070206@googlemail.com> References: <4772D239.9070206@googlemail.com> Content-Type: multipart/signed; micalg=pgp-sha1; protocol="application/pgp-signature"; boundary="=-KLRUWL68gHWJAUsiEexx" Date: Wed, 02 Jan 2008 13:40:00 -0000 Message-Id: <1199280754.30472.14.camel@gargon.hooperlab> Mime-Version: 1.0 X-Mailer: Evolution 2.10.3 (2.10.3-4.fc7) X-PMX-Information: http://technology.ohio.edu/email/filtering/ X-PMX-Spam: Gauge=IIIIIII, Probability=7%, Report='__CT 0, __CTYPE_HAS_BOUNDARY 0, __CTYPE_MULTIPART 0, __HAS_MSGID 0, __HAS_X_MAILER 0, __MIME_VERSION 0, __SANE_MSGID 0' 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: 2008-q1/txt/msg00000.txt.bz2 --=-KLRUWL68gHWJAUsiEexx Content-Type: text/plain Content-Transfer-Encoding: quoted-printable Content-length: 819 On Wed, 2007-12-26 at 23:14 +0100, Frank Reininghaus wrote: > Hi, >=20 > I've found some problems with adaptive step size control which can > lead to an infinite loop under special circumstances.=20 > In testcase1.c, the function rhs (t) yields 1 for t>=3D1.0 and 0 > otherwise.=20 I may have encountered this when using step functions for current injection into neuron models. I'd like to be included in the replies.=20 My workaround was to adjust the step size so that the step ends at the time of the current injection change. I had to tell the model to change the current through the comm pointer so that the end point of the step before the change and the start point of the step after the change have the same current level as the rest of the step. This workaround is probably too complex for beginners. --=-KLRUWL68gHWJAUsiEexx Content-Type: application/pgp-signature; name=signature.asc Content-Description: This is a digitally signed message part Content-length: 189 -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.7 (GNU/Linux) iD8DBQBHe5JuqBtEuW+gRPERAvfWAJ9W1gG4wgPWHfHstQkGTouNWxu57wCfc0Ku AWKf1RscpOCkr0VUWUs02iM= =XvdZ -----END PGP SIGNATURE----- --=-KLRUWL68gHWJAUsiEexx--