#include "../tests.h"
-std::ofstream logfile("output");
template <typename T>
bool
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
#include "../tests.h"
-std::ofstream logfile("output");
-
-
template <int dim, int fe_degree>
void
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
const unsigned int degree_p = 1;
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
namespace Assembly
{
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "create_mesh.h"
-std::ofstream logfile("output");
template <int dim>
void
#include "create_mesh.h"
-std::ofstream logfile("output");
template <int dim>
void
#include "create_mesh.h"
-std::ofstream logfile("output");
-
template <int dim>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "create_mesh.h"
-std::ofstream logfile("output");
-
template <int dim>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
// the heap (using AlignedVector) instead of allocating it on the stack. Tests
// also copy constructors of FEEvaluation.
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/aligned_vector.h>
#include <deal.II/base/utilities.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
template <int dim, int degree_p, typename VectorType>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
// We want to use the matrix-vector product provided by this function (which
#include "../tests.h"
-std::ofstream logfile("output");
-
template <int dim>
// different polynomial degrees. The mesh uses a hypercube mesh with no
// hanging nodes and no other constraints
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/function_lib.h>
#include <deal.II/fe/fe_dgq.h>
#include <deal.II/numerics/vector_tools.h>
+#include "../tests.h"
+
void
output_double_number(double input, const std::string &text)
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(2);
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename number>
class MatrixFreeTest
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename number>
class MatrixFreeTest
#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
int
main()
{
- deallog.attach(logfile);
- deallog.depth_console(0);
+ initlog();
deallog << std::setprecision(3);
{
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
template <int dim, int fe_degree, typename number>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
template <int dim,
int fe_degree,
int
main()
{
- deallog.attach(logfile);
- // need to set quite a loose tolerance because
- // FEValues approximates Hessians with finite
- // differences, which are not so accurate
- deallog << std::setprecision(3);
+ initlog();
+ deallog << std::setprecision(7);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
-
template <int dim,
int fe_degree,
int
main()
{
- deallog.attach(logfile);
- // need to set quite a loose tolerance because
- // FEValues approximates Hessians with finite
- // differences, which are not so accurate
- deallog << std::setprecision(3);
+ initlog();
+ deallog << std::setprecision(7);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "get_functions_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename Number>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename Number>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
{
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename Number>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
{
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename Number>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
{
#include "create_mesh.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree, typename Number>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
{
#include "create_mesh.h"
-std::ofstream logfile("output");
template <int dim>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
test<2>();
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
-
template <int dim,
int fe_degree,
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
-
template <int dim,
int fe_degree,
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
int
main()
{
- deallog.attach(logfile);
- deallog.depth_console(0);
-
- deallog << std::setprecision(3);
+ initlog();
{
deallog.push("2d");
// vector-valued problem (curl-curl operator which does not really make a lot
// of sense from a problem point of view, though).
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
const double global_coefficient = 0.1;
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
// vector-valued problem (div-div operator which does not really make a lot
// of sense from a problem point of view, though).
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
#include "create_mesh.h"
const double global_coefficient = 0.1;
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
template <int dim, int fe_degree>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
#include "../tests.h"
-std::ofstream logfile("output");
-
#include "create_mesh.h"
#include "matrix_vector_common.h"
// matrix for hp DoFHandler on a hyperball mesh with hanging nodes and finite
// elements orders distributed randomly.
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/function.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/fe_values.h>
+#include "../tests.h"
+
#include "matrix_vector_common.h"
// matrix-based results in this case, this test simply checks that the
// matrix-vector product runs without error
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/affine_constraints.h>
#include <deal.II/lac/vector.h>
+#include "../tests.h"
+
#include "matrix_vector_mf.h"
// have a constant-coefficient Laplacian, we can verify the implementation by
// comparing to the result with 21 quadrature points.
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/affine_constraints.h>
#include <deal.II/lac/vector.h>
+#include "../tests.h"
+
#include "matrix_vector_mf.h"
// matrix for MG DoFHandler on a hyperball mesh with no hanging nodes but
// homogeneous Dirichlet conditions
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/function.h>
#include <deal.II/multigrid/mg_matrix.h>
#include <deal.II/multigrid/mg_tools.h>
+#include "../tests.h"
+
#include "matrix_vector_common.h"
// matrix. No hanging nodes and no other
// constraints for a vector-valued problem (stokes equations).
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
#include "create_mesh.h"
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
// the Stokes equations. Like matrix_vector_stokes_onedof except for
// different constraints
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/logstream.h>
#include <deal.II/base/utilities.h>
#include <fstream>
#include <iostream>
+#include "../tests.h"
+
template <int dim, int degree_p, typename VectorType>
// between the vector components in the form of no-normal flux constraints on
// the Stokes equations.
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
template <int dim, int degree_p, typename VectorType>
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
// same as matrix_vector_stokes_noflux but no template parameter on the
// polynomial degree
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
template <int dim, typename VectorType>
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
// putting all degrees of freedom into a single DoFHandler, where the
// selection is done through FEEvaluation
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/logstream.h>
#include <deal.II/base/utilities.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
template <int dim, int degree_p, typename VectorType>
// matrix. No hanging nodes and no other
// constraints for a vector-valued problem (stokes equations).
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
#include "create_mesh.h"
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
// problem (stokes equations). Same as matrix_vector_stokes_qdg0 but now
// without the template on the degree of the element.
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <complex>
#include <iostream>
+#include "../tests.h"
+
#include "create_mesh.h"
int
main()
{
- deallog.attach(logfile);
-
+ initlog();
deallog << std::setprecision(3);
{
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
// check that FEEvaluation can be evaluated without indices initialized (and
// throws an exception when trying to read/write from/to vectors)
-#include "../tests.h"
-
-
-std::ofstream logfile("output");
-
#include <deal.II/base/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <iostream>
+#include "../tests.h"
+
// forward declare this function
template <int dim, int fe_degree>
main()
{
deal_II_exceptions::disable_abort_on_exception();
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
test<2, 1>();
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
-
template <typename MatrixType, typename Number>
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim,
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
{
deallog.push("2d");
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
main(int argc, char **argv)
{
Utilities::MPI::MPI_InitFinalize mpi_init(argc, argv, 1);
-
- if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
- {
- deallog.attach(logfile);
- deallog << std::setprecision(4);
- }
+ mpi_initlog();
+ deallog << std::setprecision(4);
// 2 blocks
{
#include "../tests.h"
-std::ofstream logfile("output");
using namespace dealii::MatrixFreeOperators;
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim, int fe_degree>
int
main()
{
- deallog.attach(logfile);
+ initlog();
deallog << std::setprecision(3);
{
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim>
#include "../tests.h"
-std::ofstream logfile("output");
template <int dim>
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_common.h"
#include "../tests.h"
#include "create_mesh.h"
-
-std::ofstream logfile("output");
-
#include "matrix_vector_faces_common.h"
// matrix free matrix-vector products for hp elements by comparing to the
// serial version
-#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include <deal.II/base/function.h>
#include <deal.II/base/template_constraints.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/fe_values.h>
+#include "../tests.h"
+
#include "create_mesh.h"
#include "matrix_vector_common.h"
#include "../tests.h"
-
-std::ofstream logfile("output");
-
#include "get_functions_common.h"