17 #ifndef __deal2__dof_print_solver_step_h
18 #define __deal2__dof_print_solver_step_h
20 #include <deal.II/base/config.h>
21 #include <deal.II/base/logstream.h>
22 #include <deal.II/lac/solver_control.h>
23 #include <deal.II/lac/vector_memory.h>
24 #include <deal.II/numerics/data_out.h>
30 DEAL_II_NAMESPACE_OPEN
54 template<
int dim,
class SOLVER,
class VECTOR = Vector<
double> >
82 const VECTOR &d)
const;
98 template<
int dim,
class SOLVER,
class VECTOR>
102 const std::string &basename)
109 template<
int dim,
class SOLVER,
class VECTOR>
114 const VECTOR &d)
const
116 out.clear_data_vectors();
117 out.add_data_vector(x,
"solution");
118 out.add_data_vector(r,
"residual");
119 out.add_data_vector(d,
"update");
121 std::ostringstream filename;
123 << std::setw(3) << std::setfill(
'0') << step
124 << out.default_suffix();
126 const std::string fname = filename.str();
128 deallog <<
"Writing file:" << fname << std::endl;
131 std::ofstream of (fname.c_str());
135 DEAL_II_NAMESPACE_CLOSE
const std::string basename
DoFPrintSolverStep(SolverControl &control, VectorMemory< VECTOR > &mem, DataOut< dim > &data_out, const std::string &basename)
virtual void print_vectors(const unsigned int step, const VECTOR &x, const VECTOR &r, const VECTOR &d) const