fastmarchingmethod.h

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2013   Borek Patzak
 *
 *
 *
 *       Czech Technical University, Faculty of Civil Engineering,
 *   Department of Structural Mechanics, 166 29 Prague, Czech Republic
 *
 *  This library 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 2.1 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 library; if not, write to the Free Software
 *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */

#ifndef fastmarchingmethod_h
#define fastmarchingmethod_h

#include "floatarray.h"
#include "mathfem.h"

#include <vector>
#include <list>
#include <queue>

namespace oofem {
class Domain;

class FastMarchingMethod
{
protected:

    enum FNM_Status_Type {
        FMM_Status_FAR,            
        FMM_Status_TRIAL,          
        FMM_Status_KNOWN,          
        FMM_Status_KNOWN_BOUNDARY, 
    };
    class FMM_DofmanRecord
    {
public:
        FNM_Status_Type status;
    };

    std :: vector< FMM_DofmanRecord >dmanRecords;
    const FloatArray *dmanValuesPtr;
    // Delegate of FMM_DofmanRecord; stored in priority queue.
    //class FMM_DofmanRecordDelegate {
    //public:FMM_DofmanRecordDelegaFMM_DofmanRecordDelegatete
    //  int id;
    //};

    struct FMM_DofmanRecordDelegate_greater {
        const FloatArray **dmanValuesPtrRef;
        FMM_DofmanRecordDelegate_greater(const FloatArray **dv) : dmanValuesPtrRef(dv)
        { }
        bool operator() (const int & p, const int & q) const
        { return ( fabs( ( * dmanValuesPtrRef )->at(p) ) > fabs( ( * dmanValuesPtrRef )->at(q) ) ); }
    };

    Domain *domain;

    std :: priority_queue< int, std :: vector< int > , FMM_DofmanRecordDelegate_greater >dmanTrialQueue;

public:
    FastMarchingMethod(Domain * d) : dmanTrialQueue( FMM_DofmanRecordDelegate_greater(& this->dmanValuesPtr) ) {
        domain = d;
    }
    ~FastMarchingMethod() { }

    void solve(FloatArray &dmanValues, const std :: list< int > &bcDofMans, double F);

    // identification
    const char *giveClassName() const { return "FastMarchingMethod"; }

protected:
    void initialize(FloatArray &dmanValues, const std :: list< int > &bcDofMans, double F);

    void updateTrialValue(FloatArray &dmanValues, int id, double F);

    int  getSmallestTrialDofMan();
};
} // end namespace oofem
#endif // fastmarchingmethod_h