-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathfractals_2.c
113 lines (104 loc) · 2.96 KB
/
fractals_2.c
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
/* ************************************************************************** */
/* */
/* ::: :::::::: */
/* fractals_2.c :+: :+: :+: */
/* +:+ +:+ +:+ */
/* By: jmousset <marvin@42.fr> +#+ +:+ +#+ */
/* +#+#+#+#+#+ +#+ */
/* Created: 2019/10/10 12:00:46 by jmousset #+# #+# */
/* Updated: 2019/10/14 18:26:25 by jmousset ### ########.fr */
/* */
/* ************************************************************************** */
#include "fractol.h"
void celtic_mandelbar(t_env *env, t_coord p)
{
t_complex old;
t_complex new;
t_complex c;
int i;
i = 0;
old.r = 0;
old.i = 0;
new.r = 0;
new.i = 0;
c.r = 1.5 * (p.x - WIDTH / 2) / (0.5 * env->zoom * WIDTH) + env->x;
c.i = 1 * (p.y - HEIGHT / 2) / (0.5 * env->zoom * HEIGHT) + env->y;
while ((i < env->max_it) && (new.r * new.r + new.i * new.i < 4))
{
old.r = new.r;
old.i = new.i;
new.r = fabs(old.r * old.r - old.i * old.i) + c.r;
new.i = -2 * old.r * old.i + c.i;
i++;
}
put_pixel(env, p, i, new);
}
void celtic_mandelbrot(t_env *env, t_coord p)
{
t_complex old;
t_complex new;
t_complex c;
int i;
i = 0;
old.r = 0;
old.i = 0;
new.r = 0;
new.i = 0;
c.r = 2 * (p.x - WIDTH / 2) / (0.5 * env->zoom * WIDTH) + env->x;
c.i = 1.5 * (p.y - HEIGHT / 2) / (0.5 * env->zoom * HEIGHT) + env->y;
while ((i < env->max_it) && (new.r * new.r + new.i * new.i < 4))
{
old.r = new.r;
old.i = new.i;
new.r = fabs(old.r * old.r - old.i * old.i) + c.r;
new.i = 2 * old.r * old.i + c.i;
i++;
}
put_pixel(env, p, i, new);
}
void buffalo(t_env *env, t_coord p)
{
t_complex old;
t_complex new;
t_complex c;
int i;
i = 0;
old.r = 0;
old.i = 0;
new.r = 0;
new.i = 0;
c.r = 2 * (p.x - WIDTH / 2) / (0.5 * env->zoom * WIDTH) + env->x;
c.i = 1.5 * (p.y - HEIGHT / 2) / (0.5 * env->zoom * HEIGHT) + env->y;
while ((i < env->max_it) && (new.r * new.r + new.i * new.i < 4))
{
old.r = new.r;
old.i = new.i;
new.r = fabs(old.r * old.r - old.i * old.i) + c.r;
new.i = 2 * fabs(old.r * old.i) + c.i;
i++;
}
put_pixel(env, p, i, new);
}
void celtic_heart(t_env *env, t_coord p)
{
t_complex old;
t_complex new;
t_complex c;
int i;
i = 0;
old.r = 0;
old.i = 0;
new.r = 0;
new.i = 0;
c.r = 2 * (p.x - WIDTH / 2) / (0.5 * env->zoom * WIDTH) + env->x;
c.i = 1.5 * (p.y - HEIGHT / 2) / (0.5 * env->zoom * HEIGHT) + env->y;
while ((i < env->max_it) && (new.r * new.r + new.i * new.i < 4))
{
old.r = new.r;
old.i = new.i;
new.r = fabs(old.r * old.r - old.i * old.i) + c.r;
new.i = 2 * fabs(old.r) * old.i + c.i;
i++;
}
put_pixel(env, p, i, new);
}