/*---------------------------------------------------------------------------------*/
/*  NOMAD - Nonlinear Optimization by Mesh Adaptive Direct search -                */
/*                                                                                 */
/*  NOMAD - version 3.9.1 has been created by                                      */
/*                 Charles Audet               - Ecole Polytechnique de Montreal   */
/*                 Sebastien Le Digabel        - Ecole Polytechnique de Montreal   */
/*                 Viviane Rochon Montplaisir - Ecole Polytechnique de Montreal   */
/*                 Christophe Tribes           - Ecole Polytechnique de Montreal   */
/*                                                                                 */
/*  The copyright of NOMAD - version 3.9.1 is owned by                             */
/*                 Sebastien Le Digabel        - Ecole Polytechnique de Montreal   */
/*                 Viviane Rochon Montplaisir - Ecole Polytechnique de Montreal   */
/*                 Christophe Tribes           - Ecole Polytechnique de Montreal   */
/*                                                                                 */
/*  NOMAD v3 has been funded by AFOSR and Exxon Mobil.                             */
/*                                                                                 */
/*  NOMAD v3 is a new version of NOMAD v1 and v2. NOMAD v1 and v2 were created     */
/*  and developed by Mark Abramson, Charles Audet, Gilles Couture, and John E.     */
/*  Dennis Jr., and were funded by AFOSR and Exxon Mobil.                          */
/*                                                                                 */
/*  Contact information:                                                           */
/*    Ecole Polytechnique de Montreal - GERAD                                      */
/*    C.P. 6079, Succ. Centre-ville, Montreal (Quebec) H3C 3A7 Canada              */
/*    e-mail: nomad@gerad.ca                                                       */
/*    phone : 1-514-340-6053 #6928                                                 */
/*    fax   : 1-514-340-5665                                                       */
/*                                                                                 */
/*  This program is free software: you can redistribute it and/or modify it        */
/*  under the terms of the GNU Lesser General Public License as published by       */
/*  the Free Software Foundation, either version 3 of the License, or (at your     */
/*  option) any later version.                                                     */
/*                                                                                 */
/*  This program is distributed in the hope that it will be useful, but WITHOUT    */
/*  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or          */
/*  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License    */
/*  for more details.                                                              */
/*                                                                                 */
/*  You should have received a copy of the GNU Lesser General Public License       */
/*  along with this program. If not, see <http://www.gnu.org/licenses/>.           */
/*                                                                                 */
/*  You can find information on the NOMAD software at www.gerad.ca/nomad           */
/*---------------------------------------------------------------------------------*/

/**
 \file   Signature.cpp
 \brief  Evaluation point signature (implementation)
 \author Sebastien Le Digabel
 \date   2010-04-22
 \see    Signature.hpp
 */
#include "Signature.hpp"
#include "SMesh.hpp"
#include "XMesh.hpp"
#include "GMesh.hpp"

/*-----------------------------------*/
/*   static members initialization   */
/*-----------------------------------*/
#ifdef MEMORY_DEBUG
int NOMAD::Signature::_cardinality     = 0;
int NOMAD::Signature::_max_cardinality = 0;
#endif

bool NOMAD::Signature::_warning_has_been_displayed=false;

/*--------------------------------------------------*/
/*                    constructor 1                 */
/*--------------------------------------------------*/
NOMAD::Signature::Signature ( int                                                n                    ,
                             const std::vector<NOMAD::bb_input_type>           & input_types          ,
                             const NOMAD::Point                                & lb                   ,
                             const NOMAD::Point                                & ub                   ,
                             const NOMAD::mesh_type                            & mesh_type            ,
                             bool                                                anisotropic_mesh     ,
                             NOMAD::Double                                       anisotropy_factor    ,
                             const NOMAD::Point                                & granularity          ,
                             const NOMAD::Point                                & initial_poll_size    ,
                             const NOMAD::Point                                & min_poll_size        ,
                             const NOMAD::Point                                & min_mesh_size        ,
                             NOMAD::Double                                     & mesh_update_basis    ,
                             NOMAD::Double                                     & poll_update_basis    ,
                             int                                               & mesh_coarsening_exponent,
                             int                                               & mesh_refining_exponent,
                             int                                                 initial_mesh_index    ,
                             const NOMAD::Point                                & scaling              ,
                             const NOMAD::Point                                & fixed_variables      ,
                             const std::vector<bool>                           & periodic_variables   ,
                             std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp>   & var_groups           ,
                             const NOMAD::Display                              & out )
:  _std  ( false ) ,
_out ( out )
{
    // Mesh index starts at 0 for (x/g)mesh and decreases as the mesh size decreases
    // The mesh index is reset to 0 as the mesh is reset.
    
    if ( mesh_type == NOMAD::GMESH )
        _mesh= new NOMAD::GMesh(anisotropic_mesh,
                                anisotropy_factor,
                                initial_poll_size,
                                min_poll_size,
                                min_mesh_size,
                                fixed_variables,
                                granularity,
                                poll_update_basis,  // poll update basis and the 2 other parameters are not needed
                                mesh_coarsening_exponent,
                                mesh_refining_exponent);
    
    else if ( mesh_type == NOMAD::XMESH )
        _mesh=new NOMAD::XMesh( anisotropic_mesh,
                               initial_poll_size,
                               min_poll_size,
                               min_mesh_size,
                               fixed_variables,
                               poll_update_basis,  // XMesh sets poll update basis (default 2)
                               mesh_coarsening_exponent,
                               mesh_refining_exponent);
    
    else if ( mesh_type == NOMAD::SMESH )
        // Mesh index can be provided for isotropic smesh and increases as the mesh size decreases
        _mesh=new NOMAD::SMesh(initial_poll_size,   // WARNING THIS IS NOT AN ERROR     new initial_poll_size ~ old initial mesh size
                               min_poll_size,
                               min_mesh_size,
                               fixed_variables,
                               mesh_update_basis, // SMesh set mesh update basis (default 4)
                               mesh_coarsening_exponent,
                               mesh_refining_exponent,
                               initial_mesh_index );
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::Signature(): a valid mesh type must be given" );
    
    init ( n              ,
          input_types        ,
          lb                 ,
          ub                 ,
          scaling            ,
          fixed_variables    ,
          granularity        ,
          periodic_variables ,
          var_groups           );
    
#ifdef MEMORY_DEBUG
    ++NOMAD::Signature::_cardinality;
    if ( NOMAD::Signature::_cardinality > NOMAD::Signature::_max_cardinality )
        ++NOMAD::Signature::_max_cardinality;
#endif
}

