gsXml.hpp

Go to the documentation of this file.

#include <sstream>
#include <gsCore/gsLinearAlgebra.h>
#include <gsIO/gsIOUtils.h>

namespace gismo {

namespace internal {

/*
template<class Object>
std::string gsXml<Object>::tag ()
{
    // Next line will produce compile-time error
    // if gsXml is not specialized for Object
    Object::Object_does_not_exist_ERROR;
    return "";
}
*/

template<class T>
gsXmlNode * makeNode( const std::string & name,
                      const gsMatrix<T> & value, gsXmlTree & data,
                      bool transposed)
{
    return data.allocate_node(rapidxml::node_element ,
                              data.allocate_string(name.c_str() ),
                              makeValue(value,data,transposed) );
}

template<class T>
char * makeValue(const gsMatrix<T> & value, gsXmlTree & data,
                 bool transposed)
{
    std::ostringstream oss;
    // Set precision
    oss << std::setprecision(data.getFloatPrecision()) << "\n";

    // Read/Write is RowMajor
    if ( transposed )
        for ( index_t j = 0; j< value.cols(); ++j)
        {
            for ( index_t i = 0; i< value.rows(); ++i)
                oss << value(i,j) <<" ";
            oss << "\n";
        }
    else
        for ( index_t i = 0; i< value.rows(); ++i)
        {
            for ( index_t j = 0; j< value.cols(); ++j)
                oss << value(i,j) <<" ";
            oss << "\n";
        }

    return data.allocate_string( oss.str().c_str() );
}

template <class T>
void getMatrixFromXml(gsXmlNode* node, unsigned const& rows,
                      unsigned const& cols, gsMatrix<T>& result,
                      const std::string& base_type_flag) {
    // Make sure that flag is in lower case for comparisons
    std::string base_type_flag_;
    base_type_flag_.reserve(base_type_flag.size());
    std::transform(base_type_flag.cbegin(), base_type_flag.cend(),
                   std::back_inserter(base_type_flag_),
                   [](unsigned char c) { return std::tolower(c); });
    if (base_type_flag_ == "ascii") {
        std::istringstream str;
        str.str(node->value());
        result.resize(rows, cols);
        for (unsigned i = 0; i < rows; ++i)  // Read is RowMajor
            for (unsigned j = 0; j < cols; ++j)
                // if ( !(str >> result(i,j) ) )
                if (!gsGetValue(str, result(i, j))) {
                    gsWarn << "XML Warning: Reading matrix of size " << rows
                           << "x" << cols << " failed.\n";
                    gsWarn << "Tag: " << node->name() << ", Matrix entry: ("
                           << i << ", " << j << ").\n";
                    return;
                }
    } else {
        // Read the node-value as the given type and cast into the requested
        // Matrix Scalar Type (T)
        result.resize(rows, cols);
        Base64::DecodeIntoGsType(node->value(), base_type_flag_, result);
    }
}

template<class T>
gsXmlNode * putMatrixToXml ( gsMatrix<T> const & mat, gsXmlTree & data, std::string name)
{
    // Create XML tree node
    gsXmlNode* new_node = internal::makeNode(name, mat, data);
    return new_node;
}

template<class T>
gsXmlNode * putSparseMatrixToXml ( gsSparseMatrix<T> const & mat,
                                   gsXmlTree & data, std::string name)
{
    typedef typename gsSparseMatrix<T>::InnerIterator cIter;

    std::ostringstream str;
    str << std::setprecision(data.getFloatPrecision());
    const index_t nCol = mat.cols();

    for (index_t j=0; j != nCol; ++j) // for all columns
        for ( cIter it(mat,j); it; ++it ) // for all non-zeros in column
        {
            // Write the matrix entry
            str <<it.index() <<" "<<j<<" " << it.value() << "\n";
        }

    // Create XML tree node
    gsXmlNode* new_node = internal::makeNode(name, str.str(), data);
    return new_node;
}

template<class T>
void getSparseEntriesFromXml ( gsXmlNode * node,
                              gsSparseEntries<T> & result )
{
    result.clear();

    std::istringstream str;
    str.str( node->value() );
    index_t r,c;
    T val;

    //while( (str >> r) && (str >> c) && (str >> val) )
    while( (str >> r) && (str >> c) && ( gsGetValue(str,val)) )
        result.add(r,c,val);
}

}// end namespace internal

}// end namespace gismo