Menu
▾
▴
Re: [Libmesh-users] Strange results for example 13 only
|
From: Rob v. T. <Rob...@wt...> - 2004-05-06 14:27:06
|
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...@wt...
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
|