Re: [SOSlib-devel] Problem with senstivities

[Rainer M. <ra...@tb...>]

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.

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
>