commons-issues mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From "Pascal Parraud (JIRA)" <j...@apache.org>
Subject [jira] Commented: (MATH-388) ODE integrator: different size needed for state vector and tolerance error vector dimension
Date Wed, 21 Jul 2010 13:29:52 GMT

    [ https://issues.apache.org/jira/browse/MATH-388?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=12890725#action_12890725
] 

Pascal Parraud commented on MATH-388:
-------------------------------------

After a few tests, the patch:
*  seems ok for DormandPrince853Integrator,
*  seems ok for  GraggBulirschStoerIntegrator, with 2 fixes :
** in GraggBulirschStoerIntegrator, y0.length changed into scale.length (or mainSetDimension)
at lines 712 & 716 ( integrate method),
** in GraggBulirschStoerStepInterpolator, currentState.length changed into scale.length at
lines 302 & 306 (estimateErrror method),
* needs to be tested for others integrators.

And, waiting for a new architecture, ODEWithJacobians & ParameterizedODE classes should
extend ExtendedFirstOrderDifferentialEquations instead of FirstOrderDifferentialEquations.

> ODE integrator: different size needed for state vector and tolerance error vector dimension
> -------------------------------------------------------------------------------------------
>
>                 Key: MATH-388
>                 URL: https://issues.apache.org/jira/browse/MATH-388
>             Project: Commons Math
>          Issue Type: Bug
>            Reporter: Pascal Parraud
>            Assignee: Luc Maisonobe
>             Fix For: 2.2
>
>         Attachments: math-388.patch
>
>
> The user should be allowed to chose a tolerance vector dimension different from the state
vector dimension.
> For example, using the FirstOrderIntegratorWithJacibians, we don't want to set some tolerance
error for the jacobian part in order to not interfere with the main state vector integration.
> This is really a problem with the dimension of the tolerance vector, there is no work-around
assigning some particular value to the tolerance vector.

-- 
This message is automatically generated by JIRA.
-
You can reply to this email to add a comment to the issue online.


Mime
View raw message