From: Rob van Tol <Rob.VanT<ol@wt...>  20040505 14:12:15
Attachments:
Message as HTML

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 WTCMCRIF "Your Gateway To Innovation". Have a look at the future of materials and manufacturing technology http://techniline.wtcm.be 
From: Rob van Tol <Rob.VanT<ol@wt...>  20040505 15:51:50
Attachments:
Message as HTML

Dear John, I have tried to compile the code again: ./configure disablempi disableparmetis disablelaspack make The code compiles, as well as the examples. However, when I try to run the examples I get the error massage: ERROR: Unrecognized solver package: 0 src/numerics/numeric_vector.C, line 66 So I am a bit lost. Kind regards, Rob van Tol. At 09:27 5/05/2004 0500, you wrote: >Rob van Tol writes: > > 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. > >Are you using laspack? We have had some problems with laspack >reporting the correct final residual in the past. > >John dr. R. van Tol Computer Simulations and Foundry Processes tel.:09 2645704  fax:09 2645848 mailto:rob.vantol@... WTCMCRIF Belgian Center of the Technological Industry Technologiepark 915, B9052 GentZwijnaarde tel: 09 264 5697 http://www.wtcm.be WTCMCRIF "Your Gateway To Innovation". Have a look at the future of materials and manufacturing technology http://techniline.wtcm.be 
From: Rob van Tol <Rob.VanT<ol@wt...>  20040505 16:31:15
Attachments:
Message as HTML

Dear John, I have PetSc installed ( version 2.2.0 ), without mpi, for which the tests run correctly. So, I think this may not be the problem. I will check the installation of PetSc again and add MPI. Would it be possible that the absence of MPI is the problem ? I would be grateful if you could run example 13 and test the results against the results on the website and then post your version of ex13.C. This way I will be sure that my version of Libmesh causes the problem and not ex13.C itself. The version of ex13.C I have, comes from: libmesh0.4.2 at Sourceforge, which is slightly different from the version on the website. It is remarkable that only example 13 gives different results. Kind regards, Rob van Tol At 11:12 5/05/2004 0500, you wrote: >Rob van Tol writes: > > Dear John, > > > > I have tried to compile the code again: > > > > ./configure disablempi disableparmetis disablelaspack > > make > > > > The code compiles, as well as the examples. However, when I try to run > the > > examples I get the error massage: >Yes, you cannot solve any problem without a linear algebra package. You MUST >have either PETSc or laspack enabled. I would recommend downloading and >installing >PETSc, it works quite well. > >John dr. R. van Tol Computer Simulations and Foundry Processes tel.:09 2645704  fax:09 2645848 mailto:rob.vantol@... WTCMCRIF Belgian Center of the Technological Industry Technologiepark 915, B9052 GentZwijnaarde tel: 09 264 5697 http://www.wtcm.be WTCMCRIF "Your Gateway To Innovation". Have a look at the future of materials and manufacturing technology http://techniline.wtcm.be 
From: Rob van Tol <Rob.VanT<ol@wt...>  20040505 17:03:04
Attachments:
Message as HTML

