rectangularize.hpp

Go to the documentation of this file.

/*
 * Copyright (C) 1998-2018 ALPS Collaboration. See COPYRIGHT.TXT
 * All rights reserved. Use is subject to license terms. See LICENSE.TXT
 * For use in publications, see ACKNOWLEDGE.TXT
 */

#ifndef ALPS_NUMERIC_RECTANGULARIZE_HPP_08679765945f4171af5c396beff05503
#define ALPS_NUMERIC_RECTANGULARIZE_HPP_08679765945f4171af5c396beff05503

#include <vector>

namespace alps {
    namespace numeric {

        template <typename T>
        void rectangularize(const T&) {}


        template <typename T>
        void rectangularize(std::vector< std::vector<T> >& vec)
        {
            std::size_t mx_size=0;
            for(std::vector<T>& val: vec) {
                // FIXME: // would be needed for 3-tensors, but not quite working:
                // FIXME: rectangularize(val);
                if (mx_size<val.size()) mx_size=val.size();
            }
            for(std::vector<T>& val: vec) {
                val.resize(mx_size);
            }
        }

    }
}

#endif /* ALPS_NUMERIC_RECTANGULARIZE_HPP_08679765945f4171af5c396beff05503 */