Dear Libmesh users,

I downloaded verion 0.4.2 of Libmesh from sourceforge and found that example 13 did not run correctly. I did not find the problem, but downloaded the version from the CVS. That version gave realistic results.
        The exact output of the program still differs from the output shown on the website. This may be is due to the iterative solver, because I have used laspack instead of petsc. It's remarkable though that the linear solver always converges in 30 steps.

        If I try to compile the CVS version without laspack ( "./configure --disable-mpi --disable-parmetis --disable-laspack") I get an compilation error, which I did not find for the 0.4.2 version. The error message can be found below.
        The problem does not occur when I use "./configure --disable-mpi --disable-parmetis" or "./configure --disable-mpi --disable-parmetis --disable-petsc". The output is in for both these versions of Limesh the same, also shown below.



Kind regards,

Rob van Tol

Compilation error:
------------------------------------------------------------------------------
Compiling C++ (in optimized mode) src/numerics/analytic_function.C...
Compiling C++ (in optimized mode) src/numerics/coupling_matrix.C...
Compiling C++ (in optimized mode) src/numerics/dense_matrix_base.C...
Compiling C++ (in optimized mode) src/numerics/dense_matrix.C...
Compiling C++ (in optimized mode) src/numerics/dense_submatrix.C...
Compiling C++ (in optimized mode) src/numerics/dense_subvector.C...
Compiling C++ (in optimized mode) src/numerics/dense_vector_base.C...
Compiling C++ (in optimized mode) src/numerics/dense_vector.C...
Compiling C++ (in optimized mode) src/numerics/distributed_vector.C...
Compiling C++ (in optimized mode) src/numerics/error_estimator.C...
Compiling C++ (in optimized mode) src/numerics/function_base.C...
Compiling C++ (in optimized mode) src/numerics/laspack_interface.C...
Compiling C++ (in optimized mode) src/numerics/laspack_matrix.C...
Compiling C++ (in optimized mode) src/numerics/laspack_vector.C...
Compiling C++ (in optimized mode) src/numerics/linear_solver_interface.C...
src/numerics/linear_solver_interface.C: In static member function `static
   AutoPtr<LinearSolverInterface<T> >
   LinearSolverInterface<T>::build(libMeshEnums::SolverPackage) [with T =
   Number]':
src/numerics/linear_solver_interface.C:74:   instantiated from here
src/numerics/linear_solver_interface.C:36: return type `struct
   AutoPtr<LinearSolverInterface<Number> >' is incomplete
src/numerics/linear_solver_interface.C:66: `ap' has incomplete type
src/numerics/linear_solver_interface.C:66: storage size of `ap' isn't known
make: *** [src/numerics/linear_solver_interface.i686-pc-linux-gnu.o] Error 1
------------------------------------------------------------------------------


Ex13 Output
------------------------------------------------------------------------------
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 "Navier-Stokes"
    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: 30, final residual: 2.31901e-06  Nonlinear conv
ergence: ||u - u_old|| = 224.418
Linear solver converged at step: 30, final residual: 3.26614e-06  Nonlinear conv
ergence: ||u - u_old|| = 0.453347
Linear solver converged at step: 30, final residual: 3.26574e-06  Nonlinear conv
ergence: ||u - u_old|| = 1.09929e-06
 Nonlinear solver converged at step 2
 Solving time step 1, time = 0.01
Linear solver converged at step: 30, final residual: 2.02832e-06  Nonlinear conv
ergence: ||u - u_old|| = 35.4333
Linear solver converged at step: 30, final residual: 2.48543e-06  Nonlinear conv
ergence: ||u - u_old|| = 0.0422274
Linear solver converged at step: 30, final residual: 2.48532e-06  Nonlinear conv
ergence: ||u - u_old|| = 5.53181e-07
 Nonlinear solver converged at step 2
 Solving time step 2, time = 0.015
