×
INTELLIGENT WORK FORUMS
FOR COMPUTER PROFESSIONALS

Contact US

Log In

Come Join Us!

Are you a
Computer / IT professional?
Join Tek-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!

*Tek-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Hello I demand your help . I'm go

Hello I demand your help . I'm go

Hello I demand your help . I'm go

(OP)
Hello
I demand your help .
I'm going to study the accretion-ejection disks using the DVODE-f90 code.
when I compile my example. the calcul stops at T=3,21.

the above message, R1 = 0.3215402479903D+01
At t = 0.3215D+01 y = 0.321568D+01 0.158443D+14 0.193170D+39



Error halt: ISTATE = -1

Can you help me how to integrate my equations system for large T (INFINITY)
you find in the attached files my source file and the dvode_f90 subroutine.

RE: Hello I demand your help . I'm go

Hello hasnaoui,

Unfortunately, I can't help you with that. Solving such systems requires some experience that I no longer have. I did something like that about 20 years ago - but haven't since.

You didn't even state what the system of 11 differential equations actually looks like. I mean, in order to be able to solve a system of ODEs, I must first have an idea of ​​what kind of system it is, i.e. I must first see it written somewhere - in an article or in a book or somewhere else.

You seem to be trying to solve the system as non-stiff, but most of the time the systems of many ODEs are stiff.
Therefore, if I were you, I would try to use the stiff options - described in the source of dvode_f90_m.f90 from line 213: https://github.com/valkenburg/vd2020/blob/master/s...
Also I would try to explore the stiffoptions.f90 program, to see what stiff options they use: https://github.com/valkenburg/vd2020/blob/master/s...

RE: Hello I demand your help . I'm go

You are not going to get any help if you demand it.

RE: Hello I demand your help . I'm go

(OP)
Dear Mikrom (programmer)
Thank you for your responses

Exactly , the 11 differential equations describe the magnetized accretion-ejection disks arround young stellar objects.
These equations are developed by Jonathan Ferreira et al.
I asked him how to resolve these equations for long cordinate space x=z/h. unfortunately, they dont response me.
best regards

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Tek-Tips Forums free from inappropriate posts.
The Tek-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Tek-Tips forums is a member-only feature.

Click Here to join Tek-Tips and talk with other members! Already a Member? Login


Close Box

Join Tek-Tips® Today!

Join your peers on the Internet's largest technical computer professional community.
It's easy to join and it's free.

Here's Why Members Love Tek-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close