FermatSelector.h

Go to the documentation of this file.

/*
 * Copyright (C) <2012> <EDF-R&D> <FRANCE>
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 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 General Public License for more details.
 * You should have received a copy of the GNU General Public License along
 * with this program; if not, write to the Free Software Foundation, Inc.,
 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 */
 
#ifndef FERMAT_SELECTOR
#define FERMAT_SELECTOR
 
#include "Selector.h"
 
template <typename T> class FermatSelector : public Selector<T>
{
public:
    FermatSelector() : Selector<T>() {}
    virtual Selector<T>* Copy()
    {
        FermatSelector* newSelector = new FermatSelector();
        newSelector->setIsDeletable(this->deletable);
        return newSelector;
    }
 
    virtual SELECTOR_RESPOND canBeInserted(T* r, unsigned long long& replace)
    {
        vec3 receptorPos(static_cast<Recepteur*>(r->getRecepteur())->getPosition());
        vec3 finalPos(r->getFinalPosition());
 
        vec3 closestPoint;
        // find closest position (closestPoint) from the receptor on the segment receptorPos-finalPos and
        // compute the trueLength (ray length from source to closestPoint)
        decimal trueLength = r->computeTrueLength(receptorPos, finalPos, closestPoint);
 
        // compute the thickness of the ray for the previsouly computed trueLength
        decimal epaisseur = r->getThickness(trueLength, false);
 
        // compute the distance from  the receptor and its closest point on the ray
        decimal closestDistance =
            static_cast<Recepteur*>(r->getRecepteur())->getPosition().distance(closestPoint);
 
        // ray should be rejected if the receptor is outside the cone of the ray (i.e thickness/2 <=
        // closestsDistance)
        if (closestDistance >= (epaisseur / 2.))
        {
            return SELECTOR_REJECT;
        }
 
        return SELECTOR_ACCEPT;
    }
    virtual void insert(T* r, unsigned long long& replace)
    {
        return;
    }
 
    virtual bool insertWithTest(T* r)
    {
        vec3 receptorPos(static_cast<Recepteur*>(r->getRecepteur())->getPosition());
        vec3 finalPos(r->getFinalPosition());
 
        vec3 closestPoint;
        // find closest position (closestPoint) from the receptor on the segment receptorPos-finalPos ray and
        // compute the trueLength (ray length from source to closestPoint)
        decimal trueLength = r->computeTrueLength(receptorPos, finalPos, closestPoint);
 
        // compute the thickness of the ray for the previsouly computed trueLength
        decimal epaisseur = r->getThickness(trueLength, false);
 
        // compute the distance from  the receptor and its closest point on the ray
        decimal closestDistance =
            static_cast<Recepteur*>(r->getRecepteur())->getPosition().distance(closestPoint);
 
        // ray should be rejected if the receptor is outside the cone of the ray (i.e thickness/2 <=
        // closestsDistance)
        if (closestDistance >= (epaisseur / 2.))
        {
            return false;
        }
 
        return true;
    }
 
    virtual const char* getSelectorName()
    {
        return typeid(this).name();
    }
 
protected:
}; // FERMAT_SELECTOR
 
#endif