Re: [SOSlib-devel] Bug in sensitivity analysis?

[Witold W. <wew...@gm...>]

Hi Rainer,

The point I wanted to make is that I am wondering why depending if I use th=
e
function globalizeParameter or not prior to integration I am getting
different results.

cheers
Eryk


On 5/10/06, Rainer Machne <ra...@tb...> wrote:
>
>
>
> > If I do not call it the output reads:
> >
> > CVode-- At t =3D 0 and h =3D 1.12418e-06, the corrector
> > convergence failed repeatedly or with |h| =3D hmin.
> >
> > 1.12418e-06  0 0 0 0 0 0 0 0
> > 1.12418e-06  10  0 0 0 0
>
>
> The first two lines mean that cvodes has failed for integration with
> sensitivity analysis. Thus results and sensitivities were not calculated.
>
> IntegratorInstance_dumpY/PSensitivities
>
> print first the current time, then the current value of the
> variable/parameter, and then the sensitivities to all constants of the
> model.
>
> I have never noticed however that upon failure the current time seems to
> be the h, which is cvodes' internal timestep. So the sensitivities for
> this time h are wrong and the output should look different. Maybe upon
> failure cvodeData the time should be reset to the last valid time.
>
>
> Rainer
>
> On Wed, 10 May 2006, Witold Wolski wrote:
>
> > Hello,
> >
> > Depending if I apply globalise parameter or not I am getting different
> > outputs in the sensitivity analysis.
> >
> > Given the MAPK.xml
> >
> > Namely the output of the functions (see attachment)
> >
> > IntegratorInstance_dumpPSensitivities(integratorInstance,p1);
> > IntegratorInstance_dumpYSensitivities(integratorInstance,y1);
> >
> > differs depedning if
> >
> > globalizeParameter(model,"V2","J1");
> >
> > was called.
> >
> > If I do not call it the output reads:
> >
> > CVode-- At t =3D 0 and h =3D 1.12418e-06, the corrector
> > convergence failed repeatedly or with |h| =3D hmin.
> >
> > 1.12418e-06  0 0 0 0 0 0 0 0
> > 1.12418e-06  10  0 0 0 0
> >
> > I am working on a very recent cvs checkout.
> >
> > cheers
> >
> > Eryk
> >
>