/*--------------------------------------------------*/
/*                    constructor 2                 */
/*--------------------------------------------------*/
NOMAD::Signature::Signature ( int                                      n                  ,
                             const std::vector<NOMAD::bb_input_type> & input_types        ,
                             const NOMAD::Point                      & initial_poll_size  ,
                             const NOMAD::Point                      & lb                 ,
                             const NOMAD::Point                      & ub                 ,
                             const std::set<NOMAD::direction_type>   & direction_types    ,
                             const std::set<NOMAD::direction_type>   & sec_poll_dir_types ,
                             const std::set<NOMAD::direction_type>   & int_poll_dir_types ,
                             const NOMAD::Display                    & out                  )
:    _std  ( false ) ,
_out  ( out )
{
    if ( static_cast<int> ( input_types.size() ) != n )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::Signature(): bad argument: input_types" );
    
    // Default mesh is isotropic xmesh
    _mesh=new NOMAD::XMesh (false,
                            initial_poll_size,
                            NOMAD::Point(),
                            NOMAD::Point(),
                            NOMAD::Point(n));
    
    
    // automatic creation of groups of variables:
    // ------------------------------------------
    std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp> var_groups;
    
    {
        std::set<int> vi_cbi;  // list of cont./bin./int. variables
        std::set<int> vi_cat;  // list of categorical variables
        
        for ( int i = 0 ; i < n ; ++i )
            if (  input_types[i] != NOMAD::CATEGORICAL )
                vi_cbi.insert(i);
            else
                vi_cat.insert(i);
        
        // creation of a group for cont./bin./int. variables:
        if ( !vi_cbi.empty() )
            var_groups.insert ( new NOMAD::Variable_Group ( vi_cbi            ,
                                                           direction_types    ,
                                                           sec_poll_dir_types ,
                                                           int_poll_dir_types ,
                                                           out                ) );
        
        // creation of a group for categorical variables:
        if ( !vi_cat.empty() )
            var_groups.insert ( new NOMAD::Variable_Group ( vi_cat            ,
                                                           direction_types    ,
                                                           sec_poll_dir_types ,
                                                           int_poll_dir_types ,
                                                           out                ) );
        
    }
    
    // init:
    // -----
    init ( n                   ,
          input_types         ,
          lb                  ,
          ub                  ,
          NOMAD::Point()      ,
          NOMAD::Point()      ,
          NOMAD::Point()      ,
          std::vector<bool>() ,
          var_groups            );
    
    // delete the temporary groups of variables:
    std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp>::iterator
    it , end = var_groups.end();
    
    for ( it = var_groups.begin() ; it != end ; ++it )
        delete *it;
    
#ifdef MEMORY_DEBUG
    ++NOMAD::Signature::_cardinality;
    if ( NOMAD::Signature::_cardinality > NOMAD::Signature::_max_cardinality )
        ++NOMAD::Signature::_max_cardinality;
#endif
}

/*--------------------------------------------------*/
/*                  copy constructor                */
/*--------------------------------------------------*/
NOMAD::Signature::Signature ( const NOMAD::Signature & s )
:  _lb                 ( s._lb                     ) ,
_ub                 ( s._ub                     ) ,
_scaling            ( s._scaling                ) ,
_fixed_variables    ( s._fixed_variables        ) ,
_granularity        ( s._granularity            ) ,
_input_types        ( s._input_types            ) ,
_all_continuous     ( s._all_continuous         ) ,
_has_categorical    ( s._has_categorical        ) ,
_periodic_variables ( s._periodic_variables        ) ,
_std                ( false                     ) ,
_feas_success_dir   ( s._feas_success_dir       ) ,
_infeas_success_dir ( s._infeas_success_dir     ) ,
_out (s._out)
{
    
    
    if ( dynamic_cast<NOMAD::SMesh*> (s._mesh) )
        _mesh = new NOMAD::SMesh (*(static_cast<NOMAD::SMesh*>(s._mesh)));
    else if ( dynamic_cast<NOMAD::XMesh*> (s._mesh) )
        _mesh = new NOMAD::XMesh (*(static_cast<NOMAD::XMesh*>(s._mesh)));
    else if ( dynamic_cast<NOMAD::GMesh*> (s._mesh) )
        _mesh = new NOMAD::GMesh (*(static_cast<NOMAD::GMesh*>(s._mesh)));
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::Signature(): copy constructor needs a valid mesh_type in the copied object" );
    
    std::list<NOMAD::Variable_Group *>::const_iterator it , end = s._var_groups.end();
    for ( it = s._var_groups.begin() ; it != end ; ++it )
        _var_groups.push_back ( new NOMAD::Variable_Group(**it) );
    
#ifdef MEMORY_DEBUG
    ++NOMAD::Signature::_cardinality;
    if ( NOMAD::Signature::_cardinality > NOMAD::Signature::_max_cardinality )
        ++NOMAD::Signature::_max_cardinality;
#endif
}

/*--------------------------------------------------*/
/*                    destructor                    */
/*--------------------------------------------------*/
NOMAD::Signature::~Signature ( void )
{
    clear();
#ifdef MEMORY_DEBUG
    --NOMAD::Signature::_cardinality;
#endif
}

/*--------------------------------------------------------------------------------*/
/*  clear (private method called by destructor and after an exception is thrown)  */
/*--------------------------------------------------------------------------------*/
void NOMAD::Signature::clear ( void )
{
    _all_continuous  = true;
    _has_categorical = false;
    _std             = false;
    reset_var_groups();
    _feas_success_dir.clear();
    _infeas_success_dir.clear();
    _lb.clear();
    _ub.clear();
    _scaling.clear();
    _fixed_variables.clear();
    _input_types.clear();
    _periodic_variables.clear();
    _granularity.clear();
    delete _mesh;
    
}

