globe/fc.cpp at main · shtuk41/globe · 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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
#include "fc.h"

CTable::CTable() : input(0), output(0),numberPoints(0)
{

}

template<typename... Args>
void CTable::Init(int n_points, Args... args)
{
   numberPoints = n_points;
   input = new double[n_points];
   output = new double[n_points];

   va_list ap;
   va_start(ap, n_points);

   double values[] = { args... };

   for(int ii = 0; ii < n_points; ii++) {
       input[ii]  = values[ii * 2];  // first value
       output[ii] = values[ii * 2 + 1];  // second value
   }
}

double CTable::operator<<(double in)
 {
    double out = output[0];

    if (in <= input[0])
    {
        return out;
    }
    else if (in > input[0] && in < input[numberPoints - 1])
    {
        for (int ii = 0; ii < numberPoints - 1; ii++)
        {
            if (in > input[ii] && in <= input[ii+1])
            {
                out = output[ii] + (in - input[ii]) / (input[ii + 1] - input[ii]) * (output[ii + 1] - output[ii]);
                return out;
            }
        }
    }
    else
    {
        out = output[numberPoints - 1];
        return out;
    }
 }

 CTable::~CTable()
 {
     if (input != 0)
     {
         delete [] input;
         input = 0;
     }

     if (output != 0)
     {
         delete [] output;
         output = 0;
     }
 }