#include <fstream>
#include <string>
-#define PRECISION 5
-
template<int dim>
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> stretch (const Point<dim> &p)
{
- Point<dim> q = p ;
- q(dim-1) *= 2. ;
+ Point<dim> q = p;
+ q(dim-1) *= 2.;
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 5.);
- GridTools::transform (&stretch<dim>, triangulation) ;
+ GridTools::transform (&stretch<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> affine (const Point<dim> &p)
{
- Point<dim> q = p ;
+ Point<dim> q = p;
if (dim>=2)
{
- q(0) = 2.*p(0) + p(1) ;
- q(1) = p(0) + 3.*p(1) ;
+ q(0) = 2.*p(0) + p(1);
+ q(1) = p(0) + 3.*p(1);
}
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 1.);
- GridTools::transform (&affine<dim>, triangulation) ;
+ GridTools::transform (&affine<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
-// check the P1NC element on a genuine quadrilateral.
+// check the P1NC element on a general quadrilateral.
#include "../tests.h"
#include <deal.II/base/logstream.h>
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> bilinear (const Point<dim> &p)
{
- Point<dim> q = p ;
+ Point<dim> q = p;
if (dim>=2)
{
- q(0) = 16.*p(0) +4.*p(1) -10.*p(0)*p(1) ;
- q(1) = 6.*p(1) +4.*p(0)*p(1) ;
+ q(0) = 16.*p(0) +4.*p(1) -10.*p(0)*p(1);
+ q(1) = 6.*p(1) +4.*p(0)*p(1);
}
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 1.);
- GridTools::transform (&bilinear<dim>, triangulation) ;
+ GridTools::transform (&bilinear<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
-
template<int dim>
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> stretch (const Point<dim> &p)
{
- Point<dim> q = p ;
- q(dim-1) *= 2. ;
+ Point<dim> q = p;
+ q(dim-1) *= 2.;
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 5.);
- GridTools::transform (&stretch<dim>, triangulation) ;
+ GridTools::transform (&stretch<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> affine (const Point<dim> &p)
{
- Point<dim> q = p ;
+ Point<dim> q = p;
if (dim>=2)
{
- q(0) = 2.*p(0) + p(1) ;
- q(1) = p(0) + 3.*p(1) ;
+ q(0) = 2.*p(0) + p(1);
+ q(1) = p(0) + 3.*p(1);
}
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 1.);
- GridTools::transform (&affine<dim>, triangulation) ;
+ GridTools::transform (&affine<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-
#include <fstream>
#include <string>
-#define PRECISION 5
template<int dim>
Point<dim> bilinear (const Point<dim> &p)
{
- Point<dim> q = p ;
+ Point<dim> q = p;
if (dim>=2)
{
- q(0) = 16.*p(0) +4.*p(1) -10.*p(0)*p(1) ;
- q(1) = 6.*p(1) +4.*p(0)*p(1) ;
+ q(0) = 16.*p(0) +4.*p(1) -10.*p(0)*p(1);
+ q(1) = 6.*p(1) +4.*p(0)*p(1);
}
- return q ;
+ return q;
}
{
Triangulation<dim> triangulation;
GridGenerator::hyper_cube (triangulation, 0., 1.);
- GridTools::transform (&bilinear<dim>, triangulation) ;
+ GridTools::transform (&bilinear<dim>, triangulation);
FE_P1NC fe;
DoFHandler<dim> dof_handler (triangulation);
main()
{
std::ofstream logfile ("output");
- deallog << std::setprecision(PRECISION);
+ deallog << std::setprecision(5);
deallog << std::fixed;
deallog.attach(logfile);
deallog.depth_console(0);
check<2>();
}
-
-
-