diff options
author | Brett Weiland <brett_weiland@gmail.com> | 2024-06-11 14:50:14 -0500 |
---|---|---|
committer | Brett Weiland <brett_weiland@gmail.com> | 2024-06-11 14:50:14 -0500 |
commit | cb69732f68c0bd46c1574de16ce1aee6f38e439b (patch) | |
tree | def1daaec81a0d4cd7b3d44b2c26b9535e07579c /live_play/kernel.c |
Diffstat (limited to 'live_play/kernel.c')
-rw-r--r-- | live_play/kernel.c | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/live_play/kernel.c b/live_play/kernel.c new file mode 100644 index 0000000..e572a48 --- /dev/null +++ b/live_play/kernel.c @@ -0,0 +1,40 @@ +//#include <math.h> +#define PI 3.141592653589793115997963468544185161590576171875 + +double cosecant_single(double a, double b) { return a / sin(b); } +double secant_single(double a, double b) { + //double killme = (double)floor(b / (PI / 2.0)) * (double)(PI / 2.0); + //return a / sin(killme); + //return a / sin((PI / 2.0)); + //return a / sin(PI); + return a / sin(-PI); +} + + +__kernel void render_frame(__global unsigned int *frame_output, __global double *mask, + double x_step, double y_step, + double x_start, double y_start, + unsigned int iterations, unsigned int escape, double ratio) { + unsigned int result; + double x = (get_global_id(0) * x_step) + x_start; + double y = (get_global_id(1) * y_step) + y_start; + size_t img_index = (get_global_id(1) * get_global_size(1)) + get_global_id(0); + + unsigned int iter; + + + + for(iter = 0; iter < iterations; iter++) { + double next_x; + double r = mask[img_index]; + next_x = (r * cosecant_single(x, y)) + ((1 - r) * secant_single(x, y)); + y = (r * cosecant_single(y, x)) + ((1 - r) * secant_single(y, x)); + x = next_x; + if((pow(x, 2) + pow(y, 2)) >= escape) break; + } + + + + frame_output[img_index] = iter; + //frame_output[img_index] = mask[img_index] * 255; +} |