Julia

[IgorErin]

Done:

• host ref x86
• halide naive

[dkurt]

Add "${binary_dir}/install/lib/libopencv_imgcodecs.so" to target_link_libraries(algos instead:

--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -92,6 +92,7 @@ target_include_directories(algos PRIVATE
 target_link_libraries(algos
   "${binary_dir}/install/lib/libopencv_core.so"
   "${binary_dir}/install/lib/libopencv_imgproc.so"
+  "${binary_dir}/install/lib/libopencv_imgcodecs.so"
 )

[IgorErin]

@dkurt Добавил наивную реализацю на halide. На моем x86 работает заметно медленнее референсной. Cейчас буду мерить, заодно подготовлю бенчмарки на завтра. Или есть задачи актуальнее?

Стала ли лучше картинка?

[dkurt]

@IgorErin, да, всё супер! Если увеличить julia_width/julia_height то побольше картинка, всё рисует правильно. Да, пока медленнее референс части, но предлагаю сейчас PR довести до зелёного состояния и мерджить, работа сделана уже большая.

[dkurt]

Suggested change

	void halide_julia(uint8_t* dst, int height, int width) {
	static const uint8_t step_bound = 255;
	static const uint8_t upper_bound = 4;
	Buffer output(dst, width, height);
	Var x, y;
	Complex c(cr, ci);
	Expr x_ranged;
	x_ranged = halide_julia_norm(x, height);
	Expr y_ranged;
	y_ranged = halide_julia_norm(y, width);
	Func julia;
	Var t;
	julia(x, y, t) = Complex(x_ranged, y_ranged);
	// loop
	RDom index(1, step_bound);
	Complex current = julia(x, y, index - 1);
	julia(x, y, index) = current * current + c;
	Expr esc_cond = Complex(julia(x, y, index)).magnitude_squared() < upper_bound;
	Tuple first_escape = argmin(esc_cond);
	// proj to result
	Func result;
	result(x, y) = cast<uint8_t>(first_escape[0]);
	result.realize(output);
	}
	void halide_julia(uint8_t* dst, int height, int width) {
	#ifdef __riscv
	// add later
	#else
	static const uint8_t step_bound = 255;
	static const uint8_t upper_bound = 4;
	Buffer output(dst, width, height);
	Var x, y;
	Complex c(cr, ci);
	Expr x_ranged;
	x_ranged = halide_julia_norm(x, height);
	Expr y_ranged;
	y_ranged = halide_julia_norm(y, width);
	Func julia;
	Var t;
	julia(x, y, t) = Complex(x_ranged, y_ranged);
	// loop
	RDom index(1, step_bound);
	Complex current = julia(x, y, index - 1);
	julia(x, y, index) = current * current + c;
	Expr esc_cond = Complex(julia(x, y, index)).magnitude_squared() < upper_bound;
	Tuple first_escape = argmin(esc_cond);
	// proj to result
	Func result;
	result(x, y) = cast<uint8_t>(first_escape[0]);
	result.realize(output);
	#endif
	}