ff3d-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [ff3d-users] Strange error during execution


From: Emilio Melero García
Subject: Re: [ff3d-users] Strange error during execution
Date: Wed, 20 Dec 2006 09:05:39 +0100

Dear Stephane,

I am very sorry, i did not realize that i was replying directly to your mail box. I will be more careful now. :-) Sorry also about the missing pov. I attach it now:

merge{
cylinder {
<0, 0.25, 0.25>, <5.3, 0.25, 0.25>, 0.205
}
cylinder {
<5.3, 0.25, 0.25>, <9, 0.25, 0.25>, 0.05
}
pigment { color rgb <1,0,0>}
}


The version is 1.0pre7 because is the most recent one availabe for MAC OS X. I will not be able to compile the CVS version until a couple of weeks. I am behind a too restrictive firewall and i need several  packages to compile ff3d. Maybe is just somekind of bug inside a too old version??

Thank  you so much for your answers

Emilio



2006/12/19, Stephane Del Pino <address@hidden>:
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


_______________________________________________
ff3d-users mailing list
address@hidden
http://lists.nongnu.org/mailman/listinfo/ff3d-users



--
###############################################
Dr. Emilio Melero García
Laboratorio de Estudios Cristalográficos, IACT
Instituto de parasitología Lopez Neyra,
P.T. Ciencias de la Salud, E-18100 Armilla,
Granada, España

Correo elec. :emilio_arroba_lec_._csic_._es
Tfno: 958181621, lab 107
###############################################
reply via email to

[Prev in Thread] Current Thread [Next in Thread]