Thread Subject: how to solve my error

I am solving a PDE by ode 15s.
I am getting error massage (below).
Warning: Failure at t=1.173441e-001. Unable to meet integration
tolerances without reducing the step size below the smallest value
allowed (2.220446e-016) at time t.
> In ode15s at 819
  In mypde_3_main1_with_PS_OP_ES at 96
??? Index exceeds matrix dimensions.

Error in ==> mypde_3_main1_with_PS_OP_ES at 106
    N(it,i,j)=y(it,(i-1)*ne+j);
 
Can you give me the idea that where is the problem?
Thanks

On 24 Okt., 15:19, "Sunipa Som" <sunipa_...@yahoo.com> wrote:
> I am solving a PDE by ode 15s.
> I am getting error massage (below).
> Warning: Failure at t=1.173441e-001.  Unable to meet integration
> tolerances without reducing the step size below the smallest value
> allowed (2.220446e-016) at time t.> In ode15s at 819
>
>   In mypde_3_main1_with_PS_OP_ES at 96
> ??? Index exceeds matrix dimensions.
>
> Error in ==> mypde_3_main1_with_PS_OP_ES at 106
>     N(it,i,j)=y(it,(i-1)*ne+j);
>
> Can you give me the idea that where is the problem?
> Thanks

1. Check your ODEs since the integrator has problems to integrate them
right at the beginning.
2. You seem to exceed the array bounds of the matrix y with the index
it and/or (i-1)*ne+j
   or the array bounds of the matrix N with the index it, i and/or j.

Not very helpful, I guess, but nothing more can be said unless you
provide more information.

Best wishes
Torsten.

Torsten <Torsten.Hennig@umsicht.fraunhofer.de> wrote in message <236b8371-8aec-4b5e-bb4d-1fb17870e2f6@u10g2000prl.googlegroups.com>...
> On 24 Okt., 15:19, "Sunipa Som" <sunipa_...@yahoo.com> wrote:
> > I am solving a PDE by ode 15s.
> > I am getting error massage (below).
> > Warning: Failure at t=1.173441e-001.  Unable to meet integration
> > tolerances without reducing the step size below the smallest value
> > allowed (2.220446e-016) at time t.> In ode15s at 819
> >
> >   In mypde_3_main1_with_PS_OP_ES at 96
> > ??? Index exceeds matrix dimensions.
> >
> > Error in ==> mypde_3_main1_with_PS_OP_ES at 106
> >     N(it,i,j)=y(it,(i-1)*ne+j);
> >
> > Can you give me the idea that where is the problem?
> > Thanks
>
> 1. Check your ODEs since the integrator has problems to integrate them
> right at the beginning.
> 2. You seem to exceed the array bounds of the matrix y with the index
> it and/or (i-1)*ne+j
> or the array bounds of the matrix N with the index it, i and/or j.
>
> Not very helpful, I guess, but nothing more can be said unless you
> provide more information.
>
> Best wishes
> Torsten.
>
Hello,
My programme is working upto t=1e-001. You see from warning it is failure at t=1.173441e-001.
Another point is, in my programming N is the particle number. In the initial distribution there is one parameter. By changing this parameter I can change total particle number.
When particle number is 10^7 or below than this, it is working. But over 10^7 it is not working.
I can not understand why it is working for 10^7 number of particles but not working over 10^7 number of particles.
Thank you.