#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.

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

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

Member Function Documentation

◆ operator bool()

libMesh::ConstCouplingAccessor::operator bool ( ) const

inline

Definition at line 144 of file coupling_matrix.h.

References _location, _mat, libMesh::CouplingMatrix::_ranges, std::max(), and libMesh::MeshTools::Subdivision::next.

                         {
     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);
   }