SparseMatrix assignment operators

#include <jem/numeric/sparse/operators.h> or
#include <jem/numeric/Sparse.h>.

In addition to the standard assignment operator, jem provides a set of overloaded computed assignment operators for sparse matrices. These are defined as follows:

template <class T>
  const SparseMatrix<T>& operator \<op\>= ( const SparseMatrix<T>& lhs,
                                            const T&               rhs )
{
  return ( lhs = lhs <op> rhs );
}


template <class T>
  const SparseMatrix<T>& operator \<op\>= ( const SparseMatrix<T>& lhs,
                                            const SparseMatrix<T>& rhs )
{
  return ( lhs = lhs <op> rhs );
}

where <op> can be one of the following operators:

• * (multiplication)
• / (division)
• % (modulo division)
• + (addition)
• - (substraction)

The computed assignment operators should be called with two arguments. The first one should be an SparseMatrix object. The second argument should be an SparseMatrix object, or a scalar. All operators return a reference to the first argument.

The computed assignment operators are defined in the header file <jem/numeric/sparse/operators.h>. They are also available through the general header file <jem/numeric/Sparse.h>.

Note that the computed assignment operators are not member functions.