for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<2,dim> hessian = fe_face_values[single_component].hessian (i,q);
- Tensor<3,dim> hessian_normal_outer_prod;
-
- outer_product(hessian_normal_outer_prod, hessian, fe_face_values.normal_vector(q));
+ Tensor<3, dim> hessian_normal_outer_prod = outer_product(
+ hessian, fe_face_values.normal_vector(q));
boundary_integral += hessian_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<2,dim> hessian = fe_face_values[single_component].hessian (i,q);
- Tensor<3,dim> hessian_normal_outer_prod;
-
- outer_product(hessian_normal_outer_prod, hessian, fe_face_values.normal_vector(q));
+ Tensor<3, dim> hessian_normal_outer_prod = outer_product(
+ hessian, fe_face_values.normal_vector(q));
boundary_integral += hessian_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<1,dim> gradient = fe_face_values[single_component].gradient (i,q);
- Tensor<2,dim> gradient_normal_outer_prod;
-
- outer_product(gradient_normal_outer_prod, gradient, fe_face_values.normal_vector(q));
+ Tensor<2, dim> gradient_normal_outer_prod = outer_product(
+ gradient, fe_face_values.normal_vector(q));
boundary_integral += gradient_normal_outer_prod * fe_face_values.JxW(q);
}
}
// ---------------------------------------------------------------------
//
-// Copyright (C) 2003 - 2014 by the deal.II authors
+// Copyright (C) 2003 - 2015 by the deal.II authors
//
// This file is part of the deal.II library.
//
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<1,dim> gradient = fe_face_values[single_component].gradient (i,q);
- Tensor<2,dim> gradient_normal_outer_prod;
-
- outer_product(gradient_normal_outer_prod, gradient, fe_face_values.normal_vector(q));
+ Tensor<2, dim> gradient_normal_outer_prod = outer_product(
+ gradient, fe_face_values.normal_vector(q));
boundary_integral += gradient_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<2,dim> hessian = fe_face_values[single_component].hessian (i,q);
- Tensor<3,dim> hessian_normal_outer_prod;
-
- outer_product(hessian_normal_outer_prod, hessian, fe_face_values.normal_vector(q));
+ Tensor<3, dim> hessian_normal_outer_prod = outer_product(
+ hessian, fe_face_values.normal_vector(q));
boundary_integral += hessian_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<1,dim> gradient = fe_face_values[single_component].gradient (i,q);
- Tensor<2,dim> gradient_normal_outer_prod;
-
- outer_product(gradient_normal_outer_prod, gradient, fe_face_values.normal_vector(q));
+ Tensor<2, dim> gradient_normal_outer_prod = outer_product(
+ gradient, fe_face_values.normal_vector(q));
boundary_integral += gradient_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<2,dim> hessian = fe_face_values[single_component].hessian (i,q);
- Tensor<3,dim> hessian_normal_outer_prod;
-
- outer_product(hessian_normal_outer_prod, hessian, fe_face_values.normal_vector(q));
+ Tensor<3, dim> hessian_normal_outer_prod = outer_product(
+ hessian, fe_face_values.normal_vector(q));
boundary_integral += hessian_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<1,dim> gradient = fe_face_values[single_component].gradient (i,q);
- Tensor<2,dim> gradient_normal_outer_prod;
-
- outer_product(gradient_normal_outer_prod, gradient, fe_face_values.normal_vector(q));
+ Tensor<2, dim> gradient_normal_outer_prod = outer_product(
+ gradient, fe_face_values.normal_vector(q));
boundary_integral += gradient_normal_outer_prod * fe_face_values.JxW(q);
}
}
for (unsigned int q=0; q<fe_face_values.n_quadrature_points; ++q)
{
Tensor<2,dim> hessian = fe_face_values[single_component].hessian (i,q);
- Tensor<3,dim> hessian_normal_outer_prod;
-
- outer_product(hessian_normal_outer_prod, hessian, fe_face_values.normal_vector(q));
+ Tensor<3, dim> hessian_normal_outer_prod = outer_product(
+ hessian, fe_face_values.normal_vector(q));
boundary_integral += hessian_normal_outer_prod * fe_face_values.JxW(q);
}
}