Dear John, I will get the latest version from the cvs tomorrow, since its getting late at my end of the world, and hope for better results. Thank you for your time and support. I'll keep you informed. Best regards, ROB At 11:53 5/05/2004 0500, you wrote: >Rob van Tol writes: > > I have PetSc installed ( version 2.2.0 ), without mpi, for which the > tests > > run correctly. So, I think this may not be the problem. I will check the > > installation of PetSc again and add MPI. Would it be possible that the > > absence of MPI is the problem ? > >I would definitely check out the latest version of the code from CVS. There >have been a lot of changes since the last official release. We are planning >a new release "any day now" :) > > > >Here are my results for ex13 using PETSc, they are similar to those on the >web page: >*************************************************************** >* 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 >Linear solver converged at step: 16, final residual: 0.121329 Nonlinear >convergence: u  u_old = 223.471 >Linear solver converged at step: 24, final residual: >0.000317615 Nonlinear convergence: u  u_old = 6.531 7 >Linear solver converged at step: 28, final residual: >2.58051e07 Nonlinear convergence: u  u_old = 0.007 2933 > Nonlinear solver converged at step 2 > Solving time step 1, time = 0.01 >Linear solver converged at step: 23, final residual: 0.0025053 Nonlinear >convergence: u  u_old = 35.4701 >Linear solver converged at step: 24, final residual: >1.47282e05 Nonlinear convergence: u  u_old = 0.046 6908 >Linear solver converged at step: 24, final residual: >1.27369e08 Nonlinear convergence: u  u_old = 0.000 27294 > Nonlinear solver converged at step 2 > Solving time step 2, time = 0.015 >Linear solver converged at step: 24, final residual: >0.000790085 Nonlinear convergence: u  u_old = 9.794 14 >Linear solver converged at step: 24, final residual: >3.41055e06 Nonlinear convergence: u  u_old = 0.020 8939 >Linear solver converged at step: 28, final residual: >2.96023e09 Nonlinear convergence: u  u_old = 0.000 135564 > Nonlinear solver converged at step 2 > Solving time step 3, time = 0.02 >Linear solver converged at step: 24, final residual: >0.000382917 Nonlinear convergence: u  u_old = 4.07167 >Linear solver converged at step: 25, final residual: >1.42401e06 Nonlinear convergence: u  u_old = 0.00672638 > Nonlinear solver converged at step 1 > Solving time step 4, time = 0.025 >Linear solver converged at step: 23, final residual: >0.000247016 Nonlinear convergence: u  u_old = 2.13266 >Linear solver converged at step: 25, final residual: >8.20034e07 Nonlinear convergence: u  u_old = 0.00772513 > Nonlinear solver converged at step 1 > Solving time step 5, time = 0.03 >Linear solver converged at step: 23, final residual: >0.000126994 Nonlinear convergence: u  u_old = 1.26206 >Linear solver converged at step: 27, final residual: >4.97111e07 Nonlinear convergence: u  u_old = 0.0057915 > Nonlinear solver converged at step 1 > Solving time step 6, time = 0.035 >Linear solver converged at step: 24, final residual: >8.03056e05 Nonlinear convergence: u  u_old = 0.805059 >Linear solver converged at step: 24, final residual: >2.93058e07 Nonlinear convergence: u  u_old = 0.00147397 > Nonlinear solver converged at step 1 > Solving time step 7, time = 0.04 >Linear solver converged at step: 24, final residual: >7.94782e05 Nonlinear convergence: u  u_old = 0.53619 >Linear solver converged at step: 23, final residual: >1.75705e07 Nonlinear convergence: u  u_old = 0.000919154 > Nonlinear solver converged at step 1 > Solving time step 8, time = 0.045 >Linear solver converged at step: 23, final residual: >6.92094e05 Nonlinear convergence: u  u_old = 0.368994 >Linear solver converged at step: 27, final residual: >9.13022e08 Nonlinear convergence: u  u_old = 0.00104694 > Nonlinear solver converged at step 1 > Solving time step 9, time = 0.05 >Linear solver converged at step: 23, final residual: >4.73186e05 Nonlinear convergence: u  u_old = 0.260159 >Linear solver converged at step: 27, final residual: >6.18452e08 Nonlinear convergence: u  u_old = 0.000762091 > Nonlinear solver converged at step 1 > Solving time step 10, time = 0.055 >Linear solver converged at step: 22, final residual: >4.01872e05 Nonlinear convergence: u  u_old = 0.187318 >Linear solver converged at step: 27, final residual: >3.72042e08 Nonlinear convergence: u  u_old = 0.000711503 > Nonlinear solver converged at step 1 > Solving time step 11, time = 0.06 >Linear solver converged at step: 22, final residual: 2.7971e05 Nonlinear >convergence: u  u_old = 0.137235 >Linear solver converged at step: 28, final residual: >2.40094e08 Nonlinear convergence: u  u_old = 0.000617341 > Nonlinear solver converged at step 1 > Solving time step 12, time = 0.065 >Linear solver converged at step: 21, final residual: >2.51357e05 Nonlinear convergence: u  u_old = 0.102129 >Linear solver converged at step: 28, final residual: >2.11058e08 Nonlinear convergence: u  u_old = 0.00086522 > Nonlinear solver converged at step 1 > Solving time step 13, time = 0.07 >Linear solver converged at step: 21, final residual: >1.78319e05 Nonlinear convergence: u  u_old = 0.07677 >Linear solver converged at step: 28, final residual: >1.60542e08 Nonlinear convergence: u  u_old = 0.000734319 > Nonlinear solver converged at step 1 > Solving time step 14, time = 0.075 >Linear solver converged at step: 21, final residual: >1.27677e05 Nonlinear convergence: u  u_old = 0.0582848 >Linear solver converged at step: 28, final residual: >1.17971e08 Nonlinear convergence: u  u_old = 0.000574086 > Nonlinear solver converged at step 1 > >  > Time: Wed May 5 11:51:51 2004 > OS: Linux > HostName: arthur > OS Release 2.4.2019.9smp > OS Version: #1 SMP Tue Jul 15 17:04:18 EDT 2003 > Machine: i686 > Username: peterson >  >  > Example 13 Performance: Alive time=53.6325, Active time=50.724 >  > Event nCalls Total Avg Percent of  > Time Time Active Time  > >  > linear solve 33 50.7240 1.537090 100.00  >  > Totals: 33 50.7240 100.00  >  > > >  > Reference count information  >  > 12SparseMatrixIdE reference count information: > Creations: 1 > Destructions: 1 > 13NumericVectorIdE reference count information: > Creations: 6 > Destructions: 6 > 21LinearSolverInterfaceIdE reference count information: > Creations: 1 > Destructions: 1 > 4Elem reference count information: > Creations: 4640 > Destructions: 4640 > 4Node reference count information: > Creations: 1681 > Destructions: 1681 > 5QBase reference count information: > Creations: 66 > Destructions: 66 > 6DofMap reference count information: > Creations: 1 > Destructions: 1 > 6FEBase reference count information: > Creations: 66 > Destructions: 66 > 6System reference count information: > Creations: 1 > Destructions: 1 >  > >*************************************************************** >* Done Running Example ./ex13 >*************************************************************** >make: Leaving directory `/home/peterson/code/libmesh/examples/ex13' dr. R. van Tol Computer Simulations and Foundry Processes tel.:09 2645704  fax:09 2645848 mailto:rob.vantol@... WTCMCRIF Belgian Center of the Technological Industry Technologiepark 915, B9052 GentZwijnaarde tel: 09 264 5697 http://www.wtcm.be WTCMCRIF "Your Gateway To Innovation". Have a look at the future of materials and manufacturing technology http://techniline.wtcm.be 
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/>; > 