#ifndef  BOOST_SERIALIZATION_UBLAS_VECTOR_HPP
#define BOOST_SERIALIZATION_UBLAS_VECTOR_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_vector.hpp: serialization for boost::numeric::ublas::vector templates

// (C) Copyright 2002 Robert Ramey - http://www.rrsd.com . 
// (C) Copyright 2006 Janek Kozicki - adapted vector.hpp for ublas_vector.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/vector.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>

// 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::vector<T>
template<class Archive, class U, class Allocator>
inline void save(
    Archive & ar,
    const boost::numeric::ublas::vector<U, Allocator> &t,
    const unsigned int /* file_version */
){
    STD::vector<U> tstl; //make temporary copy to STD::vector
    unsigned int count = t.size();
    boost::numeric::ublas::vector<U, Allocator>::const_iterator it=t.begin();
    while(count-- > 0)
        tstl.push_back(*it++);
    boost::serialization::stl::save_collection<Archive, STD::vector<U> >(
        ar, tstl
    );
}

template<class Archive, class U, class Allocator>
inline void load(
    Archive & ar,
    boost::numeric::ublas::vector<U, Allocator> &t,
    const unsigned int /* file_version */
){
    STD::vector<U> tstl;
    boost::serialization::stl::load_collection<
        Archive,
        STD::vector<U>,
        boost::serialization::stl::archive_input_seq<
            Archive, STD::vector<U> 
        >,
        boost::serialization::stl::reserve_imp<STD::vector<U> >
    >(ar, tstl);
    unsigned int count = tstl.size(),n=0;
    STD::vector<U>::const_iterator it=tstl.begin();
    t.resize(count,false);
    while(count-- > 0)
        t(n++)=*it++;
}

// 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::vector<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_VECTOR_HPP