From cb69732f68c0bd46c1574de16ce1aee6f38e439b Mon Sep 17 00:00:00 2001 From: Brett Weiland Date: Tue, 11 Jun 2024 14:50:14 -0500 Subject: restarting --- live_play/kernel.c | 40 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 40 insertions(+) create mode 100644 live_play/kernel.c (limited to 'live_play/kernel.c') 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 +#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; +} -- cgit v1.2.3