* It is assumed that the vector already
* has the right size beforehand.
*/
- void get_dof_values (const dVector &values,
- vector<double> &dof_values) const;
+ void get_dof_values (const dVector &values,
+ dVector &dof_values) const;
+
/**
* Exception
template <>
void
-DoFCellAccessor<1>::get_dof_values (const dVector &values,
- vector<double> &dof_values) const {
+DoFCellAccessor<1>::get_dof_values (const dVector &values,
+ dVector &dof_values) const {
Assert (dof_handler != 0, ExcInvalidObject());
Assert (&dof_handler->get_selected_fe() != 0, ExcInvalidObject());
Assert (dof_values.size() == dof_handler->get_selected_fe().total_dofs,
template <>
void
-DoFCellAccessor<2>::get_dof_values (const dVector &values,
- vector<double> &dof_values) const {
+DoFCellAccessor<2>::get_dof_values (const dVector &values,
+ dVector &dof_values) const {
Assert (dof_handler != 0, ExcInvalidObject());
Assert (&dof_handler->get_selected_fe() != 0, ExcInvalidObject());
Assert (dof_values.size() == dof_handler->get_selected_fe().total_dofs,
// get function values of dofs
// on this cell
- vector<double> dof_values (total_dofs, 0);
+ dVector dof_values (total_dofs);
present_cell->get_dof_values (fe_function, dof_values);
// initialize with zero
// functions
for (unsigned int point=0; point<n_quadrature_points; ++point)
for (unsigned int shape_func=0; shape_func<total_dofs; ++shape_func)
- values[point] += (dof_values[shape_func] *
+ values[point] += (dof_values(shape_func) *
shape_values[selected_dataset](shape_func, point));
};
// get function values of dofs
// on this cell
- vector<double> dof_values (total_dofs, 0);
+ dVector dof_values (total_dofs);
present_cell->get_dof_values (fe_function, dof_values);
// initialize with zero
// functions
for (unsigned int point=0; point<n_quadrature_points; ++point)
for (unsigned int shape_func=0; shape_func<total_dofs; ++shape_func)
- gradients[point] += (dof_values[shape_func] *
+ gradients[point] += (dof_values(shape_func) *
shape_gradients[shape_func][point]);
};