class InternalData : public FiniteElementBase<dim>::InternalDataBase
{
public:
- /**
- * Destructor. Needed to avoid
- * memory leaks with difference
- * quotients.
- */
- ~InternalData ();
+ /**
+ * Destructor. Needed to avoid
+ * memory leaks with difference
+ * quotients.
+ */
+ ~InternalData ();
/**
* Array with shape function values
* values for each quadrature point.
*/
std::vector<std::vector<double> > shape_values;
+
/**
* Array with shape function gradients
* in quadrature points. There is one
*/
std::vector<std::vector<Tensor<1,dim> > > shape_gradients;
- /**
- * Storage for @p{FEValues}
- * objects needed to
- * approximate second
- * derivatives.
- *
- * The ordering is @p{p+hx},
- * @p{p+hy}, @p{p+hz},
- * @p{p-hx}, @p{p-hy},
- * @p{p-hz}, where unused
- * entries in lower dimensions
- * are missing.
- */
- std::vector<FEValues<dim>*> differences;
+ /**
+ * Storage for @p{FEValues}
+ * objects needed to
+ * approximate second
+ * derivatives.
+ *
+ * The ordering is @p{p+hx},
+ * @p{p+hy}, @p{p+hz},
+ * @p{p-hx}, @p{p-hy},
+ * @p{p-hz}, where unused
+ * entries in lower dimensions
+ * are missing.
+ */
+//TODO: what good is this, the base class already has it and it seems never to be used!?
+ std::vector<FEValues<dim>*> differences;
};
/**
* Initialize all vectors to correct size.
*/
void initialize (const unsigned int n_quadrature_points,
- unsigned int n_shapes,
- const UpdateFlags flags);
+ const unsigned int n_shapes,
+ const UpdateFlags flags);
/**
* Storage type for shape values.
* constructor of the derived
* class.
*/
+//TODO: What is the meaning of the n_values_array? It is not used in the implementation
FEValuesBase (const unsigned int n_q_points,
const unsigned int dofs_per_cell,
const unsigned int n_values_array,