Linear solver converged at step: 30, final residual: 1.82014e-06  Nonlinear conv
ergence: ||u - u_old|| = 9.77147
Linear solver converged at step: 30, final residual: 2.10175e-06  Nonlinear conv
ergence: ||u - u_old|| = 0.0130765
Linear solver converged at step: 30, final residual: 2.10169e-06  Nonlinear conv
ergence: ||u - u_old|| = 4.21478e-07
 Nonlinear solver converged at step 2
 Solving time step 3, time = 0.02
Linear solver converged at step: 30, final residual: 1.76004e-06  Nonlinear conv
ergence: ||u - u_old|| = 4.06722
Linear solver converged at step: 30, final residual: 1.95086e-06  Nonlinear conv
ergence: ||u - u_old|| = 0.00557085
 Nonlinear solver converged at step 1
 Solving time step 4, time = 0.025
 
[vlr@pc02023 ex13]$ ./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 "Navier-Stokes"
    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: 30, final residual: 2.31901e-06  Nonlinear convergence: ||u - u_old|| = 224.418
Linear solver converged at step: 30, final residual: 3.26614e-06  Nonlinear convergence: ||u - u_old|| = 0.453347
Linear solver converged at step: 30, final residual: 3.26574e-06  Nonlinear convergence: ||u - u_old|| = 1.09929e-06
 Nonlinear solver converged at step 2
 Solving time step 1, time = 0.01
Linear solver converged at step: 30, final residual: 2.02832e-06  Nonlinear convergence: ||u - u_old|| = 35.4333
Linear solver converged at step: 30, final residual: 2.48543e-06  Nonlinear convergence: ||u - u_old|| = 0.0422274
Linear solver converged at step: 30, final residual: 2.48532e-06  Nonlinear convergence: ||u - u_old|| = 5.53181e-07
 Nonlinear solver converged at step 2
 Solving time step 2, time = 0.015
Linear solver converged at step: 30, final residual: 1.82014e-06  Nonlinear convergence: ||u - u_old|| = 9.77147
Linear solver converged at step: 30, final residual: 2.10175e-06  Nonlinear convergence: ||u - u_old|| = 0.0130765
Linear solver converged at step: 30, final residual: 2.10169e-06  Nonlinear convergence: ||u - u_old|| = 4.21478e-07
 Nonlinear solver converged at step 2
 Solving time step 3, time = 0.02
Linear solver converged at step: 30, final residual: 1.76004e-06  Nonlinear convergence: ||u - u_old|| = 4.06722
Linear solver converged at step: 30, final residual: 1.95086e-06  Nonlinear convergence: ||u - u_old|| = 0.00557085
 Nonlinear solver converged at step 1
 Solving time step 4, time = 0.025
Linear solver converged at step: 30, final residual: 1.74311e-06  Nonlinear convergence: ||u - u_old|| = 2.13743
Linear solver converged at step: 30, final residual: 1.87848e-06  Nonlinear convergence: ||u - u_old|| = 0.00283913
 Nonlinear solver converged at step 1
 Solving time step 5, time = 0.03
Linear solver converged at step: 30, final residual: 1.73716e-06  Nonlinear convergence: ||u - u_old|| = 1.26614
Linear solver converged at step: 30, final residual: 1.83615e-06  Nonlinear convergence: ||u - u_old|| = 0.00164791
 Nonlinear solver converged at step 1
 Solving time step 6, time = 0.035
Linear solver converged at step: 30, final residual: 1.73347e-06  Nonlinear convergence: ||u - u_old|| = 0.805055
Linear solver converged at step: 30, final residual: 1.80747e-06  Nonlinear convergence: ||u - u_old|| = 0.00105802
 Nonlinear solver converged at step 1
 Solving time step 7, time = 0.04
Linear solver converged at step: 30, final residual: 1.73007e-06  Nonlinear convergence: ||u - u_old|| = 0.536107
Linear solver converged at step: 30, final residual: 1.78627e-06  Nonlinear convergence: ||u - u_old|| = 0.000731641
 Nonlinear solver converged at step 1
 Solving time step 8, time = 0.045
