Hi, I use the class petsc_nonlinear_solver to solve a nonlinear convection-diffusion type equation in porous media. The built linear system's size is 30000 and its resolution lasts about a hour on a single processor machine (my laptop...). I run for that a hundred time steps and have in average 4 Newton iterations per time step. So a hour seems huge to me with the optimized version of Petsc. I tryed to profile my code with the Petsc options "-log_summary" and "-info". What i get out of it is that the first time step lasts about 15min representing 25% of the total time and i think its not normal. The thing is, these 15min happen at the second call of the petsc routine "VecScatterCreate()" right after the Libmesh output of the first Newton residual.
Any idea? Thanks. Stephane _________________________________________________________________ News, entertainment and everything you care about at Live.com. Get it now! http://www.live.com/getstarted.aspx ------------------------------------------------------------------------------ Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, CA -OSBC tackles the biggest issue in open source: Open Sourcing the Enterprise -Strategies to boost innovation and cut costs with open source participation -Receive a $600 discount off the registration fee with the source code: SFAD http://p.sf.net/sfu/XcvMzF8H _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users