/*--------------------------------------------------*/
/*            feasible success direction            */
/*--------------------------------------------------*/
void NOMAD::Signature::set_feas_success_dir ( const NOMAD::Direction & d )
{
    if ( d.size() != static_cast<int>( _input_types.size() ) )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::set_feas_success_dir(): bad direction" );
    _feas_success_dir = d;
}

/*--------------------------------------------------*/
/*            infeasible success direction          */
/*--------------------------------------------------*/
void NOMAD::Signature::set_infeas_success_dir ( const NOMAD::Direction & d )
{
    if ( d.size() != static_cast<int>( _input_types.size() ) )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::set_infeas_success_dir(): bad direction" );
    _infeas_success_dir = d;
}


/*--------------------------------------------------*/
/*               initializations (private)          */
/*--------------------------------------------------*/
void NOMAD::Signature::init ( int                                       n                  ,
                             const std::vector<NOMAD::bb_input_type> & input_types        ,
                             const NOMAD::Point                      & lb                 ,
                             const NOMAD::Point                      & ub                 ,
                             const NOMAD::Point                      & scaling            ,
                             const NOMAD::Point                      & fixed_variables    ,
                             const NOMAD::Point                      & granularity        ,
                             const std::vector<bool>                 & periodic_variables ,
                             std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp> & var_groups )
{
    // reset directions:
    _feas_success_dir.clear();
    _infeas_success_dir.clear();
    
    // check the dimension (n):
    if ( n <= 0 )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: n" );
    
    // bounds:
    if ( lb.empty() )
        _lb.reset ( n );
    else if ( lb.size() == n )
        _lb = lb;
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: lb" );
    if ( ub.empty() )
        _ub.reset ( n );
    else if ( ub.size() == n )
        _ub = ub;
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: ub" );
    
    // scaling:
    if ( scaling.empty() )
        _scaling.reset ( n );
    else if ( scaling.size() == n )
        _scaling = scaling;
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: scaling" );
    int i;
    for ( i = 0 ; i < n ; ++i )
        if ( _scaling[i].is_defined() && _scaling[i] == 0.0 )
            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                     "NOMAD::Signature::init(): bad argument: scaling (zero value)" );
    
    // fixed variables:
    if ( fixed_variables.empty() )
        _fixed_variables.reset ( n );
    else if ( fixed_variables.size() == n )
        _fixed_variables = fixed_variables;
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: fixed_variables" );
    
    if ( granularity.empty() )
        _granularity.reset(n,0) ;
    else if ( granularity.size() == n )
        _granularity = granularity;
    else
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: granular_variables" );
    
    // periodic variables:
    _periodic_variables = periodic_variables;
    if ( !_periodic_variables.empty() )
    {
        
        if ( static_cast<int>(periodic_variables.size()) != n )
            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                     "NOMAD::Signature::init(): bad argument: periodic_variables" );
        
        for ( i = 0 ; i < n ; ++i )
            if ( _periodic_variables[i] && ( !_lb[i].is_defined() || !_ub[i].is_defined() ) )
                throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                         "NOMAD::Signature::init(): incompatible periodic variables" );
    }
    
    // input types:
    if ( static_cast<int>(input_types.size()) != n )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): bad argument: input_types" );
    _input_types     = input_types;
    _all_continuous  = true;
    _has_categorical = false;
    
    for ( i = 0 ; i < n ; ++i )
    {
        
        if ( (_lb[i].is_defined() || _ub[i].is_defined() ) && _input_types[i] == NOMAD::CATEGORICAL && _out.get_gen_dd()>=NOMAD::NORMAL_DISPLAY && !_warning_has_been_displayed  )
        {
            _out << NOMAD::open_block("Warning:")
            << "NOMAD::Signature::init(): Providing explicit bounds for categorical variables is irrelevant. Bounds are provided implicitely by the neighbors function." << std::endl
            << NOMAD::close_block();
            _warning_has_been_displayed=true;
        }
        
        if ( _fixed_variables[i].is_defined() && ((_lb[i].is_defined() && _fixed_variables[i] < _lb[i]) || (_ub[i].is_defined() && _fixed_variables[i] > _ub[i]) ||
                                                  ( (_input_types[i] == NOMAD::INTEGER     ||   _input_types[i] == NOMAD::CATEGORICAL    )
                                                   && !_fixed_variables[i].is_integer() )  ||
                                                  ( _input_types[i] == NOMAD::BINARY && !_fixed_variables[i].is_binary() ) ) )
            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                     "NOMAD::Signature::init(): bad argument: fixed variable inconsistent with bounds or input types" );
        
        if ( _lb[i].is_defined() && _ub[i].is_defined() )
        {
            if ( _lb[i].is_defined() && _ub[i].is_defined() && _lb[i].value()> ub[i].value() )
            {
                throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                         "NOMAD::Signature::init(): bad argument: lower bound must be lower than upper bound!" );
            }
            if ( _lb[i] == _ub[i] )
                _fixed_variables[i]=_lb[i];
        }
    }
    for ( i = 0 ; i < n ; ++i )
    {
        if ( _input_types[i] == NOMAD::CATEGORICAL )
        {
            _has_categorical = true;
            _all_continuous  = false;
            break;
        }
        if ( _input_types[i] != NOMAD::CONTINUOUS )
        {
            _all_continuous = false;
            if ( _has_categorical )
                break;
        }
    }
    
    // variable groups:
    reset_var_groups();
    
    std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp>::iterator
    end = var_groups.end() , it;
    bool mod=false;
    for ( it = var_groups.begin() ; it != end ; ++it )
    {
        
        
        if ( !(*it)->check ( _fixed_variables , input_types , NULL, mod ) )
            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                     "NOMAD::Signature::init(): incompatible variable group" );
        
    }
    if ( mod )
    {
        std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp> tmp_var_groups;
        for ( it = var_groups.begin() ; it != end ; ++it )
        {
            tmp_var_groups.insert(*it);
        }
        var_groups.clear();
        var_groups=tmp_var_groups;
        tmp_var_groups.clear();
        
    }
    
    for ( it = var_groups.begin() ; it != end ; ++it )
        _var_groups.push_back( new NOMAD::Variable_Group (**it) );
    
    // mesh:
    if ( _mesh->get_initial_poll_size().size() != n                             ||
        ( _mesh->get_min_mesh_size().is_defined() && _mesh->get_min_mesh_size().size() != n) ||
        ( _mesh->get_min_poll_size().is_defined() && _mesh->get_min_poll_size().size() != n)    )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::init(): mesh arguments with different sizes" );
    
    
    
}

