oaLHSUtility.h

Go to the documentation of this file.

 
#ifndef OALHSUTILITY_H
#define OALHSUTILITY_H
 
#include <vector>
 
namespace oalhslib
{
    template <class T>
    void findUniqueColumnElements(const bclib::matrix<T> & A, std::vector<std::vector<T> > & U)
    {
        if (U.size() != A.colsize())
        {
            U = std::vector<std::vector<T> >(A.colsize());
        }
        for (typename bclib::matrix<T>::size_type i = 0; i < A.colsize(); i++)
        {
            U[i] = std::vector<T>();
            for (typename bclib::matrix<T>::const_columnwise_iterator it = A.columnwisebegin(i); it != A.columnwiseend(i); ++it)
            {
                typename std::vector<T>::iterator tempit = std::find(U[i].begin(), U[i].end(), *it);
                if (tempit ==  U[i].end())
                {
                    U[i].push_back(*it);
                }
            }
        }
    }
    
}
 
#endif  /* OALHSUTILITY_H */