* reordered when writing local to
* global).
*/
- std::vector<std::vector<double> > column_values;
+ std::vector<std::vector<value_type> > column_values;
/**
* A mutex variable used to guard access to the member
// comes from an element matrix).
for (size_type j=0; j<n_cols; ++j)
{
- double value = values[j];
+ number value = values[j];
- if (value == 0 && elide_zero_values == true)
+ if (value == number(0) && elide_zero_values == true)
continue;
const std::pair<unsigned int, size_type>
// it comes from an element matrix).
for (size_type j=0; j<n_cols; ++j)
{
- double value = values[j];
+ number value = values[j];
- if (value == 0 && elide_zero_values == true)
+ if (value == number(0) && elide_zero_values == true)
continue;
const std::pair<unsigned int, size_type>
local_matrix, col_ptr, val_ptr);
const size_type n_values = col_ptr - &cols[0];
if (n_values > 0)
- {
- /* @whattodo Note, this is a tricky one! */
- /* global_matrix.add(row, n_values, &cols[0], &vals[0], false, true); */
- Assert ((false), ExcMessage ("This function is corrupt: @whattodo"));
- }
+ global_matrix.add(row, n_values, &cols[0], &vals[0], false, true);
}
}
#endif
#ifdef DEAL_II_WITH_PETSC
-ONLY_MATRIX_FUNCTIONS(PETScWrappers::SparseMatrix);
-ONLY_MATRIX_FUNCTIONS(PETScWrappers::BlockSparseMatrix);
-ONLY_MATRIX_FUNCTIONS(PETScWrappers::MPI::SparseMatrix);
-ONLY_MATRIX_FUNCTIONS(PETScWrappers::MPI::BlockSparseMatrix);
+COMPLEX_ONLY_MATRIX_FUNCTIONS(PETScWrappers::SparseMatrix);
+COMPLEX_ONLY_MATRIX_FUNCTIONS(PETScWrappers::BlockSparseMatrix);
+COMPLEX_ONLY_MATRIX_FUNCTIONS(PETScWrappers::MPI::SparseMatrix);
+COMPLEX_ONLY_MATRIX_FUNCTIONS(PETScWrappers::MPI::BlockSparseMatrix);
#endif
#include "constraint_matrix.inst"