Blame view

emulateur/epsilon-nofrendo/apps/regression/model/logarithmic_model.cpp 2.02 KB
6663b6c9   adorian   projet complet av...
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
66
67
68
69
70
71
72
  #include "logarithmic_model.h"
  #include "../store.h"
  #include "../../poincare/include/poincare_layouts.h"
  #include <math.h>
  #include <assert.h>
  
  using namespace Poincare;
  
  namespace Regression {
  
  ExpressionLayout * LogarithmicModel::layout() {
    static ExpressionLayout * layout = nullptr;
    if (layout == nullptr) {
      const ExpressionLayout * layoutChildren[] = {
        new CharLayout('a', KDText::FontSize::Small),
        new CharLayout(Ion::Charset::MiddleDot, KDText::FontSize::Small),
        new CharLayout('l', KDText::FontSize::Small),
        new CharLayout('n', KDText::FontSize::Small),
        new CharLayout('(', KDText::FontSize::Small),
        new CharLayout('X', KDText::FontSize::Small),
        new CharLayout(')', KDText::FontSize::Small),
        new CharLayout('+', KDText::FontSize::Small),
        new CharLayout('b', KDText::FontSize::Small)
      };
      layout = new HorizontalLayout(layoutChildren, 9, false);
    }
    return layout;
  }
  
  double LogarithmicModel::evaluate(double * modelCoefficients, double x) const {
    double a = modelCoefficients[0];
    double b = modelCoefficients[1];
    return a*log(x)+b;
  }
  
  double LogarithmicModel::levelSet(double * modelCoefficients, double xMin, double step, double xMax, double y, Poincare::Context * context) {
    double a = modelCoefficients[0];
    double b = modelCoefficients[1];
    if (a == 0) {
      return NAN;
    }
    return exp((y-b)/a);
  }
  
  double LogarithmicModel::partialDerivate(double * modelCoefficients, int derivateCoefficientIndex, double x) const {
    if (derivateCoefficientIndex == 0) {
      // Derivate: ln(x)
      assert(x >0);
      return log(x);
    }
    if (derivateCoefficientIndex == 1) {
      // Derivate: 1
      return 1;
    }
    assert(false);
    return 0.0;
  }
  
  bool LogarithmicModel::dataSuitableForFit(Store * store, int series) const {
    if (!Model::dataSuitableForFit(store, series)) {
      return false;
    }
    int numberOfPairs = store->numberOfPairsOfSeries(series);
    for (int j = 0; j < numberOfPairs; j++) {
      if (store->get(series, 0, j) <= 0) {
        return false;
      }
    }
    return true;
  }
  
  }