28 #ifndef MAQUIS_ONE_MATRIX_HPP
29 #define MAQUIS_ONE_MATRIX_HPP
48 assert(cols==1 && rows==1);
73 template <
typename T2>
76 template <
typename T2>
79 inline bool empty()
const {
return false; }
84 std::pair<element_iterator, element_iterator>
elements() {
return std::make_pair(val_, NULL); }
85 std::pair<const_element_iterator, const_element_iterator>
elements()
const {
return std::make_pair(val_, NULL); }
101 template <
typename T>
104 template <
typename T>
107 c(0,0) = a(0,0) * b(0,0);
113 template <
typename T>
117 template <
typename T>
121 template <
typename T>
132 template <
typename T>
133 std::ostream& operator << (std::ostream& o, one_matrix<T>
const& m) { o << m(0,0);
return o; }
138 #endif //MAQUIS_ONE_MATRIX_HPP
std::ptrdiff_t difference_type
bool operator==(one_matrix const &rhs) const
size_type num_cols() const
const one_matrix< T > operator*(one_matrix< T > const &m1, one_matrix< T > const &m2)
one_matrix< T > & operator-=(one_matrix const &rhs)
void swap(MPSTensor< Matrix, SymmGroup > &x, MPSTensor< Matrix, SymmGroup > &y)
const_reference const_element_iterator
one_matrix(size_type rows=1, size_type cols=1, T init_value=T())
std::size_t num_rows(maquis::dmrg::one_matrix< T > const &m)
std::pair< const_element_iterator, const_element_iterator > elements() const
std::pair< element_iterator, element_iterator > elements()
value_type const & operator()(const size_type i, const size_type j) const
const one_matrix< T > operator+(one_matrix< T > m1, one_matrix< T > const &m2)
one_matrix< T > & operator/=(T2 const &t)
friend void swap(one_matrix &x, one_matrix &y)
std::size_t num_cols(maquis::dmrg::one_matrix< T > const &m)
std::size_t size_of(one_matrix< T > const &m)
size_type num_rows() const
const one_matrix< T > operator-(one_matrix< T > m1, one_matrix< T > const &m2)
reference element_iterator
T const & const_reference
one_matrix< T > & operator*=(T2 const &t)
void gemm(maquis::dmrg::one_matrix< T > const &a, maquis::dmrg::one_matrix< T > const &b, maquis::dmrg::one_matrix< T > &c)
one_matrix< T > & operator+=(one_matrix const &rhs)
value_type & operator()(const size_type i, const size_type j)