void Calc Initalfa std vector MapNode alfavector for unsigned int alfa

 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
void Calc::Initalfa()
{
std::vector<MapNode*> alfavector;
for(unsigned int i = 0; i < alfa0.size()-1; i++)
{
MapNode *alfa = new MapNode();
mpfr_div (alfa->value, alfa0[i+1]->value, alfa0[i]->value, MPFR_RNDD);
alfavector.push_back(alfa);
// std::cout<<std::endl<<"it" << i <<"alfa:";
// std::cout<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfa0[i+1]->value, MPFR_RNDD);
// std::cout<<" /"<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfa0[i]->value, MPFR_RNDD);
// std::cout<<" ="<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfa->value, MPFR_RNDD);
}
alfas.push_back(alfavector);
alfavector.clear();
for(unsigned int i = 0; i < alfas[1].size()-1; i++)
{
MapNode *alfa = new MapNode();
mpfr_neg (alfa->value, alfas[1][i+1]->value, MPFR_RNDD);
mpfr_add (alfa->value, alfa->value, alfas[1][i]->value ,MPFR_RNDD);
alfavector.push_back(alfa);
}
alfas.push_back(alfavector);
alfavector.clear();
for(unsigned int i = 3; i < alfa0.size(); i++)
{
if(i%2==0)
{
for(unsigned int j = 0; j < alfas[i-1].size()-1; j++)
{
MapNode *alfa = new MapNode();
//mpfr_neg (alfa->value, alfas[i-2][j+1]->value, MPFR_RNDD);
mpfr_sub (alfa->value,alfas[i-2][j+1]->value, alfas[i-1][j+1]->value, MPFR_RNDD);
mpfr_add (alfa->value, alfa->value, alfas[i-1][j]->value , MPFR_RNDD);
alfavector.push_back(alfa);
// std::cout<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfas[i-2][j+1]->value, MPFR_RNDD);
// std::cout<<" - ";
// mpfr_out_str(stdout, 10, 10, alfas[i-1][j+1]->value, MPFR_RNDD);
// std::cout<<" + "<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfas[i-1][j]->value, MPFR_RNDD);
// std::cout<<" = "<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfa->value, MPFR_RNDD);
// std::cout<<std::endl;
}
}
else
{
for(unsigned int j = 0; j < alfas[i-1].size()-1; j++)
{
MapNode *alfa = new MapNode();
mpfr_mul (alfa->value, alfas[i-1][j+1]->value, alfas[i-2][j+1]->value, MPFR_RNDD);
mpfr_div (alfa->value, alfa->value, alfas[i-1][j]->value , MPFR_RNDD);
alfavector.push_back(alfa);
// std::cout<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfa->value, MPFR_RNDD);
// std::cout<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfas[i-2][j+1]->value, MPFR_RNDD);
// std::cout<<" * ";
// mpfr_out_str(stdout, 10, 10, alfas[i-1][j+1]->value, MPFR_RNDD);
// std::cout<<" / "<<std::endl;
// mpfr_out_str(stdout, 10, 10, alfas[i-1][j]->value, MPFR_RNDD);
// std::cout<<" = ";
// mpfr_out_str(stdout, 10, 10, alfa->value, MPFR_RNDD);
// std::cout<<std::endl;
}
}
alfas.push_back(alfavector);
alfavector.clear();
}
}