Dear Stephane, thank you so much for your help and patience. I will certainly increase the time step. It is a very good solution even with the memory leak since i do not need to go foward in time so much (no more than 13 days, which are only 936 iterations with dt=1200).
By the way, and i understand it is not the place to ask these, and that you might be very busy to answer, the stability condition that you mention and the residual of the conjugate gradient is due to the problem i am trying to solve or is inherent to the PDE solver you have implemented?, and if it is problem dependent, could you hint me on how did you evaluate the residual and the stability?. And finally, i am trying to make my way into all this world of FEM problems and solvers, but i still have a notion that i might be resorted to implement my own very specific FEM code, any good book you would recommend to learn FEM methods and implementations? Please, as i said i understand this might be not the place to ask these things but i have felt that people in these list read the e-mails and answer when they can. So please, consider this a question to the people in the list. Thank you very much from a newbie. 2006/12/20, Stephane Del Pino <[EMAIL PROTECTED]>:
Dear Emilio, With the POV data you sent, it seems that there is a small memory leak... I will investigate and let you know. However I think you should use a bigger time-step: the scheme is unconditionaly stable and here, the initial conjugate gradient residu is about 10^(-12)... I think using the timestep dt=1200 would be better. Anymore, this would require only 500 iterations to converge... Best regards, Stephane. _______________________________________________ ff3d-users mailing list [email protected] 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 ###############################################
_______________________________________________ ff3d-users mailing list [email protected] http://lists.nongnu.org/mailman/listinfo/ff3d-users