/*--------------------------------------------------*/
/* Access to the number of categorical variables */
/*--------------------------------------------------*/
int NOMAD::Signature::get_n_categorical ( void ) const
{
    int n_cat=0;
    for (int i = 0 ; i < get_n() ; ++i )
        if ( _input_types[i] == NOMAD::CATEGORICAL )
            n_cat++;
    
    return n_cat;
}

/*--------------------------------------------------*/
/* Access to the number of categorical variables */
/*--------------------------------------------------*/
int NOMAD::Signature::get_nb_fixed_variables ( void ) const
{
    int n_fixed=0;
    for (int i = 0 ; i < get_n() ; ++i )
        if ( _fixed_variables[i].is_defined() )
            n_fixed++;
    
    return n_fixed;
}


/*--------------------------------------------------------*/
/* Access to the maximum number of points in the poll set */
/*--------------------------------------------------------*/
int NOMAD::Signature::get_max_size_poll_set ( void ) const
{
    int m=1; // There is at list one point (the poll center)
    
    NOMAD::Point delta=_mesh->get_delta();
    NOMAD::Point Delta=_mesh->get_Delta();
    
    if ( delta.is_defined() && Delta.is_defined() )
    {
        for (int i = 0 ; i < get_n() ; ++i )
            if ( ! _fixed_variables[i].is_defined() )
            {
                NOMAD::Double r;
                if ( _input_types[i] == NOMAD::CONTINUOUS )
                    r=Delta[i]/delta[i]+1;
                else if ( _input_types[i] == NOMAD::INTEGER )
                    r=Delta[i]+1;
                else if ( _input_types[i] == NOMAD::BINARY )
                    r=2;
                else if ( _input_types[i] == NOMAD::CATEGORICAL )
                    r=1;
                m= m * static_cast<int> ( r.floor().value() );
            }
    }
    
    return m;
}



/*--------------------------------------------------*/
/*                       reset                      */
/*--------------------------------------------------*/
void NOMAD::Signature::reset ( int                                       n                        ,
                              const std::vector<NOMAD::bb_input_type> & input_types              ,
                              const NOMAD::Point                      & lb                       ,
                              const NOMAD::Point                      & ub                       ,
                              const NOMAD::Point                      & scaling                  ,
                              const NOMAD::Point                      & fixed_variables          ,
                              const NOMAD::Point                      & granularity              ,
                              const std::vector<bool>                 & periodic_variables       ,
                              std::set<NOMAD::Variable_Group*,NOMAD::VG_Comp> & var_groups   )
{
    reset_mesh();
    reset_var_groups();
    init ( n                  ,
          input_types        ,
          lb                 ,
          ub                 ,
          scaling            ,
          fixed_variables    ,
          granularity        ,
          periodic_variables ,
          var_groups           );
}

/*--------------------------------------------------*/
/*                reset_var_groups (private)        */
/*--------------------------------------------------*/
void NOMAD::Signature::reset_var_groups ( void )
{
    std::list<NOMAD::Variable_Group *>::const_iterator
    end = _var_groups.end() , it;
    for ( it = _var_groups.begin() ; it != end ; ++it )
        delete *it;
    _var_groups.clear();
}

/*----------------------------------------------------------*/
/*           check the compatibility of a point             */
/*   . we only check the number of variables                */
/*   . other stuff (fixed variables, binary variables,...)  */
/*     will be checked by Eval_Point::check()               */
/*----------------------------------------------------------*/
bool NOMAD::Signature::is_compatible ( const NOMAD::Point & x ) const
{
    if ( get_n() != x.size() || _var_groups.empty() )
        return false;
    return true;
}

