SparsityPattern::Iterator - add a constructor for initalizing an iterator from an...

When trying to increment a sparsity pattern iterator:

  auto iterator = sparsity_pattern.begin();
  iterator++;

The corresponding function in LinearIndexIterator:

  template <class DerivedIterator, class AccessorType>
  inline DerivedIterator
  LinearIndexIterator<DerivedIterator, AccessorType>::operator++(int)
  {
    const DerivedIterator copy(this->accessor);
    operator+=(1);
    return copy;
  }

creates a copy that is initialized by the accessor. Thus, we need a
constructor that can create an interator from an accessor.

diff --git a/include/deal.II/lac/sparsity_pattern.h b/include/deal.II/lac/sparsity_pattern.h
index 89c507a31e1643eb615d3a19d496f12bdddcd43b..573c6c29475bbe9c6fa533bada2c941c0d76f29e 100644 (file)
--- a/include/deal.II/lac/sparsity_pattern.h
+++ b/include/deal.II/lac/sparsity_pattern.h
@@ -298,6 +298,12 @@ namespace SparsityPatternIterators
      * the zeroth row).
      */
     Iterator(const SparsityPattern *sp, const std::size_t linear_index);
+
+    /**
+     * Constructor. Create an iterator into the sparsity pattern @p sp for
+     * a given accessor.
+     */
+    Iterator(const Accessor &accessor);
   };
 } // namespace SparsityPatternIterators
 
@@ -1272,6 +1278,13 @@ namespace SparsityPatternIterators
                             const std::size_t      linear_index)
     : LinearIndexIterator<Iterator, Accessor>(Accessor(sp, linear_index))
   {}
+
+
+  inline Iterator::Iterator(const Accessor &accessor)
+    : LinearIndexIterator<Iterator, Accessor>(accessor)
+  {}
+
+
 } // namespace SparsityPatternIterators