equal
deleted
inserted
replaced
95 |
95 |
96 // deltas for each grid step |
96 // deltas for each grid step |
97 val d_x = (end.re - start.re) / W |
97 val d_x = (end.re - start.re) / W |
98 val d_y = (end.im - start.im) / H |
98 val d_y = (end.im - start.im) / H |
99 |
99 |
100 for (y <- (0 until H)) { |
100 for (y <- (0 until H).par) { |
101 for (x <- (0 until W)) { |
101 for (x <- (0 until W).par) { |
102 |
102 |
103 val c = start + |
103 val c = start + |
104 (x * d_x + y * d_y * i) |
104 (x * d_x + y * d_y * i) |
105 val iters = iterations(c, max) |
105 val iters = iterations(c, max) |
106 val col = |
106 val col = |
146 |
146 |
147 // some more computations with example 3 |
147 // some more computations with example 3 |
148 |
148 |
149 val delta = (exc2 - exc1) * 0.0333 |
149 val delta = (exc2 - exc1) * 0.0333 |
150 |
150 |
151 |
151 /* |
152 time_needed( |
152 time_needed( |
153 for (n <- (0 to 12)) |
153 for (n <- (0 to 12)) |
154 mandelbrot(exc1 + delta * n, |
154 mandelbrot(exc1 + delta * n, |
155 exc2 - delta * n, 100)) |
155 exc2 - delta * n, 100)) |
156 |
156 */ |
157 /* |
157 /* |
158 time_needed( |
158 time_needed( |
159 for (n <- (0 to 12)) |
159 for (n <- (0 to 12)) |
160 mandelbrot(exc1 + delta * n, |
160 mandelbrot(exc1 + delta * n, |
161 exc2 - delta * n, 1000)) |
161 exc2 - delta * n, 1000)) |
162 */ |
162 */ |
163 |
163 |
164 |
164 |
|
165 // Larry Paulson's example |
|
166 val exl1 = -0.74364990 + 0.13188170 * i |
|
167 val exl2 = -0.74291189 + 0.13261971 * i |
|
168 |
|
169 time_needed(mandelbrot(exl1, exl2, 1000)) |
|
170 |
|
171 |
|
172 // example by Jorgen Villadsen |
|
173 val exj1 = 0.10284 - 0.63275 * i |
|
174 val exj2 = 0.11084 - 0.64075 * i |
|
175 |
|
176 time_needed(mandelbrot(exj1, exj2, 1000)) |