On δΊ”, 2011-08-05 at 09:49 -0500, Wolfgang Bangerth wrote:
> > Then I tried to run example program step-17. Everything seemed quit
> > well. But things changed when I invoked step-18. Step-18 gave me the
> > following out put:
> >
> > ========================setp-18 output==============================
> > Timestep 1 at time 1
> > Cycle 0:
> > Number of active cells: 3712 (by partition: 3712)
> > Number of degrees of freedom: 17226 (by partition: 17226)
> > Assembling system... norm of rhs is 2.34224e+10
> > Solver converged in 104 iterations.
> > Updating quadrature point data...
> > Cycle 1:
> > Number of active cells: 12812 (by partition: 12812)
> > Number of degrees of freedom: 51738 (by partition: 51738)
> > Assembling system... norm of rhs is 2.32185e+10
> > Solver converged in 122 iterations.
> > Updating quadrature point data...
> > Moving mesh...
> >
> > Timestep 2 at time 2
> > Assembling system...[0]PETSC ERROR: --------------------- Error
> > Message ------------------------------------
> > [0]PETSC ERROR: Argument out of range!
> > [0]PETSC ERROR: Inserting a new nonzero (10379, 10377) into matrix!
>
> This is a bug that Timo Heister fixed back in February, I believe, and that
> showed up with newer versions of PETSc. If you want to run step-18, you may
> have to download the current development sources of deal.II, as explained on
> the download page.
Wolfgang Bangerth,
Your advice is great!
I downloaded the deal.II 7.1.pre code under development. Compile it
(although, a tiny error occurs and I get rid of it) and run step-18
example program. The program runs very well and all the 10 cycles are
performed. Thanks for you advice, again. And thank Timo, his work is
perfect.
However, it seems that the screen output of step-18 is not exactly the
same with the tutorial. The number of active cells and the number of
degrees of freedom as well as the norm are not the same with the
tutorial after cycle0 (only cycle 0 gets all the output same with the
tutorial). I suppose that it is all right since the program has done the
job normally.
The output is as follows:
=====output of step-18=====
huyaoyu@ubuntu:~/Downloads/deal.II-dev/examples/step-18$ ./step-18
Timestep 1 at time 1
Cycle 0:
Number of active cells: 3712 (by partition: 3712)
Number of degrees of freedom: 17226 (by partition: 17226)
Assembling system... norm of rhs is 2.34224e+10
Solver converged in 104 iterations.
Updating quadrature point data...
Cycle 1:
Number of active cells: 12812 (by partition: 12812)
Number of degrees of freedom: 51738 (by partition: 51738)
Assembling system... norm of rhs is 2.32185e+10
Solver converged in 122 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 2 at time 2
Assembling system... norm of rhs is 2.29175e+10
Solver converged in 131 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 3 at time 3
Assembling system... norm of rhs is 2.2649e+10
Solver converged in 130 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 4 at time 4
Assembling system... norm of rhs is 2.24188e+10
Solver converged in 116 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 5 at time 5
Assembling system... norm of rhs is 2.22358e+10
Solver converged in 114 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 6 at time 6
Assembling system... norm of rhs is 2.21152e+10
Solver converged in 112 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 7 at time 7
Assembling system... norm of rhs is 2.20981e+10
Solver converged in 111 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 8 at time 8
Assembling system... norm of rhs is 2.23556e+10
Solver converged in 136 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 9 at time 9
Assembling system... norm of rhs is 2.3584e+10
Solver converged in 140 iterations.
Updating quadrature point data...
Moving mesh...
Timestep 10 at time 10
Assembling system... norm of rhs is 2.61692e+10
Solver converged in 162 iterations.
Updating quadrature point data...
Moving mesh...
======end of output of step-18==========================
Beat,
HuYaoyu.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
