Re: [SOSlib-devel] Problem with senstivities

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

On 8/16/06, Rainer Machne <ra...@tb...> wrote:
> Hi Eryk,
>
> this seems to be a bug! Again, I am not in Vienna, have limited access to
> computers and can't look into the problem. It looks as if the _reset
> function doesn't reset the sensitivities to 0. I guess, this problem arose
> because events should not reset sensitivities but keep the values from
> the last time step befor the solver structures where re-initialized
> by an event.


I was looking at the reset function and it seems to be fine:

the reset function is calling I think the set function and that one is
calling the
CvodeData_initialize where the sensitivity matrix in the cvodeData_t
is set to zero.

Furthermore, what makes me think that this is not the cause of the
error is that the senstivities not only change the magnitude
(increase) but also swap the signs.
I therefor would think that the reason for the error is somewhere else.
It should be relatively easy to reproduce the error just by running
the solver in a while loop as I described and printing out the
sensitivities.


cheers




>
> Maybe you have time to look into it yourself. The function
> IntegratorInstance_reset probably calls things like CvodeData_initialize
> at it most somewhere also go into initializing the sensitivities. Both
> data->sensitivities and results->sensitivities should be set to 0 by the
> _reset function.
>
> If you don't find it yourself and maybe post the problem/solution to this
> list or even update the CVS with a solution, then I will look into it next
> week.
>
> Rainer
>
> On Tue, 15 Aug 2006, W Eryk Wolski wrote:
>
> > Hi,
> >
> > I am running the solver in a loop.
> >
> > the loop includes:
> > {
> > change of the parameters
> > running the solver
> > and a call to IntegratorInstance_reset(integratorInstance);
> > }
> >
> > What I observe and call a problem is that every new simulator run in
> > the loop generates increasing sensitivity values for the same
> > parameter values and model.
> >
> > The following output are the date stored in:
> >
> > integratorInstance->data->sensitivity
> > entries
> > integratorInstance->data->nsens, integratorInstance->data->neq
> >
> > first iteration
> >
> > parameters:  par[0] 2.000000  par[1] 8.000000  par[2] 10.000000
> > par[3] 0.100000  par[4] 2.000000
> >
> > 24.33599      -12.44311       -1.60336        0.28326 26.31235        -0.12182        0.00000 0.00000
> > -24.33599     12.44311        1.60336 -0.28326        -26.31235       0.12182 0.00000 11.93225
> > 11.93225      -4.09705        -0.54112        0.08923 8.68924 -0.04727        0.00000 -7.76276
> > -7.76276      2.05021 0.27103 -0.04458        -4.34843        0.02380 0.00000 -4.16949
> > -4.16949      2.04684 0.27009 -0.04465        -4.34081        0.02346 0.00000 -0.26712
> > -0.26712      -0.45221        -0.06041        0.00965 0.95925 -0.00565        0.00000 -1.57271
> >
> >
> > while the 10 iteration:
> >
> > -1394.82594   399.44848       59.25448        -19.98846       -2225.46440     3.11907 0.00000 0.00000
> > 1394.82594    -399.44848      -59.25448       19.98846        2225.46440      -3.11907        0.00000 -21877.33372
> > -21877.33372  6547.28853      971.83061       -324.87728      -36996.03248    52.60350        0.00000 -2206.49496
> > -2206.49496   666.20014       98.97768        -33.06429       -3791.30589     5.39430 0.00000 24083.82868
> > 24083.82868   -7213.48867     -1070.80829     357.94157       40787.33837     -57.99780       0.00000 -21953.82949
> > -21953.82949  6578.30322      976.54971       -326.41598      -37205.87555    52.90915        0.00000 12235.00898
> >
> >
> >
> > I hope that the problem statement is clear but I am happy to answer
> > more questions.
> > I am looking forward to suggestions.
> >
> > Eryk
> >
> > -------------------------------------------------------------------------
> > Using Tomcat but need to do more? Need to support web services, security?
> > Get stuff done quickly with pre-integrated technology to make your job easier
> > Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
> > http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
> > _______________________________________________
> > sbmlsolver-devel mailing list
> > sbm...@li...
> > https://lists.sourceforge.net/lists/listinfo/sbmlsolver-devel
> >
>