Mailing-List: contact dev-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: "Commons Developers List" <dev@commons.apache.org>
Received-SPF: pass (athena.apache.org: domain of phil.steitz@gmail.com
 designates 209.85.160.43 as permitted sender)
DomainKey-Signature: a=rsa-sha1; c=nofws;
        d=gmail.com; s=gamma;
        h=message-id:date:from:user-agent:mime-version:to:subject:references
         :in-reply-to:content-type:content-transfer-encoding;
        b=x+E07dL4MnwSoCN/f7I74UcpiMHxK5qhOeMxLP+SYDTHXtiPIeevZyU3J/kLCEpxLB
         kwIyKLdGZlTuctG3lEQ99PjMlK4ePMqVH/bo6pBQPfHyNnkl/0TvyQAMKF6R+VG4Ds/8
         UGf+8+yym/s+NbZGSasjA6Rkh4DcN1LZRM4Gc=
Message-ID: <4DA47985.3070903@gmail.com>
Date: Tue, 12 Apr 2011 09:10:45 -0700
From: Phil Steitz <phil.steitz@gmail.com>
User-Agent: Mozilla/5.0 (Macintosh; U; Intel Mac OS X 10.6; en-US;
 rv:1.9.2.15) Gecko/20110303 Lightning/1.0b2 Thunderbird/3.1.9
MIME-Version: 1.0
To: Commons Developers List <dev@commons.apache.org>
Subject: Re: [math] resetting relative and absolute tolerances in ODE adaptive
 step size integrators
References: 
 <379411337.1755041302622016003.JavaMail.root@spooler6-g27.priv.proxad.net>
In-Reply-To: 
 <379411337.1755041302622016003.JavaMail.root@spooler6-g27.priv.proxad.net>
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

On 4/12/11 8:26 AM, luc.maisonobe@free.fr wrote:
> Hi Phil,
>
> ----- "Phil Steitz" <phil.steitz@gmail.com> a écrit :
>
>> On 4/12/11 1:51 AM, luc.maisonobe@free.fr wrote:
>>> Hi all,
>>>
>>> I have hit a limitation of the current implementation of ODE
>> integrators with adaptive step size. For now, the tolerances that are
>> used to adjust the step size are specified only at construction time
>> and cannot be changed afterwards. However, these tolerances are highly
>> problem-dependent and in fact the dimension of the problem (which is
>> related to the dimension of the vectorial version of the tolerances)
>> is specified only at integration time, not at construction time.
>>> So I consider adding at the top-level hierarchy (abstract class
>> AdaptiveStepsizeIntegrator) a few setters to allow users to change
>> these tolerances after the integrator has been built. It seems the
>> integrators by themselves were not documented as immutable (I first
>> thought they were), so this change is probably harmless.
>>> I am going to open a Jira issue for this.
>>>
>>> Any thoughts ?
>> One natural thing to consider is to provide the tolerances at
>> integration time - i.e., to either add the vectors explicitly as
>> arguments to the integrate method or to create an ODEProblem or
>> better-named class that encapsulates the DE, tolerances, initial
>> values and time parameters.
> No, it would break the interface that is shared with fixed step size integrators.
> Preserving the compatibility with both types of integrators is very important. This
> is the same reason why step size is not set at integration time for fixed step integrator.

Got it.  What about putting the specialization in the Problem class
rather than in the integrators?  So the fixed step integrator gets a
specialized integration problem rather than exposing problem
parameters itself?  Would something like that be possible? 

Phil
> Luc
>
>> Phil
>>> Luc
>>>
>>>
>> ---------------------------------------------------------------------
>>> To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
>>> For additional commands, e-mail: dev-help@commons.apache.org
>>>
>>>
>>
>> ---------------------------------------------------------------------
>> To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
>> For additional commands, e-mail: dev-help@commons.apache.org
> ---------------------------------------------------------------------
> To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
> For additional commands, e-mail: dev-help@commons.apache.org
>
>


---------------------------------------------------------------------
To unsubscribe, e-mail: dev-unsubscribe@commons.apache.org
For additional commands, e-mail: dev-help@commons.apache.org