* matrix is of dimension
* $m \times n$.
*/
- unsigned int m () const;
+ types::global_dof_index m () const;
/**
* Return the dimension of the
* matrix is of dimension
* $m \times n$.
*/
- unsigned int n () const;
+ types::global_dof_index n () const;
/**
* Return the number of nonzero
* is allowed to store zero
* values in non-existent fields.
*/
- void set (const unsigned int i,
- const unsigned int j,
+ void set (const types::global_dof_index i,
+ const types::global_dof_index j,
const number value);
/**
* is allowed to store zero
* values in non-existent fields.
*/
- void add (const unsigned int i,
- const unsigned int j,
+ void add (const types::global_dof_index i,
+ const types::global_dof_index j,
const number value);
/**
* tailored better to a sparse matrix
* structure.
*/
- number operator () (const unsigned int i,
- const unsigned int j) const;
+ number operator () (const types::global_dof_index i,
+ const types::global_dof_index j) const;
/**
* This function is mostly like
* tailored better to a sparse matrix
* structure.
*/
- number el (const unsigned int i,
- const unsigned int j) const;
+ number el (const types::global_dof_index i,
+ const types::global_dof_index j) const;
/**
* Return the main diagonal
* involve searching for the
* right column number.
*/
- number diag_element (const unsigned int i) const;
+ number diag_element (const types::global_dof_index i) const;
/**
* Same as above, but return a
* writeable reference. You're
* sure you know what you do?
*/
- number & diag_element (const unsigned int i);
+ number & diag_element (const types::global_dof_index i);
//@}
/**
* Return the location of entry
* $(i,j)$ within the val array.
*/
- unsigned int compute_location (const unsigned int i,
- const unsigned int j) const;
+ unsigned int compute_location (const types::global_dof_index i,
+ const types::global_dof_index j) const;
// make all other sparse matrices
// friends
template <typename number>
inline
-unsigned int ChunkSparseMatrix<number>::m () const
+types::global_dof_index ChunkSparseMatrix<number>::m () const
{
Assert (cols != 0, ExcNotInitialized());
return cols->rows;
template <typename number>
inline
-unsigned int ChunkSparseMatrix<number>::n () const
+types::global_dof_index ChunkSparseMatrix<number>::n () const
{
Assert (cols != 0, ExcNotInitialized());
return cols->cols;
template <typename number>
inline
unsigned int
-ChunkSparseMatrix<number>::compute_location (const unsigned int i,
- const unsigned int j) const
+ChunkSparseMatrix<number>::compute_location (const types::global_dof_index i,
+ const types::global_dof_index j) const
{
const unsigned int chunk_size = cols->get_chunk_size();
const unsigned int chunk_index
template <typename number>
inline
-void ChunkSparseMatrix<number>::set (const unsigned int i,
- const unsigned int j,
+void ChunkSparseMatrix<number>::set (const types::global_dof_index i,
+ const types::global_dof_index j,
const number value)
{
template <typename number>
inline
-void ChunkSparseMatrix<number>::add (const unsigned int i,
- const unsigned int j,
+void ChunkSparseMatrix<number>::add (const types::global_dof_index i,
+ const types::global_dof_index j,
const number value)
{
template <typename number>
inline
-number ChunkSparseMatrix<number>::operator () (const unsigned int i,
- const unsigned int j) const
+number ChunkSparseMatrix<number>::operator () (const types::global_dof_index i,
+ const types::global_dof_index j) const
{
Assert (cols != 0, ExcNotInitialized());
AssertThrow (compute_location(i,j) != SparsityPattern::invalid_entry,
template <typename number>
inline
-number ChunkSparseMatrix<number>::el (const unsigned int i,
- const unsigned int j) const
+number ChunkSparseMatrix<number>::el (const types::global_dof_index i,
+ const types::global_dof_index j) const
{
Assert (cols != 0, ExcNotInitialized());
const unsigned int index = compute_location(i,j);
template <typename number>
inline
-number ChunkSparseMatrix<number>::diag_element (const unsigned int i) const
+number ChunkSparseMatrix<number>::diag_element (const types::global_dof_index i) const
{
Assert (cols != 0, ExcNotInitialized());
Assert (cols->optimize_diagonal(), ExcNotQuadratic());
template <typename number>
inline
-number & ChunkSparseMatrix<number>::diag_element (const unsigned int i)
+number & ChunkSparseMatrix<number>::diag_element (const types::global_dof_index i)
{
Assert (cols != 0, ExcNotInitialized());
Assert (cols->optimize_diagonal(), ExcNotQuadratic());