From: Benjamin S. Kirk <benkirk@cf...>  20040507 20:47:25

Sorry I have been quiet for the past week, I had surgery and only today feel well enough to type... I did notice the sane problem. It seems that the Laspack iterative solvers were aborting after 0 iterations in the subsequent solves of ex13. I'm still not sure *why* this is. However, I changed the Laspack interface to call the iterative solvers with a 0 initial guess. That's why the CVS version does not have this issue. I was hoping to release 0.4.3 before I went in for surgery, but it didn't happen. It should be availible soon. Ben Rob van Tol wrote: > Dear Libmesh users, > > I have compiled Libmesh, and tried several different options. The > program runs fine in all cases and the output of all examples is > reproduced except for example 13. Whatever I try, I do not seem to be > able to reproduce the output given on the website. There is only a > difference between u and u_old in the first step. In the following steps > the difference is always 0. If I look at the results with GMV, all > output files are the same. This indicates that after the first time step > a stationary state has been reached, which seems impossible to me. > > Is this a known problem ? Is there something wrong with the code in > example 13, or do I have an error in my installation ? I have included > the output of the program below. > > Your help will be appreciated, > > Kind regards, > > Rob van Tol. > > *************************************************************** > Running Example ./ex13 > *************************************************************** > > Mesh Information: > mesh_dimension()=2 > spatial_dimension()=3 > n_nodes()=1681 > n_elem()=400 > n_local_elem()=400 > n_active_elem()=400 > n_subdomains()=1 > n_processors()=1 > processor_id()=0 > > EquationSystems > n_systems()=1 > System "NavierStokes" > Type "TransientImplicit" > Variables="u" "v" "p" > Finite Element Types="0" "0" "0" > Approximation Orders="2" "2" "1" > n_dofs()=3803 > n_local_dofs()=3803 > n_constrained_dofs()=0 > n_vectors()=1 > n_parameters()=2 > Parameters: > "linear solver maximum iterations"=250 > "linear solver tolerance"=0.001 > > Solving time step 0, time = 0.005 > Nonlinear convergence: u  u_old = 224.402 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 1 > Solving time step 1, time = 0.01 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 2, time = 0.015 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 3, time = 0.02 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 4, time = 0.025 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 5, time = 0.03 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 6, time = 0.035 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 7, time = 0.04 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 8, time = 0.045 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 9, time = 0.05 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 10, time = 0.055 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 11, time = 0.06 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 12, time = 0.065 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 13, time = 0.07 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > Solving time step 14, time = 0.075 > Nonlinear convergence: u  u_old = 0 > Nonlinear solver converged at step 0 > >  >  Time: Wed May 5 15:55:39 2004 >  OS: Linux >  HostName: pc02023 >  OS Release 2.4.208 >  OS Version: #1 Thu Mar 13 17:54:28 EST 2003 >  Machine: i686 >  Username: root >  >  >  Example 13 Performance: Alive time=8.19504, Active time=5.4663 >  >  Event nCalls Total Avg Percent > of  >  Time Time Active > Time  >  >  >  >  linear solve 16 5.4663 0.341644 > 100.00  >  >  Totals: 16 5.4663 > 100.00  >  > > > *dr. R. van Tol > *Computer Simulations and Foundry Processes > tel.:+32 9 2645704  fax:+32 9 2645848 > mailto:rob.vantol@... > > *WTCMCRIF Belgian Center of the Technological Industry > *Technologiepark 915, B9052 GentZwijnaarde ( Belgium ) > tel: +32 9 264 5697 > http://www.wtcm.be <http://www.wtcm.be/>; > > *WTCMCRIF "Your Gateway To Innovation". > *Have a look at the future of materials and manufacturing technology > http://techniline.wtcm.be <http://techniline.wtcm.be/>; > 