//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<1,dim> >(fe.n_components()));
fe_values.get_function_gradients (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
{
FEValuesExtractors::Scalar single_component (c);
fe_values[single_component].get_function_gradients (fe_function,
scalar_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << scalar_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ deallog << scalar_values[q][d] << (d<dim-1 ? " " : "");
+ deallog << std::endl;
Assert ((scalar_values[q] - vector_values[q][c]).norm()
<= 1e-12 * scalar_values[q].norm(),
ExcInternalError());
DEAL::125. -61.1
DEAL::125. -40.8
DEAL::component=3
-DEAL::0.00 5.84
-DEAL::0.00 5.84
-DEAL::0.00 9.52
-DEAL::0.00 9.52
+DEAL::0 5.84
+DEAL::0 5.84
+DEAL::0 9.52
+DEAL::0 9.52
DEAL::component=4
DEAL::5.84 -6.75
DEAL::5.84 6.75
DEAL::2.37 4.73 9.46
DEAL::2.37 4.73 9.46
DEAL::component=1
-DEAL::-2.21e-12 -275. -138.
-DEAL::2.77e-12 -275. -138.
+DEAL::0 -275. -138.
+DEAL::0 -275. -138.
DEAL::-367. -275. 4.71e+03
DEAL::367. -275. 4.92e+03
DEAL::-184. 4.57e+03 -138.
DEAL::4.80e+03 4.57e+03 4.71e+03
DEAL::-4.38e+03 4.79e+03 4.92e+03
DEAL::component=2
-DEAL::-275. -8.99e-13 -138.
+DEAL::-275. 0 -138.
DEAL::-275. -459. 4.80e+03
-DEAL::-275. 1.41e-12 -138.
+DEAL::-275. 0 -138.
DEAL::-275. 459. 5.02e+03
DEAL::4.67e+03 -91.8 -138.
DEAL::4.67e+03 3.52e+03 4.80e+03
DEAL::4.88e+03 91.8 -138.
DEAL::4.88e+03 -3.09e+03 5.02e+03
DEAL::component=3
-DEAL::-275. -138. 1.36e-12
+DEAL::-275. -138. 0
DEAL::-275. 4.90e+03 -367.
DEAL::4.76e+03 -138. -184.
DEAL::4.76e+03 4.90e+03 2.23e+03
-DEAL::-275. -138. -1.36e-12
+DEAL::-275. -138. 0
DEAL::-275. 5.11e+03 367.
DEAL::4.97e+03 -138. 184.
DEAL::4.97e+03 5.11e+03 -1.81e+03
DEAL::component=4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::component=5
-DEAL::5.60 -1.84e-16 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 -1.84e-16 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 7.35e-16 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 7.35e-16 22.4
-DEAL::5.60 0.00 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 0 22.4
DEAL::component=6
-DEAL::5.60 11.2 2.53e-16
-DEAL::5.60 11.2 1.01e-15
-DEAL::5.60 11.2 1.01e-15
-DEAL::5.60 11.2 0.00
-DEAL::5.60 11.2 2.53e-16
-DEAL::5.60 11.2 1.01e-15
-DEAL::5.60 11.2 1.01e-15
-DEAL::5.60 11.2 0.00
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
+DEAL::5.60 11.2 0
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<2,dim> >(fe.n_components()));
fe_values.get_function_hessians (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
{
FEValuesExtractors::Scalar single_component (c);
fe_values[single_component].get_function_hessians (fe_function,
scalar_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << scalar_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << scalar_values[q][d][e] << (d<dim-1 || e<dim-1 ? " " : "");
+ deallog << std::endl;
Assert ((scalar_values[q] - vector_values[q][c]).norm()
<= 1e-12 * scalar_values[q].norm(),
ExcInternalError());
DEAL::FE=FESystem<2>[FE_Q<2>(1)-FE_RaviartThomas<2>(1)-FE_Nedelec<2>(0)]
DEAL::component=0
-DEAL::5.60e-11 1.30 1.30 -1.50
-DEAL::0.00 1.30 1.30 1.50
-DEAL::0.00 3.45 3.45 -3.99
-DEAL::0.00 3.45 3.45 3.99
+DEAL::0 1.30 1.30 -1.50
+DEAL::0 1.30 1.30 1.50
+DEAL::0 3.45 3.45 -3.99
+DEAL::0 3.45 3.45 3.99
DEAL::component=1
DEAL::-50.8 877. 848. -953.
DEAL::-50.8 -779. -889. -641.
DEAL::-1.43e+03 2.86e+03 2.65e+03 -2.63e+03
DEAL::-1.43e+03 -2.60e+03 -2.12e+03 -1.44e+03
DEAL::component=2
-DEAL::3.99e-09 278. 209. -236.
-DEAL::-6.72e-10 278. 209. -247.
-DEAL::-1.25e-08 374. 46.3 -122.
-DEAL::8.76e-10 374. 46.3 -1.16e+03
+DEAL::0 278. 209. -236.
+DEAL::0 278. 209. -247.
+DEAL::0 374. 46.3 -122.
+DEAL::0 374. 46.3 -1.16e+03
DEAL::component=3
-DEAL::0.00 0.00 2.69e-10 9.43
-DEAL::0.00 0.00 8.06e-10 9.43
-DEAL::0.00 0.00 -8.76e-10 25.0
-DEAL::0.00 0.00 8.76e-10 25.0
+DEAL::0 0 0 9.43
+DEAL::0 0 0 9.43
+DEAL::0 0 0 25.0
+DEAL::0 0 0 25.0
DEAL::component=4
-DEAL::2.69e-10 9.43 18.9 -21.8
-DEAL::8.06e-10 9.43 18.9 21.8
-DEAL::-8.76e-10 25.0 50.0 -57.8
-DEAL::8.76e-10 25.0 50.0 57.8
+DEAL::0 9.43 18.9 -21.8
+DEAL::0 9.43 18.9 21.8
+DEAL::0 25.0 50.0 -57.8
+DEAL::0 25.0 50.0 57.8
DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(0)]
DEAL::component=0
-DEAL::0.00 9.85e-10 5.09e-10 1.54e-09 2.44e-26 5.91e-10 3.28e-10 4.60e-10 -2.17e-09
-DEAL::0.00 9.85e-10 1.74e-09 -1.38e-09 0.00 2.53e-09 -1.30e-09 5.58e-10 1.64e-09
-DEAL::4.27e-10 -1.18e-09 1.26e-09 2.89e-09 -8.21e-10 2.89e-09 -8.54e-10 -7.88e-10 1.64e-09
-DEAL::0.00 -1.18e-09 -2.30e-10 -3.94e-10 -8.21e-10 -4.60e-10 -7.88e-10 -7.88e-10 -2.17e-09
-DEAL::0.00 5.09e-10 -1.97e-10 1.97e-10 -1.94e-26 -1.38e-09 3.28e-10 4.60e-10 3.23e-27
-DEAL::0.00 5.09e-10 -2.13e-10 -3.94e-10 1.08e-09 -3.94e-10 -1.63e-09 9.52e-10 -1.11e-26
-DEAL::-5.42e-10 4.92e-11 5.42e-10 1.97e-10 -8.21e-10 -3.94e-10 -2.36e-09 -7.88e-10 -5.18e-26
-DEAL::0.00 4.92e-11 4.93e-11 -3.94e-10 0.00 -1.38e-09 -7.88e-10 -7.88e-10 0.00
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
DEAL::component=1
-DEAL::2.57e-07 -1.50e+03 -752. -1.50e+03 3.49e-07 1.99e+04 -752. 1.99e+04 5.31e-08
-DEAL::-3.47e-07 1.50e+03 752. 1.50e+03 5.81e-07 2.07e+04 752. 2.07e+04 7.10e-07
-DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 -4.08e-07 1.99e+04 2.12e+04 1.99e+04 2.96e-06
-DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 -2.94e-07 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
-DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 -6.29e-08 1.99e+04 -752. 1.99e+04 2.25e-08
-DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 2.06e-06 2.07e+04 752. 2.07e+04 3.03e-13
-DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -8.04e-07 1.99e+04 2.12e+04 1.99e+04 -7.68e-07
-DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 2.10e-08 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
+DEAL::4.13e-07 -1.50e+03 -752. -1.50e+03 0 1.99e+04 -752. 1.99e+04 -9.53e-07
+DEAL::-1.68e-07 1.50e+03 752. 1.50e+03 1.82e-07 2.07e+04 752. 2.07e+04 -3.99e-07
+DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 8.52e-07 1.99e+04 2.12e+04 1.99e+04 3.27e-06
+DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 9.12e-07 2.07e+04 -1.94e+04 2.07e+04 -8.34e-07
+DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 2.73e-07 1.99e+04 -752. 1.99e+04 5.33e-07
+DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 9.89e-07 2.07e+04 752. 2.07e+04 0
+DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -9.17e-07 1.99e+04 2.12e+04 1.99e+04 -4.73e-07
+DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 0 2.07e+04 -1.94e+04 2.07e+04 -3.80e-07
DEAL::component=2
-DEAL::-2.06e-07 -1.88e+03 2.02e+04 -1.88e+03 1.75e-06 -376. 2.02e+04 -376. 4.37e-07
-DEAL::-3.62e-07 -1.88e+03 2.02e+04 -1.88e+03 3.76e+03 1.63e+04 2.02e+04 1.63e+04 2.93e-06
-DEAL::1.88e-08 1.88e+03 2.11e+04 1.88e+03 1.28e-06 376. 2.11e+04 376. 7.40e-07
-DEAL::5.04e-08 1.88e+03 2.11e+04 1.88e+03 3.76e+03 -1.46e+04 2.11e+04 -1.46e+04 -9.08e-07
-DEAL::-8.02e-09 1.48e+04 2.02e+04 1.48e+04 752. -376. 2.02e+04 -376. -5.17e-07
-DEAL::-2.14e-08 1.48e+04 2.02e+04 1.48e+04 -2.71e+04 1.63e+04 2.02e+04 1.63e+04 5.50e-13
-DEAL::-1.46e-06 -1.30e+04 2.11e+04 -1.30e+04 752. 376. 2.11e+04 376. -9.25e-08
-DEAL::0.00 -1.30e+04 2.11e+04 -1.30e+04 -2.71e+04 -1.46e+04 2.11e+04 -1.46e+04 0.00
+DEAL::0 -1.88e+03 2.02e+04 -1.88e+03 1.25e-07 -376. 2.02e+04 -376. 5.04e-07
+DEAL::-5.05e-07 -1.88e+03 2.02e+04 -1.88e+03 3.76e+03 1.63e+04 2.02e+04 1.63e+04 3.26e-06
+DEAL::0 1.88e+03 2.11e+04 1.88e+03 -7.91e-07 376. 2.11e+04 376. -1.34e-07
+DEAL::0 1.88e+03 2.11e+04 1.88e+03 3.76e+03 -1.46e+04 2.11e+04 -1.46e+04 -5.34e-07
+DEAL::0 1.48e+04 2.02e+04 1.48e+04 752. -376. 2.02e+04 -376. -4.41e-07
+DEAL::0 1.48e+04 2.02e+04 1.48e+04 -2.71e+04 1.63e+04 2.02e+04 1.63e+04 0
+DEAL::-1.98e-06 -1.30e+04 2.11e+04 -1.30e+04 752. 376. 2.11e+04 376. 0
+DEAL::0 -1.30e+04 2.11e+04 -1.30e+04 -2.71e+04 -1.46e+04 2.11e+04 -1.46e+04 0
DEAL::component=3
-DEAL::5.81e-08 2.06e+04 -1.50e+03 2.06e+04 3.99e-08 -752. -1.50e+03 -752. -1.93e-06
-DEAL::-2.94e-08 2.06e+04 -1.50e+03 2.06e+04 2.26e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::3.01e-06 2.06e+04 9.89e+03 2.06e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::-9.50e-07 2.06e+04 9.89e+03 2.06e+04 -3.25e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::-1.81e-07 2.15e+04 1.50e+03 2.15e+04 -1.18e-07 752. 1.50e+03 752. -2.35e-06
-DEAL::0.00 2.15e+04 1.50e+03 2.15e+04 1.87e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::-1.34e-06 2.15e+04 -8.16e+03 2.15e+04 1.09e-07 752. -8.16e+03 752. 1.50e+03
-DEAL::0.00 2.15e+04 -8.16e+03 2.15e+04 5.27e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
+DEAL::0 2.06e+04 -1.50e+03 2.06e+04 -4.31e-07 -752. -1.50e+03 -752. 4.48e-07
+DEAL::0 2.06e+04 -1.50e+03 2.06e+04 1.78e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
+DEAL::3.66e-06 2.06e+04 9.89e+03 2.06e+04 -3.93e-07 -752. 9.89e+03 -752. 1.50e+03
+DEAL::-8.57e-07 2.06e+04 9.89e+03 2.06e+04 -2.68e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
+DEAL::-1.81e-07 2.15e+04 1.50e+03 2.15e+04 -4.94e-07 752. 1.50e+03 752. 1.88e-06
+DEAL::0 2.15e+04 1.50e+03 2.15e+04 2.25e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
+DEAL::-1.25e-06 2.15e+04 -8.16e+03 2.15e+04 -1.00e-06 752. -8.16e+03 752. 1.50e+03
+DEAL::0 2.15e+04 -8.16e+03 2.15e+04 3.47e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
DEAL::component=4
-DEAL::0.00 0.00 0.00 1.58e-09 2.63e-09 2.36e-09 0.00 1.05e-08 4.31e-26
-DEAL::0.00 0.00 0.00 -2.63e-10 2.63e-09 6.30e-09 -6.30e-09 1.37e-08 2.10e-08
-DEAL::0.00 0.00 0.00 5.25e-09 -3.94e-09 6.04e-09 8.14e-09 -7.62e-09 1.58e-08
-DEAL::0.00 0.00 0.00 0.00 -3.94e-09 0.00 0.00 -1.05e-08 -5.25e-09
-DEAL::0.00 0.00 0.00 0.00 1.31e-09 -1.58e-09 0.00 4.20e-09 1.73e-26
-DEAL::0.00 0.00 0.00 0.00 1.58e-09 2.63e-10 -8.41e-09 8.41e-09 -3.94e-31
-DEAL::0.00 0.00 0.00 -1.31e-09 -1.31e-09 1.31e-09 -6.57e-09 3.94e-09 5.25e-09
-DEAL::0.00 0.00 0.00 0.00 0.00 -1.31e-09 0.00 0.00 0.00
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
DEAL::component=5
-DEAL::2.35e-09 1.58e-09 2.20e-09 -1.20e-09 4.32e-26 -7.61e-10 1.14e-08 4.33e-25 4.35e-10
-DEAL::-1.31e-09 1.58e-09 6.15e-09 0.00 0.00 6.52e-10 -2.89e-09 0.00 1.58e-08
-DEAL::4.98e-09 -2.89e-09 6.14e-09 -1.20e-09 -7.35e-26 -7.61e-10 1.82e-08 -1.05e-08 2.14e-08
-DEAL::0.00 -2.89e-09 1.09e-10 0.00 0.00 6.52e-10 -1.05e-09 -5.25e-09 -5.25e-09
-DEAL::2.17e-10 1.31e-09 -1.68e-09 1.74e-09 9.45e-26 -8.70e-10 -1.81e-10 -6.81e-26 -5.69e-09
-DEAL::0.00 2.63e-09 -2.63e-10 0.00 0.00 0.00 2.36e-09 5.25e-09 1.25e-25
-DEAL::-4.52e-11 -1.31e-09 9.42e-10 1.74e-09 6.17e-26 -8.70e-10 -1.76e-09 -1.26e-08 -4.35e-10
-DEAL::0.00 0.00 -1.58e-09 0.00 0.00 0.00 -1.05e-09 0.00 -5.25e-09
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
DEAL::component=6
-DEAL::-1.63e-10 2.42e-09 2.38e-26 4.78e-09 -1.63e-10 1.76e-26 -7.07e-10 -4.89e-10 3.15e-26
-DEAL::-4.35e-10 2.58e-09 4.20e-09 -3.37e-09 -2.17e-10 7.88e-09 -3.70e-09 2.17e-10 6.90e-25
-DEAL::1.10e-09 -3.43e-09 3.94e-09 8.10e-09 -3.17e-09 8.41e-09 -2.17e-10 -5.98e-09 6.64e-25
-DEAL::0.00 -2.89e-09 -1.31e-09 -5.25e-10 -2.63e-09 -2.63e-09 0.00 0.00 0.00
-DEAL::-1.63e-10 1.11e-09 -2.51e-27 -4.71e-10 -1.63e-10 -2.63e-09 -7.07e-10 -4.89e-10 -5.96e-26
-DEAL::-4.35e-10 1.27e-09 1.00e-26 -7.43e-10 3.08e-10 -1.08e-26 -3.70e-09 2.17e-10 -5.71e-26
-DEAL::-1.79e-09 -8.06e-10 1.31e-09 -2.41e-09 -3.17e-09 -7.87e-26 -2.17e-10 -5.98e-09 5.87e-26
-DEAL::0.00 -2.63e-10 0.00 -5.25e-10 0.00 -2.63e-09 0.00 0.00 0.00
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
+DEAL::0 0 0 0 0 0 0 0 0
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<2,dim> >(fe.n_components()));
fe_values.get_function_hessians (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
{
FEValuesExtractors::Scalar single_component (c);
fe_values[single_component].get_function_hessians (fe_function,
scalar_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << scalar_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << scalar_values[q][d][e] << (d<dim-1 || e<dim-1 ? " " : "");
+ deallog << std::endl;
+
Assert ((scalar_values[q] - vector_values[q][c]).norm()
<= 1e-12 * scalar_values[q].norm(),
ExcInternalError());
DEAL::-949. 1.90e+03 1.76e+03 -1.73e+03
DEAL::-949. -1.73e+03 -1.39e+03 -877.
DEAL::component=1
-DEAL::2.33e-09 188. 147. -148.
-DEAL::-3.58e-10 188. 147. -204.
-DEAL::-7.59e-09 174. -23.5 33.7
-DEAL::5.84e-10 174. -23.5 -967.
+DEAL::0 188. 147. -148.
+DEAL::0 188. 147. -204.
+DEAL::0 174. -23.5 33.7
+DEAL::0 174. -23.5 -967.
DEAL::FE=FE_RaviartThomas<3>(1)
DEAL::component=0
-DEAL::2.57e-07 -1.50e+03 -752. -1.50e+03 2.23e-07 1.24e+04 -752. 1.24e+04 1.11e-08
-DEAL::2.41e-07 1.50e+03 752. 1.50e+03 3.42e-07 1.32e+04 752. 1.32e+04 4.58e-07
-DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 -1.98e-07 1.24e+04 2.12e+04 1.24e+04 1.83e-06
-DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 -1.68e-07 1.32e+04 -1.94e+04 1.32e+04 -5.46e-07
-DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 -2.08e-08 1.24e+04 -752. 1.24e+04 2.25e-08
-DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 1.23e-06 1.32e+04 752. 1.32e+04 1.80e-13
-DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -5.09e-07 1.24e+04 2.12e+04 1.24e+04 -5.15e-07
-DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 2.10e-08 1.32e+04 -1.94e+04 1.32e+04 -5.46e-07
+DEAL::4.13e-07 -1.50e+03 -752. -1.50e+03 0 1.24e+04 -752. 1.24e+04 -5.75e-07
+DEAL::2.52e-07 1.50e+03 752. 1.50e+03 0 1.32e+04 752. 1.32e+04 -3.15e-07
+DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 5.58e-07 1.24e+04 2.12e+04 1.24e+04 2.05e-06
+DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 5.76e-07 1.32e+04 -1.94e+04 1.32e+04 -5.40e-07
+DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 1.89e-07 1.24e+04 -752. 1.24e+04 3.23e-07
+DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 6.23e-07 1.32e+04 752. 1.32e+04 0
+DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -5.39e-07 1.24e+04 2.12e+04 1.24e+04 -2.63e-07
+DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 0 1.32e+04 -1.94e+04 1.32e+04 -2.54e-07
DEAL::component=1
-DEAL::-1.21e-07 -1.88e+03 1.27e+04 -1.88e+03 1.29e-06 -376. 1.27e+04 -376. 2.69e-07
-DEAL::-2.36e-07 -1.88e+03 1.27e+04 -1.88e+03 3.76e+03 1.63e+04 1.27e+04 1.63e+04 1.84e-06
-DEAL::-1.05e-08 1.88e+03 1.36e+04 1.88e+03 1.03e-06 376. 1.36e+04 376. 4.46e-07
-DEAL::5.04e-08 1.88e+03 1.36e+04 1.88e+03 3.76e+03 -1.46e+04 1.36e+04 -1.46e+04 -5.72e-07
-DEAL::-8.02e-09 1.48e+04 1.27e+04 1.48e+04 752. -376. 1.27e+04 -376. -3.07e-07
-DEAL::-2.14e-08 1.48e+04 1.27e+04 1.48e+04 -2.71e+04 1.63e+04 1.27e+04 1.63e+04 4.27e-13
-DEAL::-9.39e-07 -1.30e+04 1.36e+04 -1.30e+04 752. 376. 1.36e+04 376. -5.04e-08
-DEAL::0.00 -1.30e+04 1.36e+04 -1.30e+04 -2.71e+04 -1.46e+04 1.36e+04 -1.46e+04 0.00
+DEAL::0 -1.88e+03 1.27e+04 -1.88e+03 2.93e-07 -376. 1.27e+04 -376. 2.94e-07
+DEAL::-2.95e-07 -1.88e+03 1.27e+04 -1.88e+03 3.76e+03 1.63e+04 1.27e+04 1.63e+04 2.08e-06
+DEAL::0 1.88e+03 1.36e+04 1.88e+03 -3.29e-07 376. 1.36e+04 376. 0
+DEAL::0 1.88e+03 1.36e+04 1.88e+03 3.76e+03 -1.46e+04 1.36e+04 -1.46e+04 -3.24e-07
+DEAL::0 1.48e+04 1.27e+04 1.48e+04 752. -376. 1.27e+04 -376. -2.73e-07
+DEAL::0 1.48e+04 1.27e+04 1.48e+04 -2.71e+04 1.63e+04 1.27e+04 1.63e+04 0
+DEAL::-1.26e-06 -1.30e+04 1.36e+04 -1.30e+04 752. 376. 1.36e+04 376. 0
+DEAL::0 -1.30e+04 1.36e+04 -1.30e+04 -2.71e+04 -1.46e+04 1.36e+04 -1.46e+04 0
DEAL::component=2
-DEAL::-2.60e-08 1.31e+04 -1.50e+03 1.31e+04 -2.10e-09 -752. -1.50e+03 -752. -1.05e-06
-DEAL::1.26e-08 1.31e+04 -1.50e+03 1.31e+04 1.42e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::1.83e-06 1.31e+04 9.89e+03 1.31e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::-5.72e-07 1.31e+04 9.89e+03 1.31e+04 -2.03e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::-1.10e-07 1.40e+04 1.50e+03 1.40e+04 -7.57e-08 752. 1.50e+03 752. -1.34e-06
-DEAL::0.00 1.40e+04 1.50e+03 1.40e+04 1.20e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::-8.49e-07 1.40e+04 -8.16e+03 1.40e+04 6.72e-08 752. -8.16e+03 752. 1.50e+03
-DEAL::0.00 1.40e+04 -8.16e+03 1.40e+04 3.17e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
+DEAL::0 1.31e+04 -1.50e+03 1.31e+04 -3.05e-07 -752. -1.50e+03 -752. 4.90e-07
+DEAL::0 1.31e+04 -1.50e+03 1.31e+04 1.11e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
+DEAL::2.24e-06 1.31e+04 9.89e+03 1.31e+04 -2.25e-07 -752. 9.89e+03 -752. 1.50e+03
+DEAL::-5.21e-07 1.31e+04 9.89e+03 1.31e+04 -1.67e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
+DEAL::-1.10e-07 1.40e+04 1.50e+03 1.40e+04 -3.26e-07 752. 1.50e+03 752. 1.29e-06
+DEAL::0 1.40e+04 1.50e+03 1.40e+04 1.45e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
+DEAL::-7.99e-07 1.40e+04 -8.16e+03 1.40e+04 -6.22e-07 752. -8.16e+03 752. 1.50e+03
+DEAL::0 1.40e+04 -8.16e+03 1.40e+04 2.21e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<2,dim> >(fe.n_components()));
fe_values.get_function_hessians (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
{
FEValuesExtractors::Scalar single_component (c);
fe_values[single_component].get_function_hessians (fe_function,
scalar_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << scalar_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << scalar_values[q][d][e] << (d<dim-1 || e<dim-1 ? " " : "");
+ deallog << std::endl;
+
Assert ((scalar_values[q] - vector_values[q][c]).norm()
<= 1e-12 * scalar_values[q].norm(),
ExcInternalError());
DEAL::FE=FE_Nedelec<2>(1)
DEAL::component=0
-DEAL::1.03e-25 -17.8 -30.7 14.9
-DEAL::0.00 -17.8 -30.7 234.
-DEAL::0.00 321. 247. -363.
-DEAL::0.00 321. 247. 1.02e+03
+DEAL::0 -17.8 -30.7 14.9
+DEAL::0 -17.8 -30.7 234.
+DEAL::0 321. 247. -363.
+DEAL::0 321. 247. 1.02e+03
DEAL::component=1
DEAL::7.43 -154. -143. 177.
DEAL::7.43 219. 164. 370.
DEAL::406. 1.02e+03 1.17e+03 1.53e+03
DEAL::FE=FE_Nedelec<3>(1)
DEAL::component=0
-DEAL::0.00 -1.26e+03 -714. -1.26e+03 -32.6 -13.2 -714. -13.2 -32.6
-DEAL::0.00 -1.26e+03 -714. -1.26e+03 2.57e+03 1.10e+03 -714. 1.10e+03 2.02e+03
-DEAL::-7.62e-08 1.34e+03 -820. 1.34e+03 -32.6 13.2 -820. 13.2 -32.6
-DEAL::0.00 1.34e+03 -820. 1.34e+03 2.57e+03 -1.31e+03 -820. -1.31e+03 2.24e+03
-DEAL::0.00 -1.37e+03 1.34e+03 -1.37e+03 -32.6 13.2 1.34e+03 13.2 -32.6
-DEAL::0.00 -1.37e+03 1.34e+03 -1.37e+03 2.79e+03 -1.31e+03 1.34e+03 -1.31e+03 2.02e+03
-DEAL::8.30e-08 1.45e+03 1.45e+03 1.45e+03 -32.6 -13.2 1.45e+03 -13.2 -32.6
-DEAL::0.00 1.45e+03 1.45e+03 1.45e+03 2.79e+03 1.52e+03 1.45e+03 1.52e+03 2.24e+03
+DEAL::0 -1.26e+03 -714. -1.26e+03 -32.6 -13.2 -714. -13.2 -32.6
+DEAL::0 -1.26e+03 -714. -1.26e+03 2.57e+03 1.10e+03 -714. 1.10e+03 2.02e+03
+DEAL::0 1.34e+03 -820. 1.34e+03 -32.6 13.2 -820. 13.2 -32.6
+DEAL::0 1.34e+03 -820. 1.34e+03 2.57e+03 -1.31e+03 -820. -1.31e+03 2.24e+03
+DEAL::0 -1.37e+03 1.34e+03 -1.37e+03 -32.6 13.2 1.34e+03 13.2 -32.6
+DEAL::0 -1.37e+03 1.34e+03 -1.37e+03 2.79e+03 -1.31e+03 1.34e+03 -1.31e+03 2.02e+03
+DEAL::0 1.45e+03 1.45e+03 1.45e+03 -32.6 -13.2 1.45e+03 -13.2 -32.6
+DEAL::0 1.45e+03 1.45e+03 1.45e+03 2.79e+03 1.52e+03 1.45e+03 1.52e+03 2.24e+03
DEAL::component=1
-DEAL::-32.6 -1.33e+03 -13.2 -1.33e+03 -2.47e-14 -229. -13.2 -229. -32.6
-DEAL::-32.6 1.41e+03 13.2 1.41e+03 0.00 -335. 13.2 -335. -32.6
+DEAL::-32.6 -1.33e+03 -13.2 -1.33e+03 0 -229. -13.2 -229. -32.6
+DEAL::-32.6 1.41e+03 13.2 1.41e+03 0 -335. 13.2 -335. -32.6
DEAL::2.71e+03 -1.33e+03 1.15e+03 -1.33e+03 1.23e-07 -229. 1.15e+03 -229. 1.05e+03
-DEAL::2.71e+03 1.41e+03 -1.36e+03 1.41e+03 9.06e-08 -335. -1.36e+03 -335. 1.27e+03
-DEAL::-32.6 -1.44e+03 13.2 -1.44e+03 9.40e-15 857. 13.2 857. -32.6
-DEAL::-32.6 1.51e+03 -13.2 1.51e+03 -9.30e-08 963. -13.2 963. -32.6
-DEAL::2.92e+03 -1.44e+03 -1.36e+03 -1.44e+03 8.09e-08 857. -1.36e+03 857. 1.05e+03
-DEAL::2.92e+03 1.51e+03 1.58e+03 1.51e+03 0.00 963. 1.58e+03 963. 1.27e+03
+DEAL::2.71e+03 1.41e+03 -1.36e+03 1.41e+03 0 -335. -1.36e+03 -335. 1.27e+03
+DEAL::-32.6 -1.44e+03 13.2 -1.44e+03 0 857. 13.2 857. -32.6
+DEAL::-32.6 1.51e+03 -13.2 1.51e+03 0 963. -13.2 963. -32.6
+DEAL::2.92e+03 -1.44e+03 -1.36e+03 -1.44e+03 0 857. -1.36e+03 857. 1.05e+03
+DEAL::2.92e+03 1.51e+03 1.58e+03 1.51e+03 0 963. 1.58e+03 963. 1.27e+03
DEAL::component=2
-DEAL::-32.6 -13.2 -845. -13.2 -32.6 -478. -845. -478. -1.94e-08
-DEAL::-32.6 13.2 922. 13.2 -32.6 -584. 922. -584. 4.75e-08
-DEAL::-32.6 13.2 -951. 13.2 -32.6 739. -951. 739. 8.09e-08
-DEAL::-32.6 -13.2 1.03e+03 -13.2 -32.6 845. 1.03e+03 845. -2.94e-08
-DEAL::1.73e+03 1.21e+03 -845. 1.21e+03 1.18e+03 -478. -845. -478. -5.70e-15
-DEAL::1.73e+03 -1.42e+03 922. -1.42e+03 1.40e+03 -584. 922. -584. 3.83e-14
-DEAL::1.95e+03 -1.42e+03 -951. -1.42e+03 1.18e+03 739. -951. 739. 1.33e-14
-DEAL::1.95e+03 1.63e+03 1.03e+03 1.63e+03 1.40e+03 845. 1.03e+03 845. 0.00
+DEAL::-32.6 -13.2 -845. -13.2 -32.6 -478. -845. -478. 0
+DEAL::-32.6 13.2 922. 13.2 -32.6 -584. 922. -584. 0
+DEAL::-32.6 13.2 -951. 13.2 -32.6 739. -951. 739. 0
+DEAL::-32.6 -13.2 1.03e+03 -13.2 -32.6 845. 1.03e+03 845. 0
+DEAL::1.73e+03 1.21e+03 -845. 1.21e+03 1.18e+03 -478. -845. -478. 0
+DEAL::1.73e+03 -1.42e+03 922. -1.42e+03 1.40e+03 -584. 922. -584. 0
+DEAL::1.95e+03 -1.42e+03 -951. -1.42e+03 1.18e+03 739. -951. 739. 0
+DEAL::1.95e+03 1.63e+03 1.03e+03 1.63e+03 1.40e+03 845. 1.03e+03 845. 0
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<2,dim> >(fe.n_components()));
fe_values.get_function_hessians (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
{
FEValuesExtractors::Scalar single_component (c);
fe_values[single_component].get_function_hessians (fe_function,
scalar_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << scalar_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << scalar_values[q][d][e] << (d<dim-1 || e<dim-1 ? " " : "");
+ deallog << std::endl;
+
Assert ((scalar_values[q] - vector_values[q][c]).norm()
<= 1e-12 * scalar_values[q].norm(),
ExcInternalError());
DEAL::-949. 1.90e+03 1.76e+03 -1.73e+03
DEAL::-949. -1.73e+03 -1.39e+03 -877.
DEAL::component=1
-DEAL::2.33e-09 188. 147. -148.
-DEAL::-3.58e-10 188. 147. -204.
-DEAL::-7.59e-09 174. -23.5 33.7
-DEAL::5.84e-10 174. -23.5 -967.
+DEAL::0 188. 147. -148.
+DEAL::0 188. 147. -204.
+DEAL::0 174. -23.5 33.7
+DEAL::0 174. -23.5 -967.
DEAL::FE=FESystem<3>[FE_RaviartThomas<3>(1)]
DEAL::component=0
-DEAL::2.57e-07 -1.50e+03 -752. -1.50e+03 2.23e-07 1.24e+04 -752. 1.24e+04 1.11e-08
-DEAL::2.41e-07 1.50e+03 752. 1.50e+03 3.42e-07 1.32e+04 752. 1.32e+04 4.58e-07
-DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 -1.98e-07 1.24e+04 2.12e+04 1.24e+04 1.83e-06
-DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 -1.68e-07 1.32e+04 -1.94e+04 1.32e+04 -5.46e-07
-DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 -2.08e-08 1.24e+04 -752. 1.24e+04 2.25e-08
-DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 1.23e-06 1.32e+04 752. 1.32e+04 1.80e-13
-DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -5.09e-07 1.24e+04 2.12e+04 1.24e+04 -5.15e-07
-DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 2.10e-08 1.32e+04 -1.94e+04 1.32e+04 -5.46e-07
+DEAL::4.13e-07 -1.50e+03 -752. -1.50e+03 0 1.24e+04 -752. 1.24e+04 -5.75e-07
+DEAL::2.52e-07 1.50e+03 752. 1.50e+03 0 1.32e+04 752. 1.32e+04 -3.15e-07
+DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 5.58e-07 1.24e+04 2.12e+04 1.24e+04 2.05e-06
+DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 5.76e-07 1.32e+04 -1.94e+04 1.32e+04 -5.40e-07
+DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 1.89e-07 1.24e+04 -752. 1.24e+04 3.23e-07
+DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 6.23e-07 1.32e+04 752. 1.32e+04 0
+DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -5.39e-07 1.24e+04 2.12e+04 1.24e+04 -2.63e-07
+DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 0 1.32e+04 -1.94e+04 1.32e+04 -2.54e-07
DEAL::component=1
-DEAL::-1.21e-07 -1.88e+03 1.27e+04 -1.88e+03 1.29e-06 -376. 1.27e+04 -376. 2.69e-07
-DEAL::-2.36e-07 -1.88e+03 1.27e+04 -1.88e+03 3.76e+03 1.63e+04 1.27e+04 1.63e+04 1.84e-06
-DEAL::-1.05e-08 1.88e+03 1.36e+04 1.88e+03 1.03e-06 376. 1.36e+04 376. 4.46e-07
-DEAL::5.04e-08 1.88e+03 1.36e+04 1.88e+03 3.76e+03 -1.46e+04 1.36e+04 -1.46e+04 -5.72e-07
-DEAL::-8.02e-09 1.48e+04 1.27e+04 1.48e+04 752. -376. 1.27e+04 -376. -3.07e-07
-DEAL::-2.14e-08 1.48e+04 1.27e+04 1.48e+04 -2.71e+04 1.63e+04 1.27e+04 1.63e+04 4.27e-13
-DEAL::-9.39e-07 -1.30e+04 1.36e+04 -1.30e+04 752. 376. 1.36e+04 376. -5.04e-08
-DEAL::0.00 -1.30e+04 1.36e+04 -1.30e+04 -2.71e+04 -1.46e+04 1.36e+04 -1.46e+04 0.00
+DEAL::0 -1.88e+03 1.27e+04 -1.88e+03 2.93e-07 -376. 1.27e+04 -376. 2.94e-07
+DEAL::-2.95e-07 -1.88e+03 1.27e+04 -1.88e+03 3.76e+03 1.63e+04 1.27e+04 1.63e+04 2.08e-06
+DEAL::0 1.88e+03 1.36e+04 1.88e+03 -3.29e-07 376. 1.36e+04 376. 0
+DEAL::0 1.88e+03 1.36e+04 1.88e+03 3.76e+03 -1.46e+04 1.36e+04 -1.46e+04 -3.24e-07
+DEAL::0 1.48e+04 1.27e+04 1.48e+04 752. -376. 1.27e+04 -376. -2.73e-07
+DEAL::0 1.48e+04 1.27e+04 1.48e+04 -2.71e+04 1.63e+04 1.27e+04 1.63e+04 0
+DEAL::-1.26e-06 -1.30e+04 1.36e+04 -1.30e+04 752. 376. 1.36e+04 376. 0
+DEAL::0 -1.30e+04 1.36e+04 -1.30e+04 -2.71e+04 -1.46e+04 1.36e+04 -1.46e+04 0
DEAL::component=2
-DEAL::-2.60e-08 1.31e+04 -1.50e+03 1.31e+04 -2.10e-09 -752. -1.50e+03 -752. -1.05e-06
-DEAL::1.26e-08 1.31e+04 -1.50e+03 1.31e+04 1.42e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::1.83e-06 1.31e+04 9.89e+03 1.31e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::-5.72e-07 1.31e+04 9.89e+03 1.31e+04 -2.03e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::-1.10e-07 1.40e+04 1.50e+03 1.40e+04 -7.57e-08 752. 1.50e+03 752. -1.34e-06
-DEAL::0.00 1.40e+04 1.50e+03 1.40e+04 1.20e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::-8.49e-07 1.40e+04 -8.16e+03 1.40e+04 6.72e-08 752. -8.16e+03 752. 1.50e+03
-DEAL::0.00 1.40e+04 -8.16e+03 1.40e+04 3.17e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
+DEAL::0 1.31e+04 -1.50e+03 1.31e+04 -3.05e-07 -752. -1.50e+03 -752. 4.90e-07
+DEAL::0 1.31e+04 -1.50e+03 1.31e+04 1.11e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
+DEAL::2.24e-06 1.31e+04 9.89e+03 1.31e+04 -2.25e-07 -752. 9.89e+03 -752. 1.50e+03
+DEAL::-5.21e-07 1.31e+04 9.89e+03 1.31e+04 -1.67e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
+DEAL::-1.10e-07 1.40e+04 1.50e+03 1.40e+04 -3.26e-07 752. 1.50e+03 752. 1.29e-06
+DEAL::0 1.40e+04 1.50e+03 1.40e+04 1.45e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
+DEAL::-7.99e-07 1.40e+04 -8.16e+03 1.40e+04 -6.22e-07 752. -8.16e+03 752. 1.50e+03
+DEAL::0 1.40e+04 -8.16e+03 1.40e+04 2.21e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<1,dim> >(fe.n_components()));
fe_values.get_function_gradients (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
// use a vector extractor if there
// are sufficiently many components
fe_values[vector_components].get_function_gradients (fe_function,
selected_vector_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
for (unsigned int d=0; d<dim; ++d)
{
- deallog << selected_vector_values[q][d] << std::endl;
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << selected_vector_values[q][d][e] << (e<dim-1 ? " " : "");
+ deallog << std::endl;
Assert ((selected_vector_values[q][d] - vector_values[q][c+d]).norm()
<= 1e-12 * selected_vector_values[q][d].norm(),
ExcInternalError());
DEAL::125. -40.8
DEAL::component=2
DEAL::43.2 -31.7
-DEAL::0.00 5.84
+DEAL::0 5.84
DEAL::43.2 118.
-DEAL::0.00 5.84
+DEAL::0 5.84
DEAL::125. -61.1
-DEAL::0.00 9.52
+DEAL::0 9.52
DEAL::125. -40.8
-DEAL::0.00 9.52
+DEAL::0 9.52
DEAL::component=3
-DEAL::0.00 5.84
+DEAL::0 5.84
DEAL::5.84 -6.75
-DEAL::0.00 5.84
+DEAL::0 5.84
DEAL::5.84 6.75
-DEAL::0.00 9.52
+DEAL::0 9.52
DEAL::9.52 -11.0
-DEAL::0.00 9.52
+DEAL::0 9.52
DEAL::9.52 11.0
DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(0)]
DEAL::component=0
DEAL::2.37 4.73 9.46
-DEAL::-2.21e-12 -275. -138.
-DEAL::-275. -8.99e-13 -138.
+DEAL::0 -275. -138.
+DEAL::-275. 0 -138.
DEAL::2.37 4.73 9.46
-DEAL::2.77e-12 -275. -138.
+DEAL::0 -275. -138.
DEAL::-275. -459. 4.80e+03
DEAL::2.37 4.73 9.46
DEAL::-367. -275. 4.71e+03
-DEAL::-275. 1.41e-12 -138.
+DEAL::-275. 0 -138.
DEAL::2.37 4.73 9.46
DEAL::367. -275. 4.92e+03
DEAL::-275. 459. 5.02e+03
DEAL::-4.38e+03 4.79e+03 4.92e+03
DEAL::4.88e+03 -3.09e+03 5.02e+03
DEAL::component=1
-DEAL::-2.21e-12 -275. -138.
-DEAL::-275. -8.99e-13 -138.
-DEAL::-275. -138. 1.36e-12
-DEAL::2.77e-12 -275. -138.
+DEAL::0 -275. -138.
+DEAL::-275. 0 -138.
+DEAL::-275. -138. 0
+DEAL::0 -275. -138.
DEAL::-275. -459. 4.80e+03
DEAL::-275. 4.90e+03 -367.
DEAL::-367. -275. 4.71e+03
-DEAL::-275. 1.41e-12 -138.
+DEAL::-275. 0 -138.
DEAL::4.76e+03 -138. -184.
DEAL::367. -275. 4.92e+03
DEAL::-275. 459. 5.02e+03
DEAL::4.76e+03 4.90e+03 2.23e+03
DEAL::-184. 4.57e+03 -138.
DEAL::4.67e+03 -91.8 -138.
-DEAL::-275. -138. -1.36e-12
+DEAL::-275. -138. 0
DEAL::184. 4.79e+03 -138.
DEAL::4.67e+03 3.52e+03 4.80e+03
DEAL::-275. 5.11e+03 367.
DEAL::4.88e+03 -3.09e+03 5.02e+03
DEAL::4.97e+03 5.11e+03 -1.81e+03
DEAL::component=2
-DEAL::-275. -8.99e-13 -138.
-DEAL::-275. -138. 1.36e-12
-DEAL::0.00 5.60 22.4
+DEAL::-275. 0 -138.
+DEAL::-275. -138. 0
+DEAL::0 5.60 22.4
DEAL::-275. -459. 4.80e+03
DEAL::-275. 4.90e+03 -367.
-DEAL::0.00 5.60 22.4
-DEAL::-275. 1.41e-12 -138.
+DEAL::0 5.60 22.4
+DEAL::-275. 0 -138.
DEAL::4.76e+03 -138. -184.
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::-275. 459. 5.02e+03
DEAL::4.76e+03 4.90e+03 2.23e+03
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::4.67e+03 -91.8 -138.
-DEAL::-275. -138. -1.36e-12
-DEAL::0.00 5.60 22.4
+DEAL::-275. -138. 0
+DEAL::0 5.60 22.4
DEAL::4.67e+03 3.52e+03 4.80e+03
DEAL::-275. 5.11e+03 367.
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::4.88e+03 91.8 -138.
DEAL::4.97e+03 -138. 184.
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::4.88e+03 -3.09e+03 5.02e+03
DEAL::4.97e+03 5.11e+03 -1.81e+03
-DEAL::0.00 5.60 22.4
+DEAL::0 5.60 22.4
DEAL::component=3
-DEAL::-275. -138. 1.36e-12
-DEAL::0.00 5.60 22.4
-DEAL::5.60 -1.84e-16 22.4
+DEAL::-275. -138. 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::-275. 4.90e+03 -367.
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::4.76e+03 -138. -184.
-DEAL::0.00 5.60 22.4
-DEAL::5.60 -1.84e-16 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::4.76e+03 4.90e+03 2.23e+03
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
-DEAL::-275. -138. -1.36e-12
-DEAL::0.00 5.60 22.4
-DEAL::5.60 7.35e-16 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::-275. -138. 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::-275. 5.11e+03 367.
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::4.97e+03 -138. 184.
-DEAL::0.00 5.60 22.4
-DEAL::5.60 7.35e-16 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::4.97e+03 5.11e+03 -1.81e+03
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
DEAL::component=4
-DEAL::0.00 5.60 22.4
-DEAL::5.60 -1.84e-16 22.4
-DEAL::5.60 11.2 2.53e-16
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 11.2 1.01e-15
-DEAL::0.00 5.60 22.4
-DEAL::5.60 -1.84e-16 22.4
-DEAL::5.60 11.2 1.01e-15
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 11.2 0.00
-DEAL::0.00 5.60 22.4
-DEAL::5.60 7.35e-16 22.4
-DEAL::5.60 11.2 2.53e-16
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 11.2 1.01e-15
-DEAL::0.00 5.60 22.4
-DEAL::5.60 7.35e-16 22.4
-DEAL::5.60 11.2 1.01e-15
-DEAL::0.00 5.60 22.4
-DEAL::5.60 0.00 22.4
-DEAL::5.60 11.2 0.00
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
+DEAL::0 5.60 22.4
+DEAL::5.60 0 22.4
+DEAL::5.60 11.2 0
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<2,dim> >(fe.n_components()));
fe_values.get_function_hessians (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
// use a vector extractor if there
// are sufficiently many components
fe_values[vector_components].get_function_hessians (fe_function,
selected_vector_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
for (unsigned int d=0; d<dim; ++d)
{
- deallog << selected_vector_values[q][d] << std::endl;
+ for (unsigned int e=0; e<dim; ++e)
+ for (unsigned int f=0; f<dim; ++f)
+ deallog << selected_vector_values[q][d][e][f] << (e<dim-1 && f<dim-1 ? " " : "");
+ deallog << std::endl;
Assert ((selected_vector_values[q][d] - vector_values[q][c+d]).norm()
<= 1e-12 * selected_vector_values[q][d].norm(),
ExcInternalError());
deallog.attach(logfile);
deallog.depth_console (0);
- deallog.threshold_double(1.e-7);
+ deallog.threshold_double(1.e-6);
test_hyper_sphere<2>();
test_hyper_sphere<3>();
DEAL::FE=FESystem<2>[FE_Q<2>(1)-FE_RaviartThomas<2>(1)-FE_Nedelec<2>(0)]
DEAL::component=0
-DEAL::5.60e-11 1.30 1.30 -1.50
-DEAL::-50.8 877. 848. -953.
-DEAL::0.00 1.30 1.30 1.50
-DEAL::-50.8 -779. -889. -641.
-DEAL::0.00 3.45 3.45 -3.99
-DEAL::-1.43e+03 2.86e+03 2.65e+03 -2.63e+03
-DEAL::0.00 3.45 3.45 3.99
-DEAL::-1.43e+03 -2.60e+03 -2.12e+03 -1.44e+03
+DEAL::0 1.301.30-1.50
+DEAL::-50.8 877.848.-953.
+DEAL::0 1.301.301.50
+DEAL::-50.8 -779.-889.-641.
+DEAL::0 3.453.45-3.99
+DEAL::-1.43e+03 2.86e+032.65e+03-2.63e+03
+DEAL::0 3.453.453.99
+DEAL::-1.43e+03 -2.60e+03-2.12e+03-1.44e+03
DEAL::component=1
-DEAL::-50.8 877. 848. -953.
-DEAL::3.99e-09 278. 209. -236.
-DEAL::-50.8 -779. -889. -641.
-DEAL::-6.72e-10 278. 209. -247.
-DEAL::-1.43e+03 2.86e+03 2.65e+03 -2.63e+03
-DEAL::-1.25e-08 374. 46.3 -122.
-DEAL::-1.43e+03 -2.60e+03 -2.12e+03 -1.44e+03
-DEAL::8.76e-10 374. 46.3 -1.16e+03
+DEAL::-50.8 877.848.-953.
+DEAL::0 278.209.-236.
+DEAL::-50.8 -779.-889.-641.
+DEAL::0 278.209.-247.
+DEAL::-1.43e+03 2.86e+032.65e+03-2.63e+03
+DEAL::0 374.46.3-122.
+DEAL::-1.43e+03 -2.60e+03-2.12e+03-1.44e+03
+DEAL::0 374.46.3-1.16e+03
DEAL::component=2
-DEAL::3.99e-09 278. 209. -236.
-DEAL::0.00 0.00 2.69e-10 9.43
-DEAL::-6.72e-10 278. 209. -247.
-DEAL::0.00 0.00 8.06e-10 9.43
-DEAL::-1.25e-08 374. 46.3 -122.
-DEAL::0.00 0.00 -8.76e-10 25.0
-DEAL::8.76e-10 374. 46.3 -1.16e+03
-DEAL::0.00 0.00 8.76e-10 25.0
+DEAL::0 278.209.-236.
+DEAL::0 009.43
+DEAL::0 278.209.-247.
+DEAL::0 009.43
+DEAL::0 374.46.3-122.
+DEAL::0 0025.0
+DEAL::0 374.46.3-1.16e+03
+DEAL::0 0025.0
DEAL::component=3
-DEAL::0.00 0.00 2.69e-10 9.43
-DEAL::2.69e-10 9.43 18.9 -21.8
-DEAL::0.00 0.00 8.06e-10 9.43
-DEAL::8.06e-10 9.43 18.9 21.8
-DEAL::0.00 0.00 -8.76e-10 25.0
-DEAL::-8.76e-10 25.0 50.0 -57.8
-DEAL::0.00 0.00 8.76e-10 25.0
-DEAL::8.76e-10 25.0 50.0 57.8
+DEAL::0 009.43
+DEAL::0 9.4318.9-21.8
+DEAL::0 009.43
+DEAL::0 9.4318.921.8
+DEAL::0 0025.0
+DEAL::0 25.050.0-57.8
+DEAL::0 0025.0
+DEAL::0 25.050.057.8
DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(0)]
DEAL::component=0
-DEAL::0.00 9.85e-10 5.09e-10 1.54e-09 2.44e-26 5.91e-10 3.28e-10 4.60e-10 -2.17e-09
-DEAL::2.57e-07 -1.50e+03 -752. -1.50e+03 3.49e-07 1.99e+04 -752. 1.99e+04 5.31e-08
-DEAL::-2.06e-07 -1.88e+03 2.02e+04 -1.88e+03 1.75e-06 -376. 2.02e+04 -376. 4.37e-07
-DEAL::0.00 9.85e-10 1.74e-09 -1.38e-09 0.00 2.53e-09 -1.30e-09 5.58e-10 1.64e-09
-DEAL::-3.47e-07 1.50e+03 752. 1.50e+03 5.81e-07 2.07e+04 752. 2.07e+04 7.10e-07
-DEAL::-3.62e-07 -1.88e+03 2.02e+04 -1.88e+03 3.76e+03 1.63e+04 2.02e+04 1.63e+04 2.93e-06
-DEAL::4.27e-10 -1.18e-09 1.26e-09 2.89e-09 -8.21e-10 2.89e-09 -8.54e-10 -7.88e-10 1.64e-09
-DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 -4.08e-07 1.99e+04 2.12e+04 1.99e+04 2.96e-06
-DEAL::1.88e-08 1.88e+03 2.11e+04 1.88e+03 1.28e-06 376. 2.11e+04 376. 7.40e-07
-DEAL::0.00 -1.18e-09 -2.30e-10 -3.94e-10 -8.21e-10 -4.60e-10 -7.88e-10 -7.88e-10 -2.17e-09
-DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 -2.94e-07 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
-DEAL::5.04e-08 1.88e+03 2.11e+04 1.88e+03 3.76e+03 -1.46e+04 2.11e+04 -1.46e+04 -9.08e-07
-DEAL::0.00 5.09e-10 -1.97e-10 1.97e-10 -1.94e-26 -1.38e-09 3.28e-10 4.60e-10 3.23e-27
-DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 -6.29e-08 1.99e+04 -752. 1.99e+04 2.25e-08
-DEAL::-8.02e-09 1.48e+04 2.02e+04 1.48e+04 752. -376. 2.02e+04 -376. -5.17e-07
-DEAL::0.00 5.09e-10 -2.13e-10 -3.94e-10 1.08e-09 -3.94e-10 -1.63e-09 9.52e-10 -1.11e-26
-DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 2.06e-06 2.07e+04 752. 2.07e+04 3.03e-13
-DEAL::-2.14e-08 1.48e+04 2.02e+04 1.48e+04 -2.71e+04 1.63e+04 2.02e+04 1.63e+04 5.50e-13
-DEAL::-5.42e-10 4.92e-11 5.42e-10 1.97e-10 -8.21e-10 -3.94e-10 -2.36e-09 -7.88e-10 -5.18e-26
-DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -8.04e-07 1.99e+04 2.12e+04 1.99e+04 -7.68e-07
-DEAL::-1.46e-06 -1.30e+04 2.11e+04 -1.30e+04 752. 376. 2.11e+04 376. -9.25e-08
-DEAL::0.00 4.92e-11 4.93e-11 -3.94e-10 0.00 -1.38e-09 -7.88e-10 -7.88e-10 0.00
-DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 2.10e-08 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
-DEAL::0.00 -1.30e+04 2.11e+04 -1.30e+04 -2.71e+04 -1.46e+04 2.11e+04 -1.46e+04 0.00
+DEAL::0 0 00 0 0000
+DEAL::0 -1.50e+03 -752.-1.50e+03 0 1.99e+04-752.1.99e+040
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 0 -376.2.02e+04-376.0
+DEAL::0 0 00 0 0000
+DEAL::0 1.50e+03 752.1.50e+03 0 2.07e+04752.2.07e+040
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 3.76e+03 1.63e+042.02e+041.63e+043.26e-06
+DEAL::0 0 00 0 0000
+DEAL::3.01e+03 -1.50e+03 2.12e+04-1.50e+03 0 1.99e+042.12e+041.99e+043.27e-06
+DEAL::0 1.88e+03 2.11e+041.88e+03 0 376.2.11e+04376.0
+DEAL::0 0 00 0 0000
+DEAL::3.01e+03 1.50e+03 -1.94e+041.50e+03 0 2.07e+04-1.94e+042.07e+040
+DEAL::0 1.88e+03 2.11e+041.88e+03 3.76e+03 -1.46e+042.11e+04-1.46e+040
+DEAL::0 0 00 0 0000
+DEAL::1.50e+03 2.04e+04 -752.2.04e+04 0 1.99e+04-752.1.99e+040
+DEAL::0 1.48e+04 2.02e+041.48e+04 752. -376.2.02e+04-376.0
+DEAL::0 0 00 0 0000
+DEAL::1.50e+03 -1.87e+04 752.-1.87e+04 0 2.07e+04752.2.07e+040
+DEAL::0 1.48e+04 2.02e+041.48e+04 -2.71e+04 1.63e+042.02e+041.63e+040
+DEAL::0 0 00 0 0000
+DEAL::-3.76e+04 2.04e+04 2.12e+042.04e+04 0 1.99e+042.12e+041.99e+040
+DEAL::-1.98e-06 -1.30e+04 2.11e+04-1.30e+04 752. 376.2.11e+04376.0
+DEAL::0 0 00 0 0000
+DEAL::-3.76e+04 -1.87e+04 -1.94e+04-1.87e+04 0 2.07e+04-1.94e+042.07e+040
+DEAL::0 -1.30e+04 2.11e+04-1.30e+04 -2.71e+04 -1.46e+042.11e+04-1.46e+040
DEAL::component=1
-DEAL::2.57e-07 -1.50e+03 -752. -1.50e+03 3.49e-07 1.99e+04 -752. 1.99e+04 5.31e-08
-DEAL::-2.06e-07 -1.88e+03 2.02e+04 -1.88e+03 1.75e-06 -376. 2.02e+04 -376. 4.37e-07
-DEAL::5.81e-08 2.06e+04 -1.50e+03 2.06e+04 3.99e-08 -752. -1.50e+03 -752. -1.93e-06
-DEAL::-3.47e-07 1.50e+03 752. 1.50e+03 5.81e-07 2.07e+04 752. 2.07e+04 7.10e-07
-DEAL::-3.62e-07 -1.88e+03 2.02e+04 -1.88e+03 3.76e+03 1.63e+04 2.02e+04 1.63e+04 2.93e-06
-DEAL::-2.94e-08 2.06e+04 -1.50e+03 2.06e+04 2.26e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::3.01e+03 -1.50e+03 2.12e+04 -1.50e+03 -4.08e-07 1.99e+04 2.12e+04 1.99e+04 2.96e-06
-DEAL::1.88e-08 1.88e+03 2.11e+04 1.88e+03 1.28e-06 376. 2.11e+04 376. 7.40e-07
-DEAL::3.01e-06 2.06e+04 9.89e+03 2.06e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::3.01e+03 1.50e+03 -1.94e+04 1.50e+03 -2.94e-07 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
-DEAL::5.04e-08 1.88e+03 2.11e+04 1.88e+03 3.76e+03 -1.46e+04 2.11e+04 -1.46e+04 -9.08e-07
-DEAL::-9.50e-07 2.06e+04 9.89e+03 2.06e+04 -3.25e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::1.50e+03 2.04e+04 -752. 2.04e+04 -6.29e-08 1.99e+04 -752. 1.99e+04 2.25e-08
-DEAL::-8.02e-09 1.48e+04 2.02e+04 1.48e+04 752. -376. 2.02e+04 -376. -5.17e-07
-DEAL::-1.81e-07 2.15e+04 1.50e+03 2.15e+04 -1.18e-07 752. 1.50e+03 752. -2.35e-06
-DEAL::1.50e+03 -1.87e+04 752. -1.87e+04 2.06e-06 2.07e+04 752. 2.07e+04 3.03e-13
-DEAL::-2.14e-08 1.48e+04 2.02e+04 1.48e+04 -2.71e+04 1.63e+04 2.02e+04 1.63e+04 5.50e-13
-DEAL::0.00 2.15e+04 1.50e+03 2.15e+04 1.87e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::-3.76e+04 2.04e+04 2.12e+04 2.04e+04 -8.04e-07 1.99e+04 2.12e+04 1.99e+04 -7.68e-07
-DEAL::-1.46e-06 -1.30e+04 2.11e+04 -1.30e+04 752. 376. 2.11e+04 376. -9.25e-08
-DEAL::-1.34e-06 2.15e+04 -8.16e+03 2.15e+04 1.09e-07 752. -8.16e+03 752. 1.50e+03
-DEAL::-3.76e+04 -1.87e+04 -1.94e+04 -1.87e+04 2.10e-08 2.07e+04 -1.94e+04 2.07e+04 -8.83e-07
-DEAL::0.00 -1.30e+04 2.11e+04 -1.30e+04 -2.71e+04 -1.46e+04 2.11e+04 -1.46e+04 0.00
-DEAL::0.00 2.15e+04 -8.16e+03 2.15e+04 5.27e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
+DEAL::0 -1.50e+03 -752.-1.50e+03 0 1.99e+04-752.1.99e+040
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 0 -376.2.02e+04-376.0
+DEAL::0 2.06e+04 -1.50e+032.06e+04 0 -752.-1.50e+03-752.0
+DEAL::0 1.50e+03 752.1.50e+03 0 2.07e+04752.2.07e+040
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 3.76e+03 1.63e+042.02e+041.63e+043.26e-06
+DEAL::0 2.06e+04 -1.50e+032.06e+04 1.78e-06 1.06e+04-1.50e+031.06e+043.01e+03
+DEAL::3.01e+03 -1.50e+03 2.12e+04-1.50e+03 0 1.99e+042.12e+041.99e+043.27e-06
+DEAL::0 1.88e+03 2.11e+041.88e+03 0 376.2.11e+04376.0
+DEAL::3.66e-06 2.06e+04 9.89e+032.06e+04 0 -752.9.89e+03-752.1.50e+03
+DEAL::3.01e+03 1.50e+03 -1.94e+041.50e+03 0 2.07e+04-1.94e+042.07e+040
+DEAL::0 1.88e+03 2.11e+041.88e+03 3.76e+03 -1.46e+042.11e+04-1.46e+040
+DEAL::0 2.06e+04 9.89e+032.06e+04 -2.68e-06 1.06e+049.89e+031.06e+04-1.65e+04
+DEAL::1.50e+03 2.04e+04 -752.2.04e+04 0 1.99e+04-752.1.99e+040
+DEAL::0 1.48e+04 2.02e+041.48e+04 752. -376.2.02e+04-376.0
+DEAL::0 2.15e+04 1.50e+032.15e+04 0 752.1.50e+03752.1.88e-06
+DEAL::1.50e+03 -1.87e+04 752.-1.87e+04 0 2.07e+04752.2.07e+040
+DEAL::0 1.48e+04 2.02e+041.48e+04 -2.71e+04 1.63e+042.02e+041.63e+040
+DEAL::0 2.15e+04 1.50e+032.15e+04 2.25e-06 -8.91e+031.50e+03-8.91e+033.01e+03
+DEAL::-3.76e+04 2.04e+04 2.12e+042.04e+04 0 1.99e+042.12e+041.99e+040
+DEAL::-1.98e-06 -1.30e+04 2.11e+04-1.30e+04 752. 376.2.11e+04376.0
+DEAL::-1.25e-06 2.15e+04 -8.16e+032.15e+04 -1.00e-06 752.-8.16e+03752.1.50e+03
+DEAL::-3.76e+04 -1.87e+04 -1.94e+04-1.87e+04 0 2.07e+04-1.94e+042.07e+040
+DEAL::0 -1.30e+04 2.11e+04-1.30e+04 -2.71e+04 -1.46e+042.11e+04-1.46e+040
+DEAL::0 2.15e+04 -8.16e+032.15e+04 0 -8.91e+03-8.16e+03-8.91e+03-1.65e+04
DEAL::component=2
-DEAL::-2.06e-07 -1.88e+03 2.02e+04 -1.88e+03 1.75e-06 -376. 2.02e+04 -376. 4.37e-07
-DEAL::5.81e-08 2.06e+04 -1.50e+03 2.06e+04 3.99e-08 -752. -1.50e+03 -752. -1.93e-06
-DEAL::0.00 0.00 0.00 1.58e-09 2.63e-09 2.36e-09 0.00 1.05e-08 4.31e-26
-DEAL::-3.62e-07 -1.88e+03 2.02e+04 -1.88e+03 3.76e+03 1.63e+04 2.02e+04 1.63e+04 2.93e-06
-DEAL::-2.94e-08 2.06e+04 -1.50e+03 2.06e+04 2.26e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::0.00 0.00 0.00 -2.63e-10 2.63e-09 6.30e-09 -6.30e-09 1.37e-08 2.10e-08
-DEAL::1.88e-08 1.88e+03 2.11e+04 1.88e+03 1.28e-06 376. 2.11e+04 376. 7.40e-07
-DEAL::3.01e-06 2.06e+04 9.89e+03 2.06e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::0.00 0.00 0.00 5.25e-09 -3.94e-09 6.04e-09 8.14e-09 -7.62e-09 1.58e-08
-DEAL::5.04e-08 1.88e+03 2.11e+04 1.88e+03 3.76e+03 -1.46e+04 2.11e+04 -1.46e+04 -9.08e-07
-DEAL::-9.50e-07 2.06e+04 9.89e+03 2.06e+04 -3.25e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::0.00 0.00 0.00 0.00 -3.94e-09 0.00 0.00 -1.05e-08 -5.25e-09
-DEAL::-8.02e-09 1.48e+04 2.02e+04 1.48e+04 752. -376. 2.02e+04 -376. -5.17e-07
-DEAL::-1.81e-07 2.15e+04 1.50e+03 2.15e+04 -1.18e-07 752. 1.50e+03 752. -2.35e-06
-DEAL::0.00 0.00 0.00 0.00 1.31e-09 -1.58e-09 0.00 4.20e-09 1.73e-26
-DEAL::-2.14e-08 1.48e+04 2.02e+04 1.48e+04 -2.71e+04 1.63e+04 2.02e+04 1.63e+04 5.50e-13
-DEAL::0.00 2.15e+04 1.50e+03 2.15e+04 1.87e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::0.00 0.00 0.00 0.00 1.58e-09 2.63e-10 -8.41e-09 8.41e-09 -3.94e-31
-DEAL::-1.46e-06 -1.30e+04 2.11e+04 -1.30e+04 752. 376. 2.11e+04 376. -9.25e-08
-DEAL::-1.34e-06 2.15e+04 -8.16e+03 2.15e+04 1.09e-07 752. -8.16e+03 752. 1.50e+03
-DEAL::0.00 0.00 0.00 -1.31e-09 -1.31e-09 1.31e-09 -6.57e-09 3.94e-09 5.25e-09
-DEAL::0.00 -1.30e+04 2.11e+04 -1.30e+04 -2.71e+04 -1.46e+04 2.11e+04 -1.46e+04 0.00
-DEAL::0.00 2.15e+04 -8.16e+03 2.15e+04 5.27e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
-DEAL::0.00 0.00 0.00 0.00 0.00 -1.31e-09 0.00 0.00 0.00
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 0 -376.2.02e+04-376.0
+DEAL::0 2.06e+04 -1.50e+032.06e+04 0 -752.-1.50e+03-752.0
+DEAL::0 0 00 0 0000
+DEAL::0 -1.88e+03 2.02e+04-1.88e+03 3.76e+03 1.63e+042.02e+041.63e+043.26e-06
+DEAL::0 2.06e+04 -1.50e+032.06e+04 1.78e-06 1.06e+04-1.50e+031.06e+043.01e+03
+DEAL::0 0 00 0 0000
+DEAL::0 1.88e+03 2.11e+041.88e+03 0 376.2.11e+04376.0
+DEAL::3.66e-06 2.06e+04 9.89e+032.06e+04 0 -752.9.89e+03-752.1.50e+03
+DEAL::0 0 00 0 0000
+DEAL::0 1.88e+03 2.11e+041.88e+03 3.76e+03 -1.46e+042.11e+04-1.46e+040
+DEAL::0 2.06e+04 9.89e+032.06e+04 -2.68e-06 1.06e+049.89e+031.06e+04-1.65e+04
+DEAL::0 0 00 0 0000
+DEAL::0 1.48e+04 2.02e+041.48e+04 752. -376.2.02e+04-376.0
+DEAL::0 2.15e+04 1.50e+032.15e+04 0 752.1.50e+03752.1.88e-06
+DEAL::0 0 00 0 0000
+DEAL::0 1.48e+04 2.02e+041.48e+04 -2.71e+04 1.63e+042.02e+041.63e+040
+DEAL::0 2.15e+04 1.50e+032.15e+04 2.25e-06 -8.91e+031.50e+03-8.91e+033.01e+03
+DEAL::0 0 00 0 0000
+DEAL::-1.98e-06 -1.30e+04 2.11e+04-1.30e+04 752. 376.2.11e+04376.0
+DEAL::-1.25e-06 2.15e+04 -8.16e+032.15e+04 -1.00e-06 752.-8.16e+03752.1.50e+03
+DEAL::0 0 00 0 0000
+DEAL::0 -1.30e+04 2.11e+04-1.30e+04 -2.71e+04 -1.46e+042.11e+04-1.46e+040
+DEAL::0 2.15e+04 -8.16e+032.15e+04 0 -8.91e+03-8.16e+03-8.91e+03-1.65e+04
+DEAL::0 0 00 0 0000
DEAL::component=3
-DEAL::5.81e-08 2.06e+04 -1.50e+03 2.06e+04 3.99e-08 -752. -1.50e+03 -752. -1.93e-06
-DEAL::0.00 0.00 0.00 1.58e-09 2.63e-09 2.36e-09 0.00 1.05e-08 4.31e-26
-DEAL::2.35e-09 1.58e-09 2.20e-09 -1.20e-09 4.32e-26 -7.61e-10 1.14e-08 4.33e-25 4.35e-10
-DEAL::-2.94e-08 2.06e+04 -1.50e+03 2.06e+04 2.26e-06 1.06e+04 -1.50e+03 1.06e+04 3.01e+03
-DEAL::0.00 0.00 0.00 -2.63e-10 2.63e-09 6.30e-09 -6.30e-09 1.37e-08 2.10e-08
-DEAL::-1.31e-09 1.58e-09 6.15e-09 0.00 0.00 6.52e-10 -2.89e-09 0.00 1.58e-08
-DEAL::3.01e-06 2.06e+04 9.89e+03 2.06e+04 7.78e-08 -752. 9.89e+03 -752. 1.50e+03
-DEAL::0.00 0.00 0.00 5.25e-09 -3.94e-09 6.04e-09 8.14e-09 -7.62e-09 1.58e-08
-DEAL::4.98e-09 -2.89e-09 6.14e-09 -1.20e-09 -7.35e-26 -7.61e-10 1.82e-08 -1.05e-08 2.14e-08
-DEAL::-9.50e-07 2.06e+04 9.89e+03 2.06e+04 -3.25e-06 1.06e+04 9.89e+03 1.06e+04 -1.65e+04
-DEAL::0.00 0.00 0.00 0.00 -3.94e-09 0.00 0.00 -1.05e-08 -5.25e-09
-DEAL::0.00 -2.89e-09 1.09e-10 0.00 0.00 6.52e-10 -1.05e-09 -5.25e-09 -5.25e-09
-DEAL::-1.81e-07 2.15e+04 1.50e+03 2.15e+04 -1.18e-07 752. 1.50e+03 752. -2.35e-06
-DEAL::0.00 0.00 0.00 0.00 1.31e-09 -1.58e-09 0.00 4.20e-09 1.73e-26
-DEAL::2.17e-10 1.31e-09 -1.68e-09 1.74e-09 9.45e-26 -8.70e-10 -1.81e-10 -6.81e-26 -5.69e-09
-DEAL::0.00 2.15e+04 1.50e+03 2.15e+04 1.87e-06 -8.91e+03 1.50e+03 -8.91e+03 3.01e+03
-DEAL::0.00 0.00 0.00 0.00 1.58e-09 2.63e-10 -8.41e-09 8.41e-09 -3.94e-31
-DEAL::0.00 2.63e-09 -2.63e-10 0.00 0.00 0.00 2.36e-09 5.25e-09 1.25e-25
-DEAL::-1.34e-06 2.15e+04 -8.16e+03 2.15e+04 1.09e-07 752. -8.16e+03 752. 1.50e+03
-DEAL::0.00 0.00 0.00 -1.31e-09 -1.31e-09 1.31e-09 -6.57e-09 3.94e-09 5.25e-09
-DEAL::-4.52e-11 -1.31e-09 9.42e-10 1.74e-09 6.17e-26 -8.70e-10 -1.76e-09 -1.26e-08 -4.35e-10
-DEAL::0.00 2.15e+04 -8.16e+03 2.15e+04 5.27e-07 -8.91e+03 -8.16e+03 -8.91e+03 -1.65e+04
-DEAL::0.00 0.00 0.00 0.00 0.00 -1.31e-09 0.00 0.00 0.00
-DEAL::0.00 0.00 -1.58e-09 0.00 0.00 0.00 -1.05e-09 0.00 -5.25e-09
+DEAL::0 2.06e+04 -1.50e+032.06e+04 0 -752.-1.50e+03-752.0
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 2.06e+04 -1.50e+032.06e+04 1.78e-06 1.06e+04-1.50e+031.06e+043.01e+03
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::3.66e-06 2.06e+04 9.89e+032.06e+04 0 -752.9.89e+03-752.1.50e+03
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 2.06e+04 9.89e+032.06e+04 -2.68e-06 1.06e+049.89e+031.06e+04-1.65e+04
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 2.15e+04 1.50e+032.15e+04 0 752.1.50e+03752.1.88e-06
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 2.15e+04 1.50e+032.15e+04 2.25e-06 -8.91e+031.50e+03-8.91e+033.01e+03
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::-1.25e-06 2.15e+04 -8.16e+032.15e+04 -1.00e-06 752.-8.16e+03752.1.50e+03
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 2.15e+04 -8.16e+032.15e+04 0 -8.91e+03-8.16e+03-8.91e+03-1.65e+04
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
DEAL::component=4
-DEAL::0.00 0.00 0.00 1.58e-09 2.63e-09 2.36e-09 0.00 1.05e-08 4.31e-26
-DEAL::2.35e-09 1.58e-09 2.20e-09 -1.20e-09 4.32e-26 -7.61e-10 1.14e-08 4.33e-25 4.35e-10
-DEAL::-1.63e-10 2.42e-09 2.38e-26 4.78e-09 -1.63e-10 1.76e-26 -7.07e-10 -4.89e-10 3.15e-26
-DEAL::0.00 0.00 0.00 -2.63e-10 2.63e-09 6.30e-09 -6.30e-09 1.37e-08 2.10e-08
-DEAL::-1.31e-09 1.58e-09 6.15e-09 0.00 0.00 6.52e-10 -2.89e-09 0.00 1.58e-08
-DEAL::-4.35e-10 2.58e-09 4.20e-09 -3.37e-09 -2.17e-10 7.88e-09 -3.70e-09 2.17e-10 6.90e-25
-DEAL::0.00 0.00 0.00 5.25e-09 -3.94e-09 6.04e-09 8.14e-09 -7.62e-09 1.58e-08
-DEAL::4.98e-09 -2.89e-09 6.14e-09 -1.20e-09 -7.35e-26 -7.61e-10 1.82e-08 -1.05e-08 2.14e-08
-DEAL::1.10e-09 -3.43e-09 3.94e-09 8.10e-09 -3.17e-09 8.41e-09 -2.17e-10 -5.98e-09 6.64e-25
-DEAL::0.00 0.00 0.00 0.00 -3.94e-09 0.00 0.00 -1.05e-08 -5.25e-09
-DEAL::0.00 -2.89e-09 1.09e-10 0.00 0.00 6.52e-10 -1.05e-09 -5.25e-09 -5.25e-09
-DEAL::0.00 -2.89e-09 -1.31e-09 -5.25e-10 -2.63e-09 -2.63e-09 0.00 0.00 0.00
-DEAL::0.00 0.00 0.00 0.00 1.31e-09 -1.58e-09 0.00 4.20e-09 1.73e-26
-DEAL::2.17e-10 1.31e-09 -1.68e-09 1.74e-09 9.45e-26 -8.70e-10 -1.81e-10 -6.81e-26 -5.69e-09
-DEAL::-1.63e-10 1.11e-09 -2.51e-27 -4.71e-10 -1.63e-10 -2.63e-09 -7.07e-10 -4.89e-10 -5.96e-26
-DEAL::0.00 0.00 0.00 0.00 1.58e-09 2.63e-10 -8.41e-09 8.41e-09 -3.94e-31
-DEAL::0.00 2.63e-09 -2.63e-10 0.00 0.00 0.00 2.36e-09 5.25e-09 1.25e-25
-DEAL::-4.35e-10 1.27e-09 1.00e-26 -7.43e-10 3.08e-10 -1.08e-26 -3.70e-09 2.17e-10 -5.71e-26
-DEAL::0.00 0.00 0.00 -1.31e-09 -1.31e-09 1.31e-09 -6.57e-09 3.94e-09 5.25e-09
-DEAL::-4.52e-11 -1.31e-09 9.42e-10 1.74e-09 6.17e-26 -8.70e-10 -1.76e-09 -1.26e-08 -4.35e-10
-DEAL::-1.79e-09 -8.06e-10 1.31e-09 -2.41e-09 -3.17e-09 -7.87e-26 -2.17e-10 -5.98e-09 5.87e-26
-DEAL::0.00 0.00 0.00 0.00 0.00 -1.31e-09 0.00 0.00 0.00
-DEAL::0.00 0.00 -1.58e-09 0.00 0.00 0.00 -1.05e-09 0.00 -5.25e-09
-DEAL::0.00 -2.63e-10 0.00 -5.25e-10 0.00 -2.63e-09 0.00 0.00 0.00
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
+DEAL::0 0 00 0 0000
//----------------------------------------------------------------------
// $Id$
-// Version: $Name$
+// Version: $Name$
//
-// Copyright (C) 2007, 2008, 2010 by the deal.II authors
+// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
Vector<double> fe_function(dof.n_dofs());
for (unsigned int i=0; i<dof.n_dofs(); ++i)
fe_function(i) = i+1;
-
+
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients | update_hessians);
std::vector<Tensor<1,dim> >(fe.n_components()));
fe_values.get_function_gradients (fe_function, vector_values);
-
+
for (unsigned int c=0; c<fe.n_components(); ++c)
// use a vector extractor if there
// are sufficiently many components
fe_values[vector_components].get_function_symmetric_gradients (fe_function,
selected_vector_values);
deallog << "component=" << c << std::endl;
-
+
for (unsigned int q=0; q<fe_values.n_quadrature_points; ++q)
{
- deallog << selected_vector_values[q] << std::endl;
+ for (unsigned int d=0; d<dim; ++d)
+ for (unsigned int e=0; e<dim; ++e)
+ deallog << selected_vector_values[q][d][e] << (d<dim-1 || e<dim-1 ? " " : "");
+ deallog << std::endl;
Tensor<2,dim> grad;
for (unsigned int d=0; d<dim; ++d)
grad[d] = vector_values[q][c+d];
-
+
Assert ((selected_vector_values[q] - symmetrize(grad)).norm()
<= 1e-12 * selected_vector_values[q].norm(),
ExcInternalError());
DEAL::125. -30.5 -30.5 9.52
DEAL::125. -20.4 -20.4 9.52
DEAL::component=3
-DEAL::0.00 5.84 5.84 -6.75
-DEAL::0.00 5.84 5.84 6.75
-DEAL::0.00 9.52 9.52 -11.0
-DEAL::0.00 9.52 9.52 11.0
+DEAL::0 5.84 5.84 -6.75
+DEAL::0 5.84 5.84 6.75
+DEAL::0 9.52 9.52 -11.0
+DEAL::0 9.52 9.52 11.0
DEAL::FE=FESystem<3>[FE_Q<3>(1)-FE_RaviartThomas<3>(1)-FE_Nedelec<3>(0)]
DEAL::component=0
DEAL::2.37 2.37 -133. 2.37 -275. -68.8 -133. -68.8 -138.
DEAL::2.37 2.40e+03 2.44e+03 2.40e+03 4.57e+03 2.40e+03 2.44e+03 2.40e+03 -138.
DEAL::2.37 -2.19e+03 2.44e+03 -2.19e+03 4.79e+03 916. 2.44e+03 916. 5.02e+03
DEAL::component=1
-DEAL::-2.21e-12 -275. -206. -275. -8.99e-13 -138. -206. -138. 1.36e-12
-DEAL::2.77e-12 -275. -206. -275. -459. 4.85e+03 -206. 4.85e+03 -367.
-DEAL::-367. -275. 4.73e+03 -275. 1.41e-12 -138. 4.73e+03 -138. -184.
+DEAL::0 -275. -206. -275. 0 -138. -206. -138. 0
+DEAL::0 -275. -206. -275. -459. 4.85e+03 -206. 4.85e+03 -367.
+DEAL::-367. -275. 4.73e+03 -275. 0 -138. 4.73e+03 -138. -184.
DEAL::367. -275. 4.84e+03 -275. 459. 4.96e+03 4.84e+03 4.96e+03 2.23e+03
-DEAL::-184. 4.62e+03 -206. 4.62e+03 -91.8 -138. -206. -138. -1.36e-12
+DEAL::-184. 4.62e+03 -206. 4.62e+03 -91.8 -138. -206. -138. 0
DEAL::184. 4.73e+03 -206. 4.73e+03 3.52e+03 4.96e+03 -206. 4.96e+03 367.
DEAL::4.80e+03 4.73e+03 4.84e+03 4.73e+03 91.8 -138. 4.84e+03 -138. 184.
DEAL::-4.38e+03 4.83e+03 4.95e+03 4.83e+03 -3.09e+03 5.06e+03 4.95e+03 5.06e+03 -1.81e+03
DEAL::4.97e+03 -68.8 94.6 -68.8 5.60 11.2 94.6 11.2 22.4
DEAL::4.97e+03 2.55e+03 -901. 2.55e+03 5.60 11.2 -901. 11.2 22.4
DEAL::component=4
-DEAL::0.00 5.60 14.0 5.60 -1.84e-16 16.8 14.0 16.8 2.53e-16
-DEAL::0.00 5.60 14.0 5.60 0.00 16.8 14.0 16.8 1.01e-15
-DEAL::0.00 5.60 14.0 5.60 -1.84e-16 16.8 14.0 16.8 1.01e-15
-DEAL::0.00 5.60 14.0 5.60 0.00 16.8 14.0 16.8 0.00
-DEAL::0.00 5.60 14.0 5.60 7.35e-16 16.8 14.0 16.8 2.53e-16
-DEAL::0.00 5.60 14.0 5.60 0.00 16.8 14.0 16.8 1.01e-15
-DEAL::0.00 5.60 14.0 5.60 7.35e-16 16.8 14.0 16.8 1.01e-15
-DEAL::0.00 5.60 14.0 5.60 0.00 16.8 14.0 16.8 0.00
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0
+DEAL::0 5.60 14.0 5.60 0 16.8 14.0 16.8 0