block_diag_matrix_to_string.h

//-*- mode: c++ -*-
#ifndef _BLOCK_DIAG_MATRIX_TO_STRING_H_
#define _BLOCK_DIAG_MATRIX_TO_STRING_H_

#include <string>
#include <sstream>

template <class el_type>
std::string block_diag_matrix<el_type> :: to_string() const
{
  std::ostringstream os;
#ifdef SYM_ILDL_DEBUG
  os << "Block Diagonal Matrix (" << n_rows() << ", " << n_cols() << ", " << nnz() << ")" << std::endl;

  os << "Main Diagonal Values     = " << main_diag << std::endl;
  os << "Off Diagonal (col, val)  = " << "[";
  for (int i = 0; i < n_cols(); i++) {
    if (block_size(i) == 2) {
        os << "(" << i << ", " << off_diag.find(i)->second << "), ";
        i++;
    }
  }
  os << "]";
#endif
  
  return os.str();
}

#endif