they all share the same system matrix. This can be exploited in the program.
However, given the nonstandard boundary conditions, to be able to take them into account we need to use
-the following identity
+the following %identity
@f[
\Delta u = \nabla\nabla\cdot u - \nabla\times\nabla\times u,
@f]
# The initial and final time, and the Reynolds number
set initial_time = 0.
set final_time = 10.
- set Reynolds = 1e3
+ set Reynolds = 100
end
subsection Time step data
# In this section we declare the parameters that are going to control the solution process
# for the velocity.
set max_iterations = 1000 # maximal number of iterations that GMRES must make
- set eps = 1e-6 # stopping criterion
+ set eps = 1e-8 # stopping criterion
set Krylov_size = 30 # size of the Krylov subspace to be used in GMRES
- set off_diagonals = 60 # number of off diagonals that ILU must compute
- set diag_strength = 0.01 # diagonal strengthening value
+ set off_diagonals = 70 # number of off diagonals that ILU must compute
+ set diag_strength = 0.1 # diagonal strengthening value
set update_prec = 10 # this number indicates how often the preconditioner must be updated
end
#The output frequency
-set output = 5
+set output = 50
#Finally we set the verbosity level
-set verbose = true
\ No newline at end of file
+set verbose = false
GridIn<dim> grid_in;
grid_in.attach_triangulation (triangulation);
- std::string filename = "nsbench2.inp";
+ std::string filename = "../nsbench2.inp";
std::ifstream file (filename.c_str());
Assert (file, ExcFileNotOpen (filename.c_str()));
grid_in.read_ucd (file);
try
{
RunTimeParameters::Data_Storage data;
- data.read_data ("parameter-file.prm");
+ data.read_data ("../parameter-file.prm");
deallog.depth_console (data.verbose ? 2 : 0);
NavierStokesProjection<2> test (data);
test.run (data.verbose, data.output);