Search results for query: *

Dear Programmer Thank you very much for your help The same results attached to your response are obtained using the dvode code. My objective is to integrate the stiff equations for large time (infinity). I don't know if the Lsode allows to resolve the stiff equations for large time Again, i...

Hello I need your help can you give me an example in which we solve a stiff problems with a sparce jacobian best wiches

Hello when I running the example2.f90 an error appers like: 88 | CALL DVODE_F90(FEX,NEQ,Y,T,TOUT,ITASK,ISTATE,OPTIONS, G_FCN=GEX) | 1 Error: There is no specific subroutine for the generic ‘dvode_f90’ at (1)...

Help me pleeease when I running my example using the dvode subroutine dvode.f90 and the options OPTIONS = SET_OPTS(BANDED_J=.TRUE.,RELERR=RTOL,ABSERR_vector=ATOL, & LOWER_BANDWIDTH=ML,UPPER_BANDWIDTH=MU) some errors appear like: 213 | LOWER_BANDWIDTH=ML,UPPER_BANDWIDTH=MU)...

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...

Thank you very much for your help. but when I coming to run the stiffoptions.90 using the dvode_f90 subroutine, some errors appear like: 1982 | DVODE_F90(DERIVS,NEQ,Y,T,TOUT,ITASK,ISTATE,OPTIONS,J_FCN=JACD) | 1 Error: There is...

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 =...

Hello When I try compiling the example1.f90 I get: $f955 dvode_f90_m.f90 example1.f90 -o test CALL DVODE_F90(FEX,NEQ,Y,T,TOUT,ITASK,ISTATE,OPTIONS,J_FCN=JEX) | 1 Error: There is no specific subroutine for the generic...

Thank you for this infromation. Ok, when I chose MF=22 some errors appear during the execution like DVODE-- At current T (=R1), MXSTEP (=I1) steps taken on this call before reaching TOUT In above message, I1 =...

Hello I would ask you if the dvode fortran code can generate automatically the jacobien matrix or we mut added it best regards

Deear Collegue I'm working on the stationary accretion-ejection disks using the fortran code DVODE. The calcul stops for short length from the disk surface. Can you help me how to integrate the equations to infinity using the devodesubroutine. you find the fortran code in the attached file...