/*-----------------------------------------------------*/
/*     compute the directions                          */
/*-----------------------------------------------------*/
void NOMAD::Signature::get_directions ( std::list<NOMAD::Direction>    & dirs     ,
                                       NOMAD::poll_type                poll       ,
                                       const NOMAD::Point            & poll_center   )
{
    
    NOMAD::Direction                          * pd;
    int                                         i;
    std::list<NOMAD::Direction>::const_iterator it_dir , end_dir;
    std::set<int>::const_iterator               it_vi  , end_vi;
    
    int          n = get_n();
    NOMAD::Point delta = _mesh->get_delta ( );
    NOMAD::Point Delta = _mesh->get_Delta ( );
    
    
    // Reset dir_group_index.
    //    For each signature, a variable_group has a  unique set of directions generated and a unique dir_group_index starting by zero (-1 if no dirs)
    _dir_group_index=-1;
    
    // loop on variable groups:
    std::list<NOMAD::Variable_Group*>::const_iterator end_vg = _var_groups.end() , it_vg;
    for ( it_vg = _var_groups.begin() ; it_vg != end_vg ; ++it_vg )
    {
        
        const std::set<int> & var_indexes = (*it_vg)->get_var_indexes();
        
        // get the directions on a unit nc-sphere for the current group of variables:
        std::list<NOMAD::Direction> dirs_nc;
        (*it_vg)->get_directions ( dirs_nc , poll , *_mesh );  // _mesh of all directions
        
        
        // Scaling and projection on the mesh
        end_dir = dirs_nc.end();
        if (static_cast<int>(dirs_nc.size())!=0)
            ++_dir_group_index;
        for ( it_dir = dirs_nc.begin() ; it_dir != end_dir ; ++it_dir )
        {
            
            dirs.push_back ( NOMAD::Direction ( n , 0.0 , it_dir->get_type(),_dir_group_index) );
            
            pd = &(*(--dirs.end()));
            
            end_vi = var_indexes.end();
            
            i = 0;
            NOMAD::Double normInfPoint=0;
            for ( it_vi = var_indexes.begin() ; it_vi != end_vi ; ++it_vi ,++i )
            {
                normInfPoint = NOMAD::max( normInfPoint , (*it_dir)[i].abs() );
            }
            
            
            i = 0;
            for ( it_vi = var_indexes.begin() ; it_vi != end_vi ; ++it_vi ,++i )
            {
                if ( _input_types[*it_vi] == NOMAD::BINARY )
                {
                    if ( (*it_dir)[i] < 0 )
                        (*pd)[*it_vi] = 0.0;
                    else
                        (*pd)[*it_vi] = 1.0;
                }
                // categorical variables: set direction=0:
                else if ( _input_types[*it_vi] == NOMAD::CATEGORICAL )
                    (*pd)[*it_vi] = 0.0;
                else
                {
                    
                    if ( dynamic_cast<NOMAD::GMesh *>(_mesh) )
                    {
                        if ( normInfPoint==0 )
                            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                                     "NOMAD::Signature::set_poll_directions: cannot handle an infinite norm of zero" );
                        
                        (*pd)[*it_vi] = _mesh->scale_and_project(*it_vi,(*it_dir)[i]  / normInfPoint );
                    }
                    else
                    {
                        // Case other than Gmesh:
                        
                        // Do not block Delta and
                        // pickup 0, -1, 1 for the integer variable when Delta < 1.0
                        if ( _input_types[*it_vi] == NOMAD::INTEGER && Delta[i] < 1.0 )
                        {
                            if ( (*it_dir)[i] > 1.0/3.0 )
                                (*pd)[*it_vi] =  1.0 ;
                            else if ( (*it_dir)[i] < -1.0/3.0 )
                                (*pd)[*it_vi] =  -1.0;
                            else
                                (*pd)[*it_vi] =  0.0;
                            
                        }
                        else
                        {
                            // Scaling and projection on the mesh
                            (*pd)[*it_vi] = _mesh->scale_and_project(*it_vi,(*it_dir)[i]);
                            
                            // round for integer variables:
                            if ( _input_types[*it_vi] == NOMAD::INTEGER )
                            {
                                if ( (*pd)[*it_vi] >= Delta[*it_vi]/3.0 )
                                    (*pd)[*it_vi] =  (*pd)[*it_vi].ceil();
                                else if ( (*pd)[*it_vi] <= -Delta[*it_vi]/3.0 )
                                    (*pd)[*it_vi] =  (*pd)[*it_vi].floor();
                                else
                                    (*pd)[*it_vi] =  (*pd)[*it_vi].roundd();
                            }
                        }
                    }
                }
            }
        }
    }
}


/*-----------------------------------------------------------------------------------*/
/*     compute one ORTHO 1 direction for a random group of variable                  */
/*-----------------------------------------------------------------------------------*/
void NOMAD::Signature::get_variable_group_direction ( NOMAD::Direction          & dir        ,
                                                     const NOMAD::Point         & poll_center   )
{
    
    int                                         i=0;
    std::set<int>::const_iterator               it_vi  , end_vi;
    
    NOMAD::Point delta=_mesh->get_delta ( );
    NOMAD::Point Delta=_mesh->get_Delta ( );
    
    
    // Reset dir_group_index.
    //    For each signature, a variable_group has a  unique set of directions generated and a unique dir_group_index starting by zero (-1 if no dirs)
    _dir_group_index=1;
    
    size_t n_vg=_var_groups.size();
    int var_group_id=0;
    
    // Select a variable group;
    if ( n_vg > 1 )
        var_group_id=NOMAD::RNG::rand() % n_vg;
    
    std::list<NOMAD::Variable_Group*>::const_iterator it_vg=_var_groups.begin();
    std::advance(it_vg, var_group_id);
    
    // scale with delta and resize the directions to size n;
    // also round integer and binary variables:
    const std::set<int> & var_indexes = (*it_vg)->get_var_indexes();
    NOMAD::Direction dir_nc ( static_cast<int>(var_indexes.size()) , 0.0 , NOMAD::ORTHO_1 );
    
    // get a direction INSIDE a unit nc-sphere
    if ( (*it_vg)->get_one_direction ( dir_nc ) ) // Get a direction ON a unit nc-sphere
    {
        // Randomly scale the direction to be inside the unit nc-sphere
        NOMAD::Double s = NOMAD::RNG::rand(0,1);
        
        // Make sure to have s!=0
        while ( s==0 )
            s=NOMAD::RNG::rand(0,1);
        
        NOMAD::Double normInfPoint=0;
        for ( int j=0; j < static_cast<int>(var_indexes.size()) ; ++j )
        {
            normInfPoint = NOMAD::max( normInfPoint , (dir_nc)[i].abs() );
        }
        
        end_vi = var_indexes.end();
        i = 0;
        for ( it_vi = var_indexes.begin() ; it_vi != end_vi ; ++it_vi ,++i)
        {
            if ( _input_types[*it_vi] == NOMAD::BINARY )
            {
                if ( dir_nc[i] < 0 )
                    dir_nc[*it_vi] = 0.0;
                else
                    dir_nc[*it_vi] = 1.0;
            }
            // categorical variables: set direction=0:
            else if ( _input_types[*it_vi] == NOMAD::CATEGORICAL )
                dir[*it_vi] = 0.0;
            else
            {
                
                // Scaling and projection on the mesh ---> round up to prevent getting the current incumbent (0<s<=1)
                if ( dynamic_cast<NOMAD::GMesh *>(_mesh) )
                {
                    if ( normInfPoint==0 )
                        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                                 "NOMAD::Signature::set_poll_directions: cannot handle an infinite norm of zero" );
                    
                    dir[*it_vi] = _mesh->scale_and_project(*it_vi,dir_nc[i] * s / normInfPoint );
                }
                else
                {
                    
                    // Case other than Gmesh:
                    
                    // Do not block Delta and
                    // pickup 0, -1, 1 for the integer variable when Delta < 1.0
                    if ( _input_types[*it_vi] == NOMAD::INTEGER && Delta[i] < 1.0 )
                    {
                        if ( dir_nc[i] > 1.0/3.0 )
                            dir[*it_vi] =  1.0 ;
                        else if ( dir_nc[i] < -1.0/3.0 )
                            dir[*it_vi] =  -1.0;
                        else
                            dir[*it_vi] =  0.0;
                        
                    }
                    else
                    {
                        
                        dir[*it_vi] = _mesh->scale_and_project( *it_vi, dir_nc[i] * s, true);
                        
                        // integer variables:
                        if ( _input_types[*it_vi] == NOMAD::INTEGER )
                        {
                            if ( dir[*it_vi] >= Delta[*it_vi]/3.0 )
                                dir[*it_vi] =  dir[*it_vi].ceil();
                            else if ( dir[*it_vi] <= -Delta[*it_vi]/3.0 )
                                dir[*it_vi] =  dir[*it_vi].floor();
                            else
                                dir[*it_vi] =  dir[*it_vi].round();
                        }
                    }
                }
            }
        }
    }
}


