Blame view

Modif/epsilon-master/apps/regression/model/quartic_model.cpp 3.85 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
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
  #include "quartic_model.h"
  #include "../../shared/poincare_helpers.h"
  #include <math.h>
  #include <assert.h>
  #include "../../poincare/include/poincare_layouts.h"
  
  using namespace Poincare;
  using namespace Shared;
  
  namespace Regression {
  
  ExpressionLayout * QuarticModel::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('X', KDText::FontSize::Small),
        new VerticalOffsetLayout(
            new CharLayout('4', KDText::FontSize::Small),
            VerticalOffsetLayout::Type::Superscript,
            false),
        new CharLayout('+', KDText::FontSize::Small),
        new CharLayout('b', KDText::FontSize::Small),
        new CharLayout(Ion::Charset::MiddleDot, KDText::FontSize::Small),
        new CharLayout('X', KDText::FontSize::Small),
        new VerticalOffsetLayout(
            new CharLayout('3', KDText::FontSize::Small),
            VerticalOffsetLayout::Type::Superscript,
            false),
        new CharLayout('+', KDText::FontSize::Small),
        new CharLayout('c', KDText::FontSize::Small),
        new CharLayout(Ion::Charset::MiddleDot, KDText::FontSize::Small),
        new CharLayout('X', KDText::FontSize::Small),
        new VerticalOffsetLayout(
            new CharLayout('2', KDText::FontSize::Small),
            VerticalOffsetLayout::Type::Superscript,
            false),
        new CharLayout('+', KDText::FontSize::Small),
        new CharLayout('d', KDText::FontSize::Small),
        new CharLayout(Ion::Charset::MiddleDot, KDText::FontSize::Small),
        new CharLayout('X', KDText::FontSize::Small),
        new CharLayout('+', KDText::FontSize::Small),
        new CharLayout('e', KDText::FontSize::Small),
      };
      layout = new HorizontalLayout(layoutChildren, 20, false);
    }
    return layout;
  }
  
  Expression * QuarticModel::simplifiedExpression(double * modelCoefficients, Poincare::Context * context) {
    double a = modelCoefficients[0];
    double b = modelCoefficients[1];
    double c = modelCoefficients[2];
    double d = modelCoefficients[3];
    double e = modelCoefficients[4];
    Expression * ax4Expression = new Multiplication(
        new Decimal(a),
        new Power(
          new Symbol('x'),
          new Decimal(4.0),
          false),
        false);
    Expression * bx3Expression = new Multiplication(
        new Decimal(b),
        new Power(
          new Symbol('x'),
          new Decimal(3.0),
          false),
        false);
    Expression * cx2Expression = new Multiplication(
      new Decimal(c),
      new Power(
        new Symbol('x'),
        new Decimal(2.0),
        false),
      false);
    Expression * dxExpression = new Multiplication(
      new Decimal(d),
      new Symbol('x'),
      false);
    Expression * eExpression = new Decimal(e);
    Expression * const operands[] = {ax4Expression, bx3Expression, cx2Expression, dxExpression, eExpression};
    Expression * result = new Addition(operands, 5, false);
    PoincareHelpers::Simplify(&result, *context);
    return result;
  }
  
  double QuarticModel::evaluate(double * modelCoefficients, double x) const {
    double a = modelCoefficients[0];
    double b = modelCoefficients[1];
    double c = modelCoefficients[2];
    double d = modelCoefficients[3];
    double e = modelCoefficients[4];
    return a*x*x*x*x+b*x*x*x+c*x*x+d*x+e;
  }
  
  double QuarticModel::partialDerivate(double * modelCoefficients, int derivateCoefficientIndex, double x) const {
    if (derivateCoefficientIndex == 0) {
      // Derivate: x^4
      return x*x*x*x;
    }
    if (derivateCoefficientIndex == 1) {
      // Derivate: x^3
      return x*x*x;
    }
    if (derivateCoefficientIndex == 2) {
      // Derivate: x^2
      return x*x;
    }
    if (derivateCoefficientIndex == 3) {
      // Derivate: x
      return x;
    }
    if (derivateCoefficientIndex == 4) {
      // Derivate: 1
      return 1;
    }
    assert(false);
    return 0.0;
  }
  
  }