}
// Apply inverse diagonal
this->inverse_vmult(block, x_cell, b_cell);
+ #ifdef DEBUG
+ for (unsigned int i=0;i<x_cell.size();++i)
+ {
+ Assert(numbers::is_finite(x_cell(i)), ExcNumberNotFinite());
+ }
+ #endif
// Store in result vector
row=additional_data->block_list.begin(block);
- for (unsigned int row_cell=0; row_cell<bs; ++row_cell, ++row)
+ for (size_type row_cell=0; row_cell<bs; ++row_cell, ++row)
dst(row->column()) = prev(row->column()) + additional_data->relaxation * x_cell(row_cell);
}
}
*/
void assemble(MATRIX &G,
const FullMatrix<double> &M,
- const std::vector<unsigned int> &i1,
- const std::vector<unsigned int> &i2);
+ const std::vector<types::global_dof_index> &i1,
+ const std::vector<types::global_dof_index> &i2);
/**
- * Assemble a single matrix
- * into a global matrix.
+ * Assemble a single matrix into a global matrix.
*/
void assemble(MATRIX &G,
const FullMatrix<double> &M,
void assemble_up(MATRIX &G,
const FullMatrix<double> &M,
- const std::vector<unsigned int> &i1,
- const std::vector<unsigned int> &i2,
+ const std::vector<types::global_dof_index> &i1,
+ const std::vector<types::global_dof_index> &i2,
const unsigned int level = numbers::invalid_unsigned_int);
/**
- * Assemble a single matrix
- * into a global matrix.
+ * Assemble a single matrix into a global matrix.
*/
void assemble_down(MATRIX &G,