Re: [ff3d-users] Strange error during execution

[Stephane Del Pino]

Dear Emilio,

Please, reply to the list.

You did not join the '3layer.pov' file. I created one and run the example.
I did not noticed any memory leak with the current version of ff3d. Which 
version are you using: 'ff3d -v'
If you are using the cvs-release, I have no explanation and I will try to 
investigate a bit more...

Best regards,
Stephane.

Le mardi 19 décembre 2006 09:08, Emilio Melero García a écrit :
> Dear Stephane, here is the text of the file:
>
> vector a = (0,0,0);
> vector b = (9,0.5,0.5);
> vector n = (200,13,13);
> mesh M = structured(n,a,b);
>
> double D = 1.1E-5;
> double dt = 12;
> double numiter = 43200;
> double numpoints = 1024;
>
> scene S = pov("3layer.pov");
> domain Omega = domain(S,inside(<1,0,0>));
>
> function conc0=2*one(x<=3.8)+0.5*one(3.8<x && x<=5.3);
>
> for (double i = 1;i<= numiter; i++)
> {
>   solve(conc) in Omega by M
>     {
>       test(v)
>     int(conc*v+dt*D*grad(conc)*grad(v))
>     -int(conc0*v)=0;
>       //conc=0 on M;
>     };
>   conc0 = conc;
>   cout << "Iteration number: ";
>   cout << i <<"\n";
> }
>
> double xpos=0;
> for (double m=1;m<=numpoints;m++)
>   {
>     cout << conc(xpos,0.25,0.25) <<"\n";
>     xpos=xpos+9/numpoints;
>   }
>
> Thanks for  your reply
>
> Best regards
>
> Emilio