FasterCard/FastMarching.hpp at main · cesare-corrado/FasterCard · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
#ifndef FMARCHING_HPP
#define FMARCHING_HPP

#include<stdio.h>
#include<iostream>
#include<string>
#include<vector>
#include<map>
#include<set>
#include "Mesh.hpp"
#include "DataEikonal.hpp"
#include "Eikonal.hpp"


class FastMarching: public Eikonal
{
      public:
            FastMarching(const Mesh * theMesh, const DataEikonal * theData);
            ~FastMarching();
            void iterate(const double & integration_time);
            void evaluateActivationTimes();

      private:
            void initializeAuxVar();
            void updateSubset(const PointConnect & neighborhood, const PointConnect & subset,const size_t & pn);
            void updateNeighPoint(const IndexVector & pointList);
            double timeij(const IndexVector & pointList);
            double ceij(const IndexVector & pointList);
            RealVector y_roots(const double & a, const double &b, const double &c ,const double &d, const double &e, const double & f,const double & DT01, const double & DT02);
            RealVector f1_x_roots_y(const double & a, const double &b, const double &c ,const double &d, const double &e, const double & f,const double & DT01, const double & DT02, const double & yr);
            RealVector f1_y_roots_x(const double & a, const double &b, const double &c ,const double &d, const double &e, const double & f,const double & DT01, const double & DT02, const double & xr);

            RealVector f2_x_roots_y(const double & a, const double &b, const double &c ,const double &d, const double &e, const double & f,const double & DT01, const double & DT02, const double & yr);
            RealVector f2_y_roots_x(const double & a, const double &b, const double &c ,const double &d, const double &e, const double & f,const double & DT01, const double & DT02, const double & xr);


            double compute_delta_4ord(const double & a,const double & b,const double & c,const double & d,const double & e);
            double compute_D_4ord(const double & a,const double & b,const double & c,const double & d,const double & e);
            double compute_delta_0_4ord(const double & a,const double & b,const double & c,const double & d,const double & e);
            double compute_delta_1_4ord(const double & a,const double & b,const double & c,const double & d,const double & e);
            // FMM has 4 states: trial, known, changed and unknown
            // It requires 2 queues: one for trial, the other for changes
            QueueType _changed;
            BoolVector _ischanged;
            queuedPtrType _changedPtr;
};

#endif