Assert (false, ExcNotImplemented());
};
}
- else if (shape_index < first_quad_index)
+ else if (shape_index < this->first_quad_index)
// ok, dof is on a line
{
const unsigned int line_index
- = (shape_index - first_line_index) / this->dofs_per_line;
+ = (shape_index - this->first_line_index) / this->dofs_per_line;
Assert (line_index < GeometryInfo<dim>::lines_per_cell,
ExcInternalError());
else
Assert (false, ExcNotImplemented());
}
- else if (shape_index < first_hex_index)
+ else if (shape_index < this->first_hex_index)
// dof is on a quad
{
const unsigned int quad_index
- = (shape_index - first_quad_index) / this->dofs_per_quad;
+ = (shape_index - this->first_quad_index) / this->dofs_per_quad;
Assert (static_cast<signed int>(quad_index) <
static_cast<signed int>(GeometryInfo<dim>::quads_per_cell),
ExcInternalError());
// get the global index of the first
// non-diagonal element in this row
const unsigned int rowstart
- = get_sparsity_pattern().get_rowstart_indices()[row] + 1;
+ = this->get_sparsity_pattern().get_rowstart_indices()[row] + 1;
number * const diagonal_element = &this->global_entry(rowstart-1);
number rowsum = 0;
for (unsigned int global_index=rowstart;
global_index<rowstart+rowlength; ++global_index)
- rowsum += std::fabs(global_entry(global_index));
+ rowsum += std::fabs(this->global_entry(global_index));
- *diagonal_element += get_strengthen_diagonal (rowsum, row) * rowsum;
+ *diagonal_element += this->get_strengthen_diagonal (rowsum, row) *
+ rowsum;
}
}
inner_sums.resize (this->m());
// precalc sum(j=k+1, N, a[k][j]))
- for(unsigned int row=0; row<m(); row++) {
+ for(unsigned int row=0; row<this->m(); row++) {
inner_sums[row] = get_rowsum(row);
}