/*----------------------------------------------------------------*/
/*  get just one direction for a given mesh (used by VNS search)  */
/*----------------------------------------------------------------*/
void NOMAD::Signature::get_one_direction ( NOMAD::Direction & dir , int ell)
{
    int                           i;
    std::set<int>::const_iterator it_vi  , end_vi;
    
    // get delta_m (mesh size parameter):
    int          n = get_n();
    NOMAD::Point delta=_mesh->get_delta ( );
    NOMAD::Point Delta=_mesh->get_Delta ( );
    
    dir.reset    ( n , 0.0 );
    dir.set_type ( NOMAD::UNDEFINED_DIRECTION );
    
    // The mesh indices are modified and will be reset properly once the one direction is obtained
    const NOMAD::Point old_mesh_indices = _mesh->get_mesh_indices();
    NOMAD::Point modified_mesh_indices( n , NOMAD::Double(ell) );
    _mesh->set_mesh_indices( modified_mesh_indices );
    
    // loop on variable groups:
    std::list<NOMAD::Variable_Group*>::const_iterator end_vg = _var_groups.end() , it_vg;
    for ( it_vg = _var_groups.begin() ; it_vg != end_vg ; ++it_vg )
    {
        
        const std::set<int> & var_indexes = (*it_vg)->get_var_indexes();
        
        NOMAD::Direction dir_nc ( static_cast<int>(var_indexes.size()) , 0.0 , NOMAD::UNDEFINED_DIRECTION );
        
        // get a direction on a unit nc-sphere
        if ( (*it_vg)->get_one_direction ( dir_nc ) )
        {
            
            NOMAD::Double normInfPoint=0;
            for ( int j=0; j < static_cast<int>(var_indexes.size()) ; ++j )
            {
                normInfPoint = NOMAD::max( normInfPoint , (dir_nc)[j].abs() );
            }
            
            
            // scale with delta and round integer and binary variables:
            end_vi = var_indexes.end();
            i      = 0;
            for ( it_vi = var_indexes.begin() ; it_vi != end_vi ; ++it_vi, ++i )
            {
                if ( _input_types[*it_vi] == NOMAD::BINARY )
                {
                    if ( dir_nc[i] < 0 )
                        dir[*it_vi] = 0.0;
                    else
                        dir[*it_vi] = 1.0;
                }
                // categorical variables: set direction=0:
                else if ( _input_types[*it_vi] == NOMAD::CATEGORICAL )
                    dir[*it_vi] = 0.0;
                else
                {
                    if ( dynamic_cast<NOMAD::GMesh *>(_mesh) )
                    {
                        if ( normInfPoint==0 )
                            throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,                                                          "NOMAD::Signature::set_poll_directions: cannot handle an infinite norm of zero" );
                        
                        dir[*it_vi] = _mesh->scale_and_project(*it_vi , dir_nc[i] / normInfPoint );
                    }
                    else
                    {
                        // Case other than Gmesh:
                        
                        // Do not block Delta and
                        // pickup 0, -1, 1 for the integer variable when Delta < 1.0
                        if ( _input_types[*it_vi] == NOMAD::INTEGER && Delta[i] < 1.0 )
                        {
                            if ( dir_nc[i] > 1.0/3.0 )
                                dir[*it_vi] =  1.0 ;
                            else if ( dir_nc[i] < -1.0/3.0 )
                                dir[*it_vi] =  -1.0;
                            else
                                dir[*it_vi] =  0.0;
                            
                        }
                        else
                        {
                            
                            
                            dir[*it_vi]=_mesh->scale_and_project(*it_vi , dir_nc[i]);
                            
                            // integer variables:
                            if ( _input_types[*it_vi] == NOMAD::INTEGER )
                            {
                                if ( dir[*it_vi] >= Delta[*it_vi]/3.0 )
                                    dir[*it_vi] = dir[*it_vi].ceil() ;
                                else if ( dir [*it_vi] <= -Delta[*it_vi]/3.0 )
                                    dir[*it_vi] = dir[*it_vi].floor();
                                else
                                    dir[*it_vi] = dir[*it_vi].round();
                            }
                        }
                    }
                }
            }
        }
    }
    
    // Reset the mesh indices to their previous values
    _mesh->set_mesh_indices( old_mesh_indices );
}

/*----------------------------------*/
/*             scaling              */
/*    (done before an evaluation)   */
/*----------------------------------*/
void NOMAD::Signature::scale ( NOMAD::Point & x )
{
    int n = get_n();
    if ( n != x.size() )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::scale(x): x.size() != signature.size()" );
    NOMAD::Double sci;
    for ( int i = 0 ; i < n ; ++i )
    {
        sci = _scaling[i];
        if ( sci.is_defined() )
            x[i] *= sci;
    }
}