Linear solver converged at step: 30, final residual: 1.72673e-06  Nonlinear convergence: ||u - u_old|| = 0.368847
Linear solver converged at step: 30, final residual: 1.76992e-06  Nonlinear convergence: ||u - u_old|| = 0.000531848
 Nonlinear solver converged at step 1
 Solving time step 9, time = 0.05
Linear solver converged at step: 30, final residual: 1.72358e-06  Nonlinear convergence: ||u - u_old|| = 0.260101
Linear solver converged at step: 30, final residual: 1.75706e-06  Nonlinear convergence: ||u - u_old|| = 0.000398898
 Nonlinear solver converged at step 1
 Solving time step 10, time = 0.055
Linear solver converged at step: 30, final residual: 1.72076e-06  Nonlinear convergence: ||u - u_old|| = 0.18707
Linear solver converged at step: 30, final residual: 1.74686e-06  Nonlinear convergence: ||u - u_old|| = 0.000304901
 Nonlinear solver converged at step 1
 Solving time step 11, time = 0.06
Linear solver converged at step: 30, final residual: 1.71832e-06  Nonlinear convergence: ||u - u_old|| = 0.136781
Linear solver converged at step: 30, final residual: 1.73876e-06  Nonlinear convergence: ||u - u_old|| = 0.000235734
 Nonlinear solver converged at step 1
 Solving time step 12, time = 0.065
Linear solver converged at step: 30, final residual: 1.71627e-06  Nonlinear convergence: ||u - u_old|| = 0.101432
Linear solver converged at step: 30, final residual: 1.73232e-06  Nonlinear convergence: ||u - u_old|| = 0.000183548
 Nonlinear solver converged at step 1
 Solving time step 13, time = 0.07
Linear solver converged at step: 30, final residual: 1.71458e-06  Nonlinear convergence: ||u - u_old|| = 0.07614
Linear solver converged at step: 30, final residual: 1.72721e-06  Nonlinear convergence: ||u - u_old|| = 0.000143559
 Nonlinear solver converged at step 1
 Solving time step 14, time = 0.075
Linear solver converged at step: 30, final residual: 1.7132e-06  Nonlinear convergence: ||u - u_old|| = 0.0577538
Linear solver converged at step: 30, final residual: 1.72316e-06  Nonlinear convergence: ||u - u_old|| = 0.000112617
 Nonlinear solver converged at step 1
 
 ----------------------------------------------------------------------------
| Time:           Thu May  6 16:16:06 2004
| OS:             Linux
| HostName:       pc02023
| OS Release      2.4.20-8
| OS Version:     #1 Thu Mar 13 17:54:28 EST 2003
| Machine:        i686
| Username:       vlr
 ----------------------------------------------------------------------------
 ----------------------------------------------------------------------------
| Example 13 Performance: Alive time=26.6235, Active time=23.3177
 ----------------------------------------------------------------------------
| Event                         nCalls  Total       Avg         Percent of   |
|                                       Time        Time        Active Time  |
|----------------------------------------------------------------------------|
|                                                                            |
| linear solve                  33      23.3177     0.706598    100.00       |
 ----------------------------------------------------------------------------
| Totals:                       33      23.3177                 100.00       |
 ----------------------------------------------------------------------------

------------------------------------------------------------------------------

dr. R. van Tol
Computer Simulations and Foundry Processes
tel.:09 2645704 - fax:09 2645848
mailto:rob.vantol@wtcm.be

WTCM-CRIF Belgian Center of the Technological Industry
Technologiepark 915, B-9052 Gent-Zwijnaarde
tel: 09 264 5697
http://www.wtcm.be

WTCM-CRIF "Your Gateway To Innovation".
Have a look at the future of materials and manufacturing technology
http://techniline.wtcm.be