Nektar::StorageSmvBsr< T >::const_iterator Class Reference

#include <StorageSmvBsr.hpp>

Collaboration diagram for Nektar::StorageSmvBsr< T >::const_iterator:

Classes
struct	IterType

Public Member Functions
	const_iterator (MatrixStorage matType, IndexType begin, IndexType end, IndexType blkDim, const DataVectorType &val, const IndexVectorType &indx, const IndexVectorType &pntr)
	const_iterator (const const_iterator &src)
	~const_iterator ()
const_iterator	operator++ (int)
const_iterator &	operator++ ()
const IterType &	operator* ()
const IterType *	operator-> ()
const bool	operator== (const const_iterator &rhs)
const bool	operator!= (const const_iterator &rhs)

Private Member Functions
void	forward ()
CoordType	storageIndexToFullCoord (IndexType storageIndex)

Private Attributes
MatrixStorage	m_matType
IterType	m_iter
IndexType	m_begin
IndexType	m_end
IndexType	m_blkDim
const DataVectorType &	m_val
const IndexVectorType &	m_indx
const IndexVectorType &	m_pntr

Detailed Description

template<typename T>
class Nektar::StorageSmvBsr< T >::const_iterator

Definition at line 85 of file StorageSmvBsr.hpp.

Constructor & Destructor Documentation

template<typename DataType >

Nektar::StorageSmvBsr< DataType >::const_iterator::const_iterator	(	MatrixStorage	matType,
		IndexType	begin,
		IndexType	end,
		IndexType	blkDim,
		const DataVectorType &	val,
		const IndexVectorType &	indx,
		const IndexVectorType &	pntr
	)

Definition at line 57 of file StorageSmvBsr.cpp.

References Nektar::StorageSmvBsr< T >::begin(), Nektar::StorageSmvBsr< T >::const_iterator::IterType::first, Nektar::NekConstants::kNekSparseNonZeroTol, Nektar::StorageSmvBsr< T >::const_iterator::m_iter, Nektar::StorageSmvBsr< T >::const_iterator::m_val, Nektar::StorageSmvBsr< T >::const_iterator::IterType::nnzindex, Nektar::Array< OneD, const DataType >::num_elements(), Nektar::StorageSmvBsr< T >::const_iterator::IterType::second, Nektar::StorageSmvBsr< T >::const_iterator::IterType::storageindex, and Nektar::StorageSmvBsr< T >::const_iterator::storageIndexToFullCoord().

                                                    :
        m_matType(matType),
        m_iter(),
        m_begin(begin),
        m_end(end),
        m_blkDim(blkDim),
        m_val(val),
        m_indx(indx),
        m_pntr(pntr)
    {
        m_iter.nnzindex = begin;

        // Here we rewind the iterator to the 'begin'-th
        // nonzero double value. Consecutive nonzeros are
        // to lie within dense block and when it's done,
        // jump to the next dense block.
        // NB: m_val stores some explicit zeros that need to be skept.
        if (begin < end)
        {
            // determine offset of 'begin'-th nonzero value

            m_iter.storageindex = 0;
            m_iter.nnzindex     = 0;
            while(m_iter.nnzindex < begin)
            {
                // explicit zero?
                if (m_val[m_iter.storageindex] > NekConstants::kNekSparseNonZeroTol)
                {
                    m_iter.nnzindex++;
                }
                m_iter.storageindex++;
                if (m_iter.storageindex >= m_val.num_elements())
                {
                    std::cout << "const_iterator: 'begin' out stored values bounds" << std::endl;
                    throw 1;
                }
            }
            m_iter.second   = m_val[m_iter.storageindex];

            CoordType c = storageIndexToFullCoord(m_iter.storageindex);
            m_iter.first.first  = c.first;
            m_iter.first.second = c.second;
        }