mirror/yacreader
1
0
Fork 0
mirror of https://github.com/YACReader/yacreader synced 2025-06-04 01:28:55 -04:00
Code Issues Actions Packages Projects Releases Wiki Activity

More portable code

Browse Source
This commit is contained in:
Luis Ángel San Martín Rodríguez 2024-08-31 15:20:37 +02:00
parent eacbd5f499
commit 9fc3ab4981
1 changed files with 30 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
image_processing/resize_image.cpp
View File

@ -195,17 +195,45 @@ QVector<QRgb> processRow(int y, int targetWidth, int targetHeight, const QImage
inline float64x2_t lanczosKernelNEON(const float64x2_t &x, int a) inline float64x2_t lanczosKernelNEON(const float64x2_t &x, int a)
{ {
// Load constants
float64x2_t zero = vdupq_n_f64(0.0); float64x2_t zero = vdupq_n_f64(0.0);
float64x2_t one = vdupq_n_f64(1.0); float64x2_t one = vdupq_n_f64(1.0);
float64x2_t a_val = vdupq_n_f64(a); float64x2_t a_val = vdupq_n_f64(a);
// Convert to scalar arrays
double x_array[2];
vst1q_f64(x_array, x);
double a_val_array[2] = { static_cast<double>(a), static_cast<double>(a) };
// Compute sin(x * pi)
float64x2_t pix = vmulq_f64(vdupq_n_f64(M_PI), x); float64x2_t pix = vmulq_f64(vdupq_n_f64(M_PI), x);
float64x2_t sin_pix = vsin_f64(pix); double pix_array[2];
float64x2_t sin_pix_a = vsin_f64(vdivq_f64(pix, a_val)); vst1q_f64(pix_array, pix);
double sin_pix_array[2];
sin_pix_array[0] = std::sin(pix_array[0]);
sin_pix_array[1] = std::sin(pix_array[1]);
float64x2_t sin_pix = vld1q_f64(sin_pix_array);
// Compute sin(x * pi / a)
float64x2_t pix_div_a = vdivq_f64(pix, a_val);
double pix_div_a_array[2];
vst1q_f64(pix_div_a_array, pix_div_a);
double sin_pix_div_a_array[2];
sin_pix_div_a_array[0] = std::sin(pix_div_a_array[0]);
sin_pix_div_a_array[1] = std::sin(pix_div_a_array[1]);
float64x2_t sin_pix_a = vld1q_f64(sin_pix_div_a_array);
// Compute Lanczos kernel
float64x2_t numerator = vmulq_f64(vmulq_f64(a_val, sin_pix), sin_pix_a); float64x2_t numerator = vmulq_f64(vmulq_f64(a_val, sin_pix), sin_pix_a);
float64x2_t denominator = vmulq_f64(pix, pix); float64x2_t denominator = vmulq_f64(pix, pix);
float64x2_t result = vdivq_f64(numerator, denominator); float64x2_t result = vdivq_f64(numerator, denominator);
// Handle the case where x is zero
uint64x2_t mask = vceqq_f64(x, zero); uint64x2_t mask = vceqq_f64(x, zero);
result = vbslq_f64(mask, one, result); result = vbslq_f64(mask, one, result);
return result; return result;
} }