This accessor class allows simple access to CouplingMatrix values. More...

#include <coupling_matrix.h>

Inheritance diagram for libMesh::ConstCouplingAccessor:

Public Member Functions
	ConstCouplingAccessor (std::size_t loc_in, const CouplingMatrix &mat_in)

	operator bool () const

Protected Attributes
std::size_t	_location

const CouplingMatrix &	_mat

Detailed Description

This accessor class allows simple access to CouplingMatrix values.

Definition at line 134 of file coupling_matrix.h.

Constructor & Destructor Documentation

◆ ConstCouplingAccessor()

libMesh::ConstCouplingAccessor::ConstCouplingAccessor	(	std::size_t	loc_in,
		const CouplingMatrix &	mat_in
	)

inline

Definition at line 137 of file coupling_matrix.h.

                                                       :
    _location(loc_in), _mat(mat_in)
  {
    libmesh_assert_less(_location, _mat.size() * _mat.size());
  }

References _location, _mat, and libMesh::CouplingMatrix::size().

Member Function Documentation

◆ operator bool()

libMesh::ConstCouplingAccessor::operator bool ( ) const

inline

Definition at line 144 of file coupling_matrix.h.

  {
    const std::size_t max_size = std::numeric_limits<std::size_t>::max();
 
    // Find the range that might contain i,j
    // lower_bound isn't *quite* what we want
    CouplingMatrix::rc_type::const_iterator lb = std::upper_bound
      (_mat._ranges.begin(), _mat._ranges.end(),
       std::make_pair(_location, max_size));
    if (lb!=_mat._ranges.begin())
      --lb;
    else
      lb=_mat._ranges.end();
 
    // If no range might contain i,j then it's 0
    if (lb == _mat._ranges.end())
      return false;
 
    const std::size_t lastloc  = lb->second;
 
#ifdef DEBUG
    const std::size_t firstloc = lb->first;
    libmesh_assert_less_equal(firstloc, lastloc);
    libmesh_assert_less_equal(firstloc, _location);
 
    CouplingMatrix::rc_type::const_iterator next = lb;
    next++;
    if (next != _mat._ranges.end())
      {
        // Ranges should be sorted and should not touch
        libmesh_assert_greater(next->first, lastloc+1);
      }
#endif
 
    return (lastloc >= _location);
  }