-
Notifications
You must be signed in to change notification settings - Fork 1
/
elliptic_test.cpp
70 lines (54 loc) · 1.28 KB
/
elliptic_test.cpp
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
#include <cmath>
#include <iostream>
#include <cstdlib>
#include <string>
#include "complex_plot.hpp"
#include "elliptic_function01.hpp"
using namespace std;
class elliptic_test01 : public complex_plot {
protected:
elliptic_function E;
double f;
cplx transform(cplx z) {
return E.compute_raw(z)*f;
}
public:
elliptic_test01() : E(exp(cplx(0, 0.5*M_PI/3)),
exp(cplx(0, -0.5*M_PI/3)))
{
E.n_ = 4;
f = 1.0/abs(E.a_);
set_title("Elliptic function approximation");
bmp_name = "ell_rat02.bmp";
width_set_window(1024, 1024, 16);
run(2);
}
};
class elliptic_test02 : public complex_plot {
protected:
elliptic_function E;
double f;
cplx transform(cplx z) {
return E.compute_raw(z)*f;
}
public:
elliptic_test02() : E(exp(cplx(0, 0.0*M_PI/3)),
3.0*exp(cplx(0, 0.5*M_PI/3)))
{
E.n_ = 6;
f = 1.0/abs(E.a_);
set_title("Elliptic function approximation");
bmp_name = "ell_rat03.bmp";
width_set_window(1024, 1024, 64);
run(2);
}
};
int main(int argc, char **argv) {
try {
elliptic_test01 E;
} catch (std::string s) {
cerr << "Error: " << s << endl;
return 1;
}
return 0;
}