int julia(int x, int y) { const float scale = 1.5F; float jx = scale * (float)(DIM / 2 - x) / (DIM / 2); float jy = scale * (float)(DIM / 2 - y) / (DIM / 2); cuComplex c = new cuComplex(-0.8F, 0.156F); cuComplex a = new cuComplex(jx, jy); int i = 0; for (i = 0; i < 200; i++) { a = a * a + c; if (a.magnitude2() > 1000) { return(0); } } return(1); }
int julia(int x, int y) { const float scale = 1.5F; float jx = scale * (float)(DIM / 2 - x) / (DIM / 2); float jy = scale * (float)(DIM / 2 - y) / (DIM / 2); cuComplex c = new cuComplex(-0.8F, 0.156F); cuComplex a = new cuComplex(jx, jy); int i = 0; for (i = 0; i < 200; i++) { a = a * a + c; if (a.magnitude2() > 1000) return 0; } return 1; }