const std::vector<std::shared_ptr<const DataType> > qpd = data_storage->get_data(cell);
- const unsigned int m = matrix_data.m();
const unsigned int n = matrix_data.n();
std::vector<double> single_qp_data(n);
- Assert (qpd.size() == m,
- ExcDimensionMismatch(qpd.size(),m));
+ Assert (qpd.size() == matrix_data.m(),
+ ExcDimensionMismatch(qpd.size(), matrix_data.m()));
for (unsigned int q = 0; q < qpd.size(); q++)
{
qpd[q]->pack_values(single_qp_data);
std::vector<std::shared_ptr<DataType> > qpd = data_storage->get_data(cell);
- const unsigned int m = values_at_qp.m();
const unsigned int n = values_at_qp.n();
std::vector<double> single_qp_data(n);
- Assert(qpd.size() == m,
- ExcDimensionMismatch(qpd.size(),m));
+ Assert(qpd.size() == values_at_qp.m(),
+ ExcDimensionMismatch(qpd.size(), values_at_qp.m()));
for (unsigned int q = 0; q < qpd.size(); q++)
{
const FEValuesBase<dim> &fetest,
double factor = 1.)
{
- unsigned int t_comp = (dim==3) ? dim : 1;
const unsigned int n_dofs = fe.dofs_per_cell;
const unsigned int t_dofs = fetest.dofs_per_cell;
AssertDimension(fe.get_fe().n_components(), dim);
- AssertDimension(fetest.get_fe().n_components(), t_comp);
+ // There should be the right number of components (3 in 3D, otherwise 1)
+ // for the curl.
+ AssertDimension(fetest.get_fe().n_components(), (dim == 3) ? dim : 1);
AssertDimension(M.m(), t_dofs);
AssertDimension(M.n(), n_dofs);