[SciPy-User] odeint - Excess work done on this call

[Farid Afandiyev]

Hi Nicky,

Correct, flow will stop after pressures are equalized. I would like it to
stop integration at that point, i.e. after pressure equilibrium is reached
odeint should stop calling derivative function.

Thanks,

Farid.

On Tue, Oct 14, 2014 at 5:05 PM, nicky van foreest <vanforeest at gmail.com>
wrote:

> Hi Farid,
>
> Please see below.
>
> [[  4.00000000e+03   2.64000000e-01   2.00000000e+03   3.30000000e-01]
>>  [  3.49242034e+03   2.30499743e-01   2.20303186e+03   3.63500257e-01]
>>  [  3.09580400e+03   2.04323064e-01   2.36167840e+03   3.89676936e-01]
>>  [  2.81015098e+03   1.85469965e-01   2.47593961e+03   4.08530035e-01]
>>  [  2.63546127e+03   1.73940444e-01   2.54581549e+03   4.20059556e-01]
>>  [  2.57173487e+03   1.69734501e-01   2.57130605e+03   4.24265499e-01]
>>  [  2.57142857e+03   1.69714286e-01   2.57142857e+03   4.24285714e-01]
>>
>
> It appears that your P1 - P2 = 0 at this line. It might mean that the
> pressure in both vessels is the same at that moment in time, and then the
> flow should stop, right?  So, what should the integrator tell you after
> this point in time?
>
> Nicky
>
>
>>  [  0.00000000e+00   0.00000000e+00   0.00000000e+00   0.00000000e+00]
>>  [  0.00000000e+00   0.00000000e+00   0.00000000e+00   0.00000000e+00]
>>  [  0.00000000e+00   0.00000000e+00   0.00000000e+00   0.00000000e+00]]
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20141014/09e16202/attachment.html>