Forums

F.A.Q. F.A.Q.    Register Register    Login Login    Home Home
Search Search
SBML Discussions » sbml-discuss » Finding Stable Points in a system of ODEs
Show: Today's Posts  :: Message Navigator
| Subscribe to topic 
Return to the default flat view Create a new topic Submit Reply
AuthorTopic
Stefan.Hoops


Posts: 170
Registered:
December 2006
Re: Finding Stable Points in a system of ODEs 07 May '12 08:43 Go to previous message

Hello Ayesha,

The COPASI error message indicates that COPASI is trying its last
attempt in finding a steady-state by backwards integration in time. This
very often fails with the error message you see.

Since Newton and forward integration did not find a steady-state please
have a look at your time course to get an idea of whether your system
actually has a steady-state. If on of your variables keeps increasing
your system does not have a steady-state and thus COPASI will not be
able to find it. You may also reduce resolution in steady-state task
which defaults to 1e-9 to e.g. 1e-6.

I could help you with details of your model if you send it to me.

Thanks,
Stefan


On 5/7/12 10:34 AM, ayesha wrote:
>
> Thank you Herbert, Luca and Stephan for your suggestions.
>
> I have already tried COPASI, but it could not solve the system (it think it is too stiff).
>
>
>> EXCEPTION 2012-05-07 09:39:27<
> CTrajectoryMethod (6): Deterministic integration failed. LSODA reported:
> DLSODA- At T(=R1) and step size H(=R2), the error
> test failed repeatedly or with ABS(H) = HMIN
> In above message, R1 = '85.1645', R2 = '1.61925e-09'
>
>
> Please see result for indications of numerical instability.
>> EXCEPTION 2012-05-07 09:39:28<
> CTrajectoryMethod (25): Invalid state at time '-0.010000'.
>> WARNING 2012-05-07 09:39:29<
> CEigen (2): Failed to compute Eigen value with index '70'.
>> WARNING 2012-05-07 09:39:30<
> CEigen (2): Failed to compute Eigen value with index '69'.
>
>
> I don't know how far I should go in integration. I have set the current iteration limit to be very high though (1e+06).
>
> Therefore, just knowing if the system has any stable point or not would be useful.
>
> Thanks.
> Ayesha
> ____________________________________________________________
> To manage your sbml-discuss list subscription, visit
> https://utils.its.caltech.edu/mailman/listinfo/sbml-discuss
>
> For a web interface to the sbml-discuss mailing list, visit
> http://sbml.org/Forums/
>
> For questions or feedback about the sbml-discuss list,
> contact sbml-team@caltech.edu


--
Stefan Hoops, Ph.D.
Senior Project Associate
Virginia Bioinformatics Institute - 0477
Virginia Tech
Bioinformatics Facility II
Blacksburg, Va 24061, USA

Phone: (540) 231-1799
Fax: (540) 231-2606
Email: shoops@vbi.vt.edu
____________________________________________________________
To manage your sbml-discuss list subscription, visit
https://utils.its.caltech.edu/mailman/listinfo/sbml-discuss

For a web interface to the sbml-discuss mailing list, visit
http://sbml.org/Forums/

For questions or feedback about the sbml-discuss list,
contact sbml-team@caltech.edu

      

SubjectPosterDate
Read Message   Finding Stable Points in a system of ODEs ayesha04 May '12 13:19
Read Message   Re: Finding Stable Points in a system of ODEs Stefan.Hoops04 May '12 14:07
Read Message   Re: Finding Stable Points in a system of ODEs eroeurbano04 May '12 18:01
Read Message   Re: Finding Stable Points in a system of ODEs Herbert M Sauro04 May '12 19:54
Read Message   Re: Finding Stable Points in a system of ODEs ayesha07 May '12 06:34
Read Message   Re: Finding Stable Points in a system of ODEs  Stefan.Hoops07 May '12 08:43
Previous Topic:Release of libSBML-5.5.0
Next Topic:Finding Stable Points in a system of ODEs
Go to forum:
-=] Back to Top [=-

Powered by FUDforum. (Copyright Advanced Internet Designs Inc.)

Please use our issue tracking system for any questions or suggestions about this website.