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
|
#include <complex.h>
#include <stdio.h>
#include <time.h>
#include <stdlib.h>
struct kompleksno {
long double r;
long double i;
};
struct kompleksno množi (struct kompleksno a, struct kompleksno b) {
struct kompleksno r;
r.r = (a.r * b.r) - (a.i * b.i /* ii = -1 */);
r.i = (a.i * b.r) + (a.r * b.i); /* zakon o distributivnosti */
return r;
}
struct kompleksno seštej (struct kompleksno a, struct kompleksno b) {
a.r += b.r;
a.i += b.i;
return a;
}
int main (int argc, char ** argv) {
long double complex c = 0.001 + 0.0001 * I;
struct kompleksno k = { .r = 0.001, .i = 0.0001 };
long double complex a = -0.5 + 0.55 * I;
struct kompleksno s = { .r = -0.5, .i = 0.55 };
int število = argc > 1 ? atoi(argv[1]) : 10000;
clock_t cikli = clock();
for (int i = 0; i < število; i++) {
c = c*c + a;
if (creall(c) > 0.01 && cimagl(c) < -0.02)
c += 0.00003;
}
fprintf(stderr, "rezultat complex.h: %Lf+%Lfi\n", creall(c), cimagl(c));
fprintf(stderr, " porabljen čas: %ld ciklov\n", clock()-cikli);
cikli = clock();
for (int i = 0; i < število; i++) {
k = seštej(množi(k, k), s);
if (k.r > 0.01 && k.i < -0.02)
k.r += 0.00003;
}
fprintf(stderr, "rezultat na roke: %Lf+%Lfi\n", k.r, k.i);
fprintf(stderr, " porabljen čas: %ld ciklov\n", clock()-cikli);
return 0;
}
|
