gsPointGrid.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsLinearAlgebra.h>
#include <gsTensor/gsGridIterator.h>

namespace gismo {


/* **************** Utility functions **************** */

template<typename T>  // gsMatrix<T> ab
gsVector<unsigned> uniformSampleCount (const gsVector<T>& lower, 
                                       const gsVector<T>& upper, 
                                       int numPoints = 1000);

template<typename T>
void uniformIntervals(const gsVector<T>& lower, const gsVector<T>& upper, 
                      std::vector< std::vector<T> >& intervals, int numIntervals = 1000);


/* **************** Uniform grids described by limits/corners **************** */


template<class T>  // gsMatrix<T> ab
gsMatrix<T> gsPointGrid( gsVector<T> const & a, gsVector<T> const & b, 
                         gsVector<unsigned> const & np );

// Specialization of the arguments for the 1D case
template<class T> inline
gsMatrix<T> gsPointGrid( T const & t1, T const & t2, unsigned const & n = 100)
{
    gsMatrix<T,1,2> ab(1,2); ab << t1, t2; 
    gsGridIterator<T,CUBE,1> pt(ab,n);
    gsMatrix<T> rvo(ab.rows(), pt.numPoints() );
    for(index_t c = 0; pt; ++pt, ++c)
        rvo.col(c) = *pt;
    return rvo;
}

template<typename T>  // todo: remove, replace by next one
gsMatrix<T>uniformPointGrid(const gsVector<T>& lower, // note: structure lost
                            const gsVector<T>& upper,
                            int numPoints = 1000);

template<class T> inline
gsMatrix<T> gsPointGrid(gsMatrix<T> const & ab, int numPoints)
{   // Note: structure lost
    gsGridIterator<T,CUBE> pt(ab,numPoints);
    gsMatrix<T> rvo(ab.rows(), pt.numPoints() );
    for(index_t c = 0; pt; ++pt, ++c)
        rvo.col(c) = *pt;
    return rvo;
}


/* **************** Cartesian grids **************** */


template<class T, class CwiseContainer> inline
void gsPointGrid(CwiseContainer const & cwise, gsMatrix<T>& res)
{
    gsGridIterator<T,CWISE> pt(cwise);
    res.resize(cwise.size(), pt.numPoints() );
    for(index_t c = 0; pt; ++pt, ++c)
        res.col(c) = *pt;
}

template<class T, class CwiseContainer> inline
gsMatrix<T> gsPointGrid(CwiseContainer const & cwise)
{
    gsMatrix<T> rvo;
    gsPointGrid(cwise, rvo);
    return rvo;
}


} // namespace gismo


#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsPointGrid.hpp)
#endif