/*----------------------------------*/
/*             unscaling            */
/*    (done after an evaluation)    */
/*----------------------------------*/
void NOMAD::Signature::unscale ( NOMAD::Point & x )
{
    int n = get_n();
    if ( n != x.size() )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::unscale(x): x.size() != signature.size()" );
    NOMAD::Double sci;
    for ( int i = 0 ; i < n ; ++i )
    {
        sci = _scaling[i];
        if ( sci.is_defined() )
            x[i] /= sci;
    }
}

/*-----------------------------------------------------------------------*/
/*                             snap to bounds                            */
/*-----------------------------------------------------------------------*/
/*  . returns true if x has been modified                                */
/*  . supposes that treat_periodic_variables() has already been invoked  */
/*    (if periodic variables have been treated, then bounds are          */
/*     satisfied and there is no need to snap anymore                    */
/*-----------------------------------------------------------------------*/
bool NOMAD::Signature::snap_to_bounds ( NOMAD::Point & x , NOMAD::Direction * direction )
{
    int n = get_n();
    if ( n != x.size() )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::snap_to_bounds(x): x.size() != signature.size()" );
    
    bool modified        = false;
    bool no_periodic_var = _periodic_variables.empty();
    
    for ( int i = 0 ; i < n ; ++i )
        
        if ( no_periodic_var || !_periodic_variables[i] )
        {
            
            const NOMAD::Double & ubi = _ub[i];
            NOMAD::Double       & xi  = x[i];
            if ( ubi.is_defined() && xi > ubi )
            {
                if (direction)
                    (*direction)[i] += ubi - xi;
                xi       = ubi;
                modified = true;
            }
            const NOMAD::Double & lbi = _lb[i];
            if ( lbi.is_defined() && xi < lbi )
            {
                if (direction)
                    (*direction)[i] += (lbi - xi);
                xi       = lbi;
                modified = true;
            }
        }
    
    return modified;
}

/*--------------------------------------------------*/
/*             treat the periodic variables         */
/*         returns true if x has been modified      */
/*--------------------------------------------------*/
bool NOMAD::Signature::treat_periodic_variables ( NOMAD::Point            & x       ,
                                                 const NOMAD::Direction *  old_dir ,
                                                 NOMAD::Direction       *& new_dir   )
{
    if ( _periodic_variables.empty() )
        return false;
    
    int n = get_n();
    if ( n != x.size() )
        throw NOMAD::Signature::Signature_Error ( "Signature.cpp" , __LINE__ , *this ,
                                                 "NOMAD::Signature::treat_periodic_variables(x): x.size() != signature.size()" );
    
    new_dir = ( old_dir ) ? new NOMAD::Direction (*old_dir) : NULL;
    
    bool modified = false;
    
    for ( int i = 0 ; i < n ; ++i )
    {
        
        NOMAD::bb_input_type bbit = _input_types[i];
        
        
        if ( _periodic_variables[i] && !_fixed_variables[i].is_defined() &&
            ( bbit == NOMAD::CONTINUOUS || bbit == NOMAD::INTEGER ) )
        {
            
            const NOMAD::Double & ubi = _ub[i];
            const NOMAD::Double & lbi = _lb[i];
            NOMAD::Double       & xi  = x[i];
            
            bool chk = false;
            
            NOMAD::Double new_x = xi;
            while ( new_x > ubi )
            {
                new_x += lbi - ubi;
                chk    = true;
            }
            
            if ( !chk )
            {
                while ( new_x < lbi )
                {
                    new_x += ubi - lbi;
                    chk    = true;
                }
            }
            
            if ( chk )
            {
                
                if ( bbit == NOMAD::INTEGER )
                    new_x = new_x.round();
                
                if (new_dir)
                    (*new_dir)[i] += new_x - xi;
                
                x[i]     = new_x;
                modified = true;
            }
        }
    }
    
    return modified;
}

/*--------------------------------------------------*/
/*                      display                     */
/*--------------------------------------------------*/
void NOMAD::Signature::display ( const NOMAD::Display & out ) const
{
    if ( _std )
        out << "(standard signature)" << std::endl;
    
    // dimension:
    out << "n                 : "  << get_n() << std::endl;
    
    // bounds:
    out << "lb                : ";
    if ( _lb.is_defined() )
        out << "( " << _lb << ")";
    else
        out << "none";
    out << std::endl;
    out << "ub                : ";
    if ( _ub.is_defined() )
        out << "( " << _ub << ")";
    else
        out << "none";
    out << std::endl;
    
    // scaling:
    out << "scaling           : ";
    if ( _scaling.is_defined() )
        out << "( " << _scaling << ")";
    else
        out << "none";
    out << std::endl;
    
    // fixed variables:
    out << "fixed variables   : ";
    if ( _fixed_variables.is_defined() )
        out << "( " << _fixed_variables << ")";
    else
        out << "none";
    out << std::endl;
    
    // input types:
    out << "input types       : (" << _input_types << " )" << std::endl;
    
    // periodic variables:
    out << "periodic variables: ";
    if ( _periodic_variables.empty() )
        out << "none" << std::endl;
    else
    {
        size_t pvs = _periodic_variables.size();
        out << "{";
        for ( size_t k = 0 ; k < pvs ; ++k )
            out << _periodic_variables[k] << " ";
        out << "}" << std::endl;
    }
    
    // success directions:
    out << "feas. succ. dir.  : ";
    if ( _feas_success_dir.is_defined() )
        out << _feas_success_dir << std::endl;
    else
        out << "none";
    out << std::endl
    << "infeas. succ. dir.: ";
    if ( _infeas_success_dir.is_defined() )
        out << _infeas_success_dir;
    else
        out << "none";
    out << std::endl;
    
    // variable groups:
    out << NOMAD::open_block ( "variable groups" );
    int i = 0;
    std::list<NOMAD::Variable_Group *>::const_iterator end = _var_groups.end() , it;
    for ( it = _var_groups.begin() ; it != end ; ++it )
    {
        out << NOMAD::open_block ( "group #" + NOMAD::itos(i++) )
        << **it << NOMAD::close_block();
    }
    out.close_block();
    
    // mesh:
    out << NOMAD::open_block ( "mesh" ) ;
    out << "initial poll size: ( " << _mesh->get_initial_poll_size() << " )" << std::endl;
    out << "initial mesh size: ( " << _mesh->get_initial_mesh_size() << " )" << std::endl;
    out << "minimum mesh size: ( " << _mesh->get_min_mesh_size() << " )" << std::endl;
    out << "minimum poll size: ( " << _mesh->get_min_poll_size() << " )" << std::endl;
    out << NOMAD::close_block();
}

