* mapping of cells at the
* boundary.
*/
- TensorProductPolynomials<dim>* tensor_pols;
+ TensorProductPolynomials<dim> *tensor_pols;
/**
* Number of the Qp tensor
/**
* Needed for inner faces.
*/
+//TODO: can we make this variable static?
StraightBoundary<dim> straight_boundary;
/**
* Used to test the covariant
* transformation.
*/
+//TODO: why have two ways to compute? if they both work, choose one and remove the other
bool alternative_normals_computation;
/**
*
* The default value is false.
*/
+//TODO: remove use_mapping_q_on_all_cells as it is set to false in the constructor and never set again
bool use_mapping_q_on_all_cells;
};
tensor_pols(0),
n_shape_functions(2),
renumber(0),
+//TODO: why have two ways to compute? if they both work, choose one and remove the other
alternative_normals_computation(false),
+//TODO: remove use_mapping_q_on_all_cells as it is set to false in the constructor and never set again
use_mapping_q_on_all_cells(false)
{}
-template<typename number> static
-number power(number x, unsigned int y)
+template<typename number>
+static number power(const number x, const unsigned int y)
{
number value=1;
for (unsigned int i=0; i<y; ++i)
template<int dim>
-MappingQ<dim>::MappingQ (unsigned int p):
+MappingQ<dim>::MappingQ (const unsigned int p):
laplace_on_quad_vector(0),
laplace_on_hex_vector(0),
degree(p),
tensor_pols(0),
n_shape_functions(0),
renumber(0),
+//TODO: why have two ways to compute? if they both work, choose one and remove the other
alternative_normals_computation(false),
+//TODO: remove use_mapping_q_on_all_cells as it is set to false in the constructor and never set again
use_mapping_q_on_all_cells(false)
{
// Construct the tensor product