fn_unique.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */

 
#ifndef BELFEM_FN_UNIQUE_HPP
#define BELFEM_FN_UNIQUE_HPP
 
#ifdef BELFEM_BLAZE
#include <memory>
#include <vector>
#include <algorithm>
#endif
 
#include "cl_Vector.hpp"
 
namespace belfem
{
//------------------------------------------------------------------------------

    template < typename T >
    void
    unique( Vector< T > & aVector )
    {
#ifdef BELFEM_ARMADILLO
        // call armadillo interface
        aVector.vector_data() = unique( aVector.vector_data() );
#elif  BELFEM_BLAZE
 
        // get length of vector
        std::size_t tLength = aVector.length();
 
        // get pointer to raw data
        T * tData = aVector.data();
 
        // sort data
        std::sort( tData, tData + tLength );
 
        // make vector unique and resize
        aVector.vector_data().resize(
                std::distance( tData,
                    std::unique( tData, tData + tLength ) ),
                true );
#endif
    }
 
//------------------------------------------------------------------------------
}
 
#endif //BELFEM_FN_UNIQUE_HPP