/*-------------------------------------------*/
/*           comparison operator '<'         */
/*-------------------------------------------*/
/*  (success directions are not considered)  */
/*-------------------------------------------*/
bool NOMAD::Signature::operator < ( const NOMAD::Signature & s ) const
{
    if ( this == &s )
        return false;
    
    // standard signature: not checked: standard and non-standard signatures can be ==
    // ( this is tested in Parameters::check() )
    
    // dimension:
    // ----------
    int  n = static_cast<int>(_lb.size());
    int sn = static_cast<int>(s._lb.size());
    
    if ( n < sn )
        return true;
    if ( sn < n )
        return false;
    
    // variable groups:
    // ----------------
    size_t nvg1 = _var_groups.size();
    size_t nvg2 = s._var_groups.size();
    if ( nvg1 != nvg2 )
        return (nvg1 < nvg2);
    
    std::list<NOMAD::Variable_Group*>::const_iterator
    it1 =   _var_groups.begin() ,
    it2 = s._var_groups.begin() ,
    end =   _var_groups.end();
    
    while ( it1 != end )
    {
        if ( **it1 < **it2 )
            return true;
        if ( **it2 < **it1 )
            return false;
        ++it1;
        ++it2;
    }
    
    // first check on the periodic variables:
    // --------------------------------------
    bool p1_empty = _periodic_variables.empty();
    bool p2_empty = s._periodic_variables.empty();
    if ( p1_empty != p2_empty )
        return p1_empty;
    
    // first check on the mesh:
    // ------------------------
    bool                 chkm          = _mesh->get_min_mesh_size().is_defined();
    bool                 chkp          = _mesh->get_min_poll_size().is_defined();
    
    bool                 s_chkm        = s._mesh->get_min_mesh_size().is_defined();
    bool                 s_chkp        = s._mesh->get_min_poll_size().is_defined();
    
    if ( _mesh->get_initial_mesh_size() != s._mesh->get_initial_mesh_size() &&
        _mesh->get_min_mesh_size() != s._mesh->get_min_mesh_size() &&
        _mesh->get_min_poll_size() != s._mesh->get_min_poll_size() )
    {
        if ( chkm != s_chkm )
            return !chkm;
        if ( chkp != s_chkp )
            return !chkp;
    }
    
    /*---------------------------*/
    /*  loop on all coordinates  */
    /*---------------------------*/
    for ( int i = 0 ; i < n ; ++i )
    {
        
        // input types:
        // ------------
        if ( _input_types[i] < s._input_types[i] )
            return true;
        if ( s._input_types[i] < _input_types[i] )
            return false;
        
        // bounds:
        // -------
        if ( _lb[i].comp_with_undef ( s._lb[i] ) )
            return true;
        if ( s._lb[i].comp_with_undef ( _lb[i] ) )
            return false;
        if ( _ub[i].comp_with_undef ( s._ub[i] ) )
            return true;
        if ( s._ub[i].comp_with_undef ( _ub[i] ) )
            return false;
        
        // scaling:
        // --------
        if ( _scaling[i].comp_with_undef ( s._scaling[i] ) )
            return true;
        if ( s._scaling[i].comp_with_undef ( _scaling[i] ) )
            return false;
        
        // fixed variables:
        // ----------------
        if ( _fixed_variables[i].comp_with_undef ( s._fixed_variables[i] ) )
            return true;
        if ( s._fixed_variables[i].comp_with_undef ( _fixed_variables[i] ) )
            return false;
        
        // periodic variables:
        // -------------------
        if ( !p1_empty && _periodic_variables[i] != s._periodic_variables[i] )
            return _periodic_variables[i];
        
        // granularity:
        // -------------------
        if ( _granularity[i].comp_with_undef ( s._granularity[i] ) )
            return true;
        if ( s._granularity[i].comp_with_undef ( _granularity[i] ) )
            return false;
        
        
        // mesh:
        // -----
        if ( _mesh->get_initial_mesh_size() != s._mesh->get_initial_mesh_size() &&
            _mesh->get_min_mesh_size() != s._mesh->get_min_mesh_size() &&
            _mesh->get_min_poll_size() != s._mesh->get_min_poll_size() )
        {
            if ( _mesh->get_initial_mesh_size()[i].comp_with_undef ( s._mesh->get_initial_mesh_size()[i] ) )
                return true;
            if ( s._mesh->get_initial_mesh_size()[i].comp_with_undef ( _mesh->get_initial_mesh_size()[i] ) )
                return false;
            
            
            if ( chkm )
            {
                if ( _mesh->get_min_mesh_size()[i].comp_with_undef ( s._mesh->get_min_mesh_size()[i] ) )
                    return true;
                if ( s._mesh->get_min_mesh_size()[i].comp_with_undef ( _mesh->get_min_mesh_size()[i] ) )
                    return false;
            }
            if ( chkp )
            {
                if ( _mesh->get_min_poll_size()[i].comp_with_undef ( s._mesh->get_min_poll_size()[i] ) )
                    return true;
                if ( s._mesh->get_min_poll_size()[i].comp_with_undef ( _mesh->get_min_poll_size()[i] ) )
                    return false;
            }
        }
    }
    
    // both signatures are equal:
    return false;
}
