#ifndef  BOOST_SERIALIZATION_UBLAS_MATRIX_HPP
#define BOOST_SERIALIZATION_UBLAS_MATRIX_HPP

// MS compatible compilers support #pragma once
#if defined(_MSC_VER) && (_MSC_VER >= 1020)
# pragma once
#endif

/////////1/////////2/////////3/////////4/////////5/////////6/////////7/////////8
// ublas_matrix.hpp: serialization for boost::numeric::ublas::matrix templates

// (C) Copyright 2002 Robert Ramey - http://www.rrsd.com . 
// (C) Copyright 2006 Janek Kozicki - adapted vector.hpp for ublas_matrix.hpp
// Use, modification and distribution is subject to the Boost Software
// License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)

//  See http://www.boost.org for updates, documentation, and revision history.

#include <boost/numeric/ublas/matrix.hpp>
#include <vector>

#include <boost/config.hpp>
#include <boost/detail/workaround.hpp>

#include <boost/serialization/collections_save_imp.hpp>
#include <boost/serialization/collections_load_imp.hpp>
#include <boost/serialization/split_free.hpp>
#include <boost/serialization/vector.hpp>

// function specializations must be defined in the appropriate
// namespace - boost::serialization
#if defined(__SGI_STL_PORT) || defined(_STLPORT_VERSION)
#define STD _STLP_STD
#else
#define STD std
#endif

namespace boost { 
namespace serialization {

/////////1/////////2/////////3/////////4/////////5/////////6/////////7/////////8
// boost::numeric::ublas::matrix<T>
template<class Archive, class U, class Allocator>
inline void save(
    Archive & ar,
    const boost::numeric::ublas::matrix<U, Allocator> &t,
    const unsigned int /* file version */
){
    STD::vector<STD::vector<U> > tstl; //make temporary copy
    unsigned int count1 = t.size1();
    unsigned int count2 = t.size2();
    tstl.resize(count1);
    for(int i=0 ; i<count1 ; ++i)
    {
        tstl[i].resize(count2);
        for(int j=0 ; j<count2 ; ++j)
	    tstl[i][j] = t(i,j);
    }
    ar & tstl;
}

template<class Archive, class U, class Allocator>
inline void load(
    Archive & ar,
    boost::numeric::ublas::matrix<U, Allocator> &t,
    const unsigned int /* file_version */
){
    STD::vector<STD::vector<U> > tstl;
    ar & tstl;
    unsigned int count1 = tstl.size();
    if(count1 == 0)
    {
	    t.resize(0,0);
	    return;
    }
    unsigned int count2 = tstl[0].size();
    t.resize(count1,count2);
    for(int i=0 ; i<count1 ; ++i)
        for(int j=0 ; j<count2 ; ++j)
            t(i,j) = tstl[i][j];
}

// split non-intrusive serialization function member into separate
// non intrusive save/load member functions
template<class Archive, class U, class Allocator>
inline void serialize(
    Archive & ar,
    boost::numeric::ublas::matrix<U, Allocator> & t,
    const unsigned int file_version
){
    boost::serialization::split_free(ar, t, file_version);
}

} // serialization
} // namespace boost

#undef STD

#endif // BOOST_SERIALIZATION_UBLAS_MATRIX_HPP