1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra. 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2012 Desire NUENTSA WAKAM <desire.nuentsa_wakam@inria.fr> 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef EIGEN_BROWSE_MATRICES_H 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define EIGEN_BROWSE_MATRICES_H 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathnamespace Eigen { 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathenum { 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath SPD = 0x100, 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath NonSymmetric = 0x0 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath/** 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * @brief Iterator to browse matrices from a specified folder 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * This is used to load all the matrices from a folder. 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * The matrices should be in Matrix Market format 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * It is assumed that the matrices are named as matname.mtx 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * and matname_SPD.mtx if the matrix is Symmetric and positive definite (or Hermitian) 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * The right hand side vectors are loaded as well, if they exist. 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * They should be named as matname_b.mtx. 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Note that the right hand side for a SPD matrix is named as matname_SPD_b.mtx 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Sometimes a reference solution is available. In this case, it should be named as matname_x.mtx 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Sample code 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \code 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \endcode 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * \tparam Scalar The scalar type 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate <typename Scalar> 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathclass MatrixMarketIterator 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath public: 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Matrix<Scalar,Dynamic,1> VectorType; 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef SparseMatrix<Scalar,ColMajor> MatrixType; 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath public: 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixMarketIterator(const std::string folder):m_sym(0),m_isvalid(false),m_matIsLoaded(false),m_hasRhs(false),m_hasrefX(false),m_folder(folder) 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_folder_id = opendir(folder.c_str()); 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (!m_folder_id){ 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_isvalid = false; 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cerr << "The provided Matrix folder could not be opened \n\n"; 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath abort(); 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Getnextvalidmatrix(); 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ~MatrixMarketIterator() 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_folder_id) closedir(m_folder_id); 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline MatrixMarketIterator& operator++() 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_matIsLoaded = false; 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasrefX = false; 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasRhs = false; 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Getnextvalidmatrix(); 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return *this; 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline operator bool() const { return m_isvalid;} 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** Return the sparse matrix corresponding to the current file */ 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline MatrixType& matrix() 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Read the matrix 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_matIsLoaded) return m_mat; 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string matrix_file = m_folder + "/" + m_matname + ".mtx"; 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if ( !loadMarket(m_mat, matrix_file)) 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_matIsLoaded = false; 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return m_mat; 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_matIsLoaded = true; 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_sym != NonSymmetric) 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { // Store the upper part of the matrix. It is needed by the solvers dealing with nonsymmetric matrices ?? 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType B; 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath B = m_mat; 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_mat = B.template selfadjointView<Lower>(); 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return m_mat; 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** Return the right hand side corresponding to the current matrix. 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * If the rhs file is not provided, a random rhs is generated 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline VectorType& rhs() 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Get the right hand side 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_hasRhs) return m_rhs; 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string rhs_file; 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath rhs_file = m_folder + "/" + m_matname + "_b.mtx"; // The pattern is matname_b.mtx 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasRhs = Fileexists(rhs_file); 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_hasRhs) 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_rhs.resize(m_mat.cols()); 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasRhs = loadMarketVector(m_rhs, rhs_file); 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (!m_hasRhs) 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Generate a random right hand side 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (!m_matIsLoaded) this->matrix(); 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_refX.resize(m_mat.cols()); 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_refX.setRandom(); 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_rhs = m_mat * m_refX; 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasrefX = true; 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasRhs = true; 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return m_rhs; 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath /** Return a reference solution 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * If it is not provided and if the right hand side is not available 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * then refX is randomly generated such that A*refX = b 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * where A and b are the matrix and the rhs. 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath * Note that when a rhs is provided, refX is not available 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath */ 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline VectorType& refX() 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Check if a reference solution is provided 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_hasrefX) return m_refX; 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string lhs_file; 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath lhs_file = m_folder + "/" + m_matname + "_x.mtx"; 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasrefX = Fileexists(lhs_file); 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_hasrefX) 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_refX.resize(m_mat.cols()); 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_hasrefX = loadMarketVector(m_refX, lhs_file); 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return m_refX; 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline std::string& matname() { return m_matname; } 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline int sym() { return m_sym; } 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline bool hasRhs() {return m_hasRhs; } 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline bool hasrefX() {return m_hasrefX; } 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath protected: 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath inline bool Fileexists(std::string file) 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::ifstream file_id(file.c_str()); 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (!file_id.good() ) 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return false; 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath else 166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath file_id.close(); 168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return true; 169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath void Getnextvalidmatrix( ) 173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_isvalid = false; 175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Here, we return with the next valid matrix in the folder 176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath while ( (m_curs_id = readdir(m_folder_id)) != NULL) { 177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_isvalid = false; 178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string curfile; 179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath curfile = m_folder + "/" + m_curs_id->d_name; 180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Discard if it is a folder 181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if (m_curs_id->d_type == DT_DIR) continue; //FIXME This may not be available on non BSD systems 182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// struct stat st_buf; 183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// stat (curfile.c_str(), &st_buf); 184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// if (S_ISDIR(st_buf.st_mode)) continue; 185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Determine from the header if it is a matrix or a right hand side 187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool isvector,iscomplex; 188c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(!getMarketHeader(curfile,m_sym,iscomplex,isvector)) continue; 189c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(isvector) continue; 190c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 191c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Get the matrix name 192c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string filename = m_curs_id->d_name; 193c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_matname = filename.substr(0, filename.length()-4); 194c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 195c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // Find if the matrix is SPD 196c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath size_t found = m_matname.find("SPD"); 197c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if( (found!=std::string::npos) && (m_sym != NonSymmetric) ) 198c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_sym = SPD; 199c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 200c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath m_isvalid = true; 201c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath break; 202c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 203c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 204c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int m_sym; // Symmetry of the matrix 205c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixType m_mat; // Current matrix 206c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VectorType m_rhs; // Current vector 207c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath VectorType m_refX; // The reference solution, if exists 208c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string m_matname; // Matrix Name 209c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool m_isvalid; 210c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool m_matIsLoaded; // Determine if the matrix has already been loaded from the file 211c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool m_hasRhs; // The right hand side exists 212c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath bool m_hasrefX; // A reference solution is provided 213c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::string m_folder; 214c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath DIR * m_folder_id; 215c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath struct dirent *m_curs_id; 216c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 217c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 218c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 219c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} // end namespace Eigen 220c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 221c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 222