jackaudio
/
jack2
mirror of https://github.com/jackaudio/jack2
1
0
Fork 0
Code Releases 45 Activity
jack2 codebase
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
3601 Commits
8 Branches
255MB
Tree: a82f3f2fb4
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'a82f3f2fb4'
${ noResults }
jack2/common/memops.c

1172 lines
32KB
Raw Blame History 

  1. /*

  2.     Copyright (C) 2000 Paul Davis 

  3. 

  4.     This program is free software; you can redistribute it and/or modify

  5.     it under the terms of the GNU General Public License as published by

  6.     the Free Software Foundation; either version 2 of the License, or

  7.     (at your option) any later version.

  8. 

  9.     This program is distributed in the hope that it will be useful,

  10.     but WITHOUT ANY WARRANTY; without even the implied warranty of

  11.     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

  12.     GNU General Public License for more details.

  13. 

  14.     You should have received a copy of the GNU General Public License

  15.     along with this program; if not, write to the Free Software

  16.     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

  17. 

  18. */

  19. 

  20. #define _ISOC9X_SOURCE  1

  21. #define _ISOC99_SOURCE  1

  22. 

  23. #define __USE_ISOC9X    1

  24. #define __USE_ISOC99    1

  25. 

  26. #include <stdio.h>

  27. #include <string.h>

  28. #include <math.h>

  29. #include <memory.h>

  30. #include <stdlib.h>

  31. #include <stdint.h>

  32. #include <limits.h>

  33. #ifdef __linux__

  34. #include <endian.h>

  35. #endif

  36. #include "memops.h"

  37. 

  38. #if defined (__SSE2__) && !defined (__sun__)

  39. #include <emmintrin.h>

  40. #ifdef __SSE4_1__

  41. #include <smmintrin.h>

  42. #endif

  43. #endif

  44. 

  45. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  46. #include <arm_neon.h>

  47. #endif

  48. 

  49. /* Notes about these *_SCALING values.

  50. 

  51.    the MAX_<N>BIT values are floating point. when multiplied by

  52.    a full-scale normalized floating point sample value (-1.0..+1.0)

  53.    they should give the maxium value representable with an integer

  54.    sample type of N bits. Note that this is asymmetric. Sample ranges 

  55.    for signed integer, 2's complement values are -(2^(N-1) to +(2^(N-1)-1)

  56. 

  57.    Complications

  58.    -------------

  59.    If we use +2^(N-1) for the scaling factors, we run into a problem:

  60. 

  61.    if we start with a normalized float value of -1.0, scaling

  62.    to 24 bits would give -8388608 (-2^23), which is ideal.

  63.    But with +1.0, we get +8388608, which is technically out of range.

  64. 

  65.    We never multiply a full range normalized value by this constant,

  66.    but we could multiply it by a positive value that is close enough to +1.0

  67.    to produce a value > +(2^(N-1)-1.

  68. 

  69.    There is no way around this paradox without wasting CPU cycles to determine

  70.    which scaling factor to use (i.e. determine if its negative or not,

  71.    use the right factor).

  72. 

  73.    So, for now (October 2008) we use 2^(N-1)-1 as the scaling factor.

  74. */

  75. 

  76. #define SAMPLE_24BIT_SCALING  8388607

  77. #define SAMPLE_16BIT_SCALING  32767

  78. 

  79. /* these are just values to use if the floating point value was out of range

  80.    

  81.    advice from Fons Adriaensen: make the limits symmetrical

  82.  */

  83. 

  84. #define SAMPLE_24BIT_MAX  8388607  

  85. #define SAMPLE_24BIT_MIN  -8388607 

  86. #define SAMPLE_24BIT_MAX_F  8388607.0f  

  87. #define SAMPLE_24BIT_MIN_F  -8388607.0f 

  88. 

  89. #define SAMPLE_16BIT_MAX  32767

  90. #define SAMPLE_16BIT_MIN  -32767

  91. #define SAMPLE_16BIT_MAX_F  32767.0f

  92. #define SAMPLE_16BIT_MIN_F  -32767.0f

  93. 

  94. /* these mark the outer edges of the range considered "within" range

  95.    for a floating point sample value. values outside (and on the boundaries) 

  96.    of this range will be clipped before conversion; values within this 

  97.    range will be scaled to appropriate values for the target sample

  98.    type.

  99. */

  100. 

  101. #define NORMALIZED_FLOAT_MIN -1.0f

  102. #define NORMALIZED_FLOAT_MAX  1.0f

  103. 

  104. /* define this in case we end up on a platform that is missing

  105.    the real lrintf functions

  106. */

  107. 

  108. #define f_round(f) lrintf(f)

  109. 

  110. #define float_16(s, d)\

  111. 	if ((s) <= NORMALIZED_FLOAT_MIN) {\

  112. 		(d) = SAMPLE_16BIT_MIN;\

  113. 	} else if ((s) >= NORMALIZED_FLOAT_MAX) {\

  114. 		(d) = SAMPLE_16BIT_MAX;\

  115. 	} else {\

  116. 		(d) = f_round ((s) * SAMPLE_16BIT_SCALING);\

  117. 	}

  118. 

  119. /* call this when "s" has already been scaled (e.g. when dithering)

  120.  */

  121. 

  122. #define float_16_scaled(s, d)\

  123.         if ((s) <= SAMPLE_16BIT_MIN_F) {\

  124. 		(d) = SAMPLE_16BIT_MIN_F;\

  125. 	} else if ((s) >= SAMPLE_16BIT_MAX_F) {	\

  126. 		(d) = SAMPLE_16BIT_MAX;\

  127. 	} else {\

  128. 	        (d) = f_round ((s));\

  129. 	}

  130. 

  131. #define float_24u32(s, d) \

  132. 	if ((s) <= NORMALIZED_FLOAT_MIN) {\

  133. 		(d) = SAMPLE_24BIT_MIN;\

  134. 	} else if ((s) >= NORMALIZED_FLOAT_MAX) {\

  135. 		(d) = SAMPLE_24BIT_MAX;\

  136. 	} else {\

  137. 		(d) = f_round ((s) * SAMPLE_24BIT_SCALING);\

  138. 	}

  139. 

  140. 

  141. #define float_24(s, d) \

  142. 	if ((s) <= NORMALIZED_FLOAT_MIN) {\

  143. 		(d) = SAMPLE_24BIT_MIN;\

  144. 	} else if ((s) >= NORMALIZED_FLOAT_MAX) {\

  145. 		(d) = SAMPLE_24BIT_MAX;\

  146. 	} else {\

  147. 		(d) = f_round ((s) * SAMPLE_24BIT_SCALING);\

  148. 	}

  149. 

  150. 

  151. #if defined (__SSE2__) && !defined (__sun__)

  152. 

  153. /* generates same as _mm_set_ps(1.f, 1.f, 1f., 1f) but faster  */

  154. static inline __m128 gen_one(void)

  155. {

  156.     volatile __m128i x = { 0 }; /* shut up, GCC */

  157.     __m128i ones = _mm_cmpeq_epi32(x, x);

  158.     return (__m128)_mm_slli_epi32 (_mm_srli_epi32(ones, 25), 23);

  159. }

  160. 

  161. static inline __m128 clip(__m128 s, __m128 min, __m128 max)

  162. {

  163.     return _mm_min_ps(max, _mm_max_ps(s, min));

  164. }

  165. 

  166. static inline __m128i float_24_sse(__m128 s)

  167. {

  168.     const __m128 upper_bound = gen_one(); /* NORMALIZED_FLOAT_MAX */

  169.     const __m128 lower_bound = _mm_sub_ps(_mm_setzero_ps(), upper_bound);

  170. 

  171.     __m128 clipped = clip(s, lower_bound, upper_bound);

  172.     __m128 scaled = _mm_mul_ps(clipped, _mm_set1_ps(SAMPLE_24BIT_SCALING));

  173.     return _mm_cvtps_epi32(scaled);

  174. }

  175. #endif

  176. 

  177. 

  178. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  179. 

  180. static inline float32x4_t clip(float32x4_t s, float32x4_t min, float32x4_t max)

  181. {

  182. 	return vminq_f32(max, vmaxq_f32(s, min));

  183. }

  184. 

  185. static inline int32x4_t float_24_neon(float32x4_t s)

  186. {

  187. 	const float32x4_t upper_bound = vdupq_n_f32(NORMALIZED_FLOAT_MAX);

  188. 	const float32x4_t lower_bound = vdupq_n_f32(NORMALIZED_FLOAT_MIN);

  189. 

  190. 	float32x4_t clipped = clip(s, lower_bound, upper_bound);

  191. 	float32x4_t scaled = vmulq_f32(clipped, vdupq_n_f32(SAMPLE_24BIT_SCALING));

  192. 	return vcvtq_s32_f32(scaled);

  193. }

  194. 

  195. static inline int16x4_t float_16_neon(float32x4_t s)

  196. {

  197. 	const float32x4_t upper_bound = vdupq_n_f32(NORMALIZED_FLOAT_MAX);

  198. 	const float32x4_t lower_bound = vdupq_n_f32(NORMALIZED_FLOAT_MIN);

  199. 

  200. 	float32x4_t clipped = clip(s, lower_bound, upper_bound);

  201. 	float32x4_t scaled = vmulq_f32(clipped, vdupq_n_f32(SAMPLE_16BIT_SCALING));

  202. 	return vmovn_s32(vcvtq_s32_f32(scaled));

  203. }

  204. #endif

  205. 

  206. /* Linear Congruential noise generator. From the music-dsp list

  207.  * less random than rand(), but good enough and 10x faster 

  208.  */

  209. static unsigned int seed = 22222;

  210. 

  211. static inline unsigned int fast_rand() {

  212. 	seed = (seed * 196314165) + 907633515;

  213. 	return seed;

  214. }

  215. 

  216. /* functions for native float sample data */

  217. 

  218. void sample_move_floatLE_sSs (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip) {

  219. 	while (nsamples--) {

  220. 		*dst = *((float *) src);

  221. 		dst++;

  222. 		src += src_skip;

  223. 	}

  224. }

  225. 

  226. void sample_move_dS_floatLE (char *dst, jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state) {

  227. 	while (nsamples--) {

  228. 		*((float *) dst) = *src;

  229. 		dst += dst_skip;

  230. 		src++;

  231. 	}

  232. }

  233. 

  234. /* NOTES on function naming:

  235. 

  236.    foo_bar_d<TYPE>_s<TYPE>

  237. 

  238.    the "d<TYPE>" component defines the destination type for the operation

  239.    the "s<TYPE>" component defines the source type for the operation

  240. 

  241.    TYPE can be one of:

  242.    

  243.    S      - sample is a jack_default_audio_sample_t, currently (October 2008) a 32 bit floating point value

  244.    Ss     - like S but reverse endian from the host CPU

  245.    32u24  - sample is an signed 32 bit integer value, but data is in lower 24 bits only

  246.    32u24s - like 32u24 but reverse endian from the host CPU

  247.    24     - sample is an signed 24 bit integer value

  248.    24s    - like 24 but reverse endian from the host CPU

  249.    16     - sample is an signed 16 bit integer value

  250.    16s    - like 16 but reverse endian from the host CPU

  251. 

  252.    For obvious reasons, the reverse endian versions only show as source types.

  253. 

  254.    This covers all known sample formats at 16 bits or larger.

  255. */   

  256. 

  257. /* functions for native integer sample data */

  258. 

  259. void sample_move_d32u24_sSs (char *dst, jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)

  260. {

  261. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  262. 	unsigned long unrolled = nsamples / 4;

  263. 	nsamples = nsamples & 3;

  264. 

  265. 	while (unrolled--) {

  266. 		float32x4_t samples = vld1q_f32(src);

  267. 		int32x4_t converted = float_24_neon(samples);

  268. 		converted = vreinterpretq_s32_u8(vrev32q_u8(vreinterpretq_u8_s32(converted)));

  269. 

  270. 		switch(dst_skip) {

  271. 			case 4:

  272. 				vst1q_s32((int32_t*)dst, converted);

  273. 				break;

  274. 			default:

  275. 				vst1q_lane_s32((int32_t*)(dst),            converted, 0);

  276. 				vst1q_lane_s32((int32_t*)(dst+dst_skip),   converted, 1);

  277. 				vst1q_lane_s32((int32_t*)(dst+2*dst_skip), converted, 2);

  278. 				vst1q_lane_s32((int32_t*)(dst+3*dst_skip), converted, 3);

  279. 				break;

  280. 		}

  281. 		dst += 4*dst_skip;

  282. 		src+= 4;

  283. 	}

  284. #endif

  285. 

  286. 	int32_t z;

  287. 

  288. 	while (nsamples--) {

  289. 

  290. 		float_24u32 (*src, z);

  291. 

  292. #if __BYTE_ORDER == __LITTLE_ENDIAN

  293. 		dst[0]=(char)(z>>24);

  294. 		dst[1]=(char)(z>>16);

  295. 		dst[2]=(char)(z>>8);

  296. 		dst[3]=(char)(z);

  297. #elif __BYTE_ORDER == __BIG_ENDIAN

  298. 		dst[0]=(char)(z);

  299. 		dst[1]=(char)(z>>8);

  300. 		dst[2]=(char)(z>>16);

  301. 		dst[3]=(char)(z>>24);

  302. #endif

  303. 		dst += dst_skip;

  304. 		src++;

  305. 	}

  306. }	

  307. 

  308. void sample_move_d32u24_sS (char *dst, jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)

  309. {

  310. #if defined (__SSE2__) && !defined (__sun__)

  311. 	__m128 int_max = _mm_set1_ps(SAMPLE_24BIT_MAX_F);

  312. 	__m128 int_min = _mm_sub_ps(_mm_setzero_ps(), int_max);

  313. 	__m128 factor = int_max;

  314. 

  315. 	unsigned long unrolled = nsamples / 4;

  316. 	nsamples = nsamples & 3;

  317. 

  318. 	while (unrolled--) {

  319. 		__m128 in = _mm_load_ps(src);

  320. 		__m128 scaled = _mm_mul_ps(in, factor);

  321. 		__m128 clipped = clip(scaled, int_min, int_max);

  322. 

  323. 		__m128i y = _mm_cvttps_epi32(clipped);

  324. 

  325. #ifdef __SSE4_1__

  326. 		*(int32_t*)dst              = _mm_extract_epi32(y, 0);

  327. 		*(int32_t*)(dst+dst_skip)   = _mm_extract_epi32(y, 1);

  328. 		*(int32_t*)(dst+2*dst_skip) = _mm_extract_epi32(y, 2);

  329. 		*(int32_t*)(dst+3*dst_skip) = _mm_extract_epi32(y, 3);

  330. #else

  331. 		__m128i shuffled1 = _mm_shuffle_epi32(y, _MM_SHUFFLE(0, 3, 2, 1));

  332. 		__m128i shuffled2 = _mm_shuffle_epi32(y, _MM_SHUFFLE(1, 0, 3, 2));

  333. 		__m128i shuffled3 = _mm_shuffle_epi32(y, _MM_SHUFFLE(2, 1, 0, 3));

  334. 

  335. 		_mm_store_ss((float*)dst, (__m128)y);

  336. 

  337. 		_mm_store_ss((float*)(dst+dst_skip), (__m128)shuffled1);

  338. 		_mm_store_ss((float*)(dst+2*dst_skip), (__m128)shuffled2);

  339. 		_mm_store_ss((float*)(dst+3*dst_skip), (__m128)shuffled3);

  340. #endif

  341. 		dst += 4*dst_skip;

  342. 

  343. 		src+= 4;

  344. 	}

  345. 

  346. 	while (nsamples--) {

  347. 		__m128 in = _mm_load_ss(src);

  348. 		__m128 scaled = _mm_mul_ss(in, factor);

  349. 		__m128 clipped = _mm_min_ss(int_max, _mm_max_ss(scaled, int_min));

  350. 

  351. 		int y = _mm_cvttss_si32(clipped);

  352. 		*((int *) dst) = y;

  353. 

  354. 		dst += dst_skip;

  355. 		src++;

  356. 	}

  357. 

  358. #elif defined (__ARM_NEON__) || defined (__ARM_NEON)

  359. 	unsigned long unrolled = nsamples / 4;

  360. 	nsamples = nsamples & 3;

  361. 

  362. 	while (unrolled--) {

  363. 		float32x4_t samples = vld1q_f32(src);

  364. 		int32x4_t converted = float_24_neon(samples);

  365. 

  366. 		switch(dst_skip) {

  367. 			case 4:

  368. 				vst1q_s32((int32_t*)dst, converted);

  369. 				break;

  370. 			default:

  371. 				vst1q_lane_s32((int32_t*)(dst),            converted, 0);

  372. 				vst1q_lane_s32((int32_t*)(dst+dst_skip),   converted, 1);

  373. 				vst1q_lane_s32((int32_t*)(dst+2*dst_skip), converted, 2);

  374. 				vst1q_lane_s32((int32_t*)(dst+3*dst_skip), converted, 3);

  375. 				break;

  376. 		}

  377. 		dst += 4*dst_skip;

  378. 

  379. 		src+= 4;

  380. 	}

  381. #endif

  382. 

  383. #if !defined (__SSE2__)

  384. 	while (nsamples--) {

  385. 		float_24u32 (*src, *((int32_t*) dst));

  386. 		dst += dst_skip;

  387. 		src++;

  388. 	}

  389. #endif

  390. }	

  391. 

  392. void sample_move_dS_s32u24s (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip)

  393. {

  394. 	const jack_default_audio_sample_t scaling = 1.0 / (SAMPLE_24BIT_SCALING << 8);

  395. 

  396. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  397. 	float32x4_t factor = vdupq_n_f32(scaling);

  398. 	unsigned long unrolled = nsamples / 4;

  399. 	while (unrolled--) {

  400. 		int32x4_t src128;

  401. 		switch(src_skip)

  402. 		{

  403. 			case 4:

  404. 				src128 = vld1q_s32((int32_t*)src);

  405. 				break;

  406. 			case 8:

  407. 				src128 = vld2q_s32((int32_t*)src).val[0];

  408. 				break;

  409. 			default:

  410. 				src128 = vld1q_lane_s32((int32_t*)src,              src128, 0);

  411. 				src128 = vld1q_lane_s32((int32_t*)(src+src_skip),   src128, 1);

  412. 				src128 = vld1q_lane_s32((int32_t*)(src+2*src_skip), src128, 2);

  413. 				src128 = vld1q_lane_s32((int32_t*)(src+3*src_skip), src128, 3);

  414. 				break;

  415. 		}

  416. 		src128 = vreinterpretq_s32_u8(vrev32q_u8(vreinterpretq_u8_s32(src128)));

  417. 		/* sign extension - left shift will be reverted by scaling */

  418. 		int32x4_t shifted = vshlq_n_s32(src128, 8);

  419. 		float32x4_t as_float = vcvtq_f32_s32(shifted);

  420. 		float32x4_t divided = vmulq_f32(as_float, factor);

  421. 		vst1q_f32(dst, divided);

  422. 

  423. 		src += 4*src_skip;

  424. 		dst += 4;

  425. 	}

  426. 	nsamples = nsamples & 3;

  427. #endif

  428. 

  429. 	/* ALERT: signed sign-extension portability !!! */

  430. 

  431. 	while (nsamples--) {

  432. 		int x;

  433. #if __BYTE_ORDER == __LITTLE_ENDIAN

  434. 		x = (unsigned char)(src[0]);

  435. 		x <<= 8;

  436. 		x |= (unsigned char)(src[1]);

  437. 		x <<= 8;

  438. 		x |= (unsigned char)(src[2]);

  439. 		x <<= 8;

  440. 		x |= (unsigned char)(src[3]);

  441. #elif __BYTE_ORDER == __BIG_ENDIAN

  442. 		x = (unsigned char)(src[3]);

  443. 		x <<= 8;

  444. 		x |= (unsigned char)(src[2]);

  445. 		x <<= 8;

  446. 		x |= (unsigned char)(src[1]);

  447. 		x <<= 8;

  448. 		x |= (unsigned char)(src[0]);

  449. #endif

  450. 		/* sign extension - left shift will be reverted by scaling */

  451. 		*dst = (x << 8) * scaling;

  452. 		dst++;

  453. 		src += src_skip;

  454. 	}

  455. }	

  456. 

  457. void sample_move_dS_s32u24 (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip)

  458. {

  459. 	const jack_default_audio_sample_t scaling = 1.0 / (SAMPLE_24BIT_SCALING << 8);

  460. 

  461. #if defined (__SSE2__) && !defined (__sun__)

  462. 	unsigned long unrolled = nsamples / 4;

  463. 	__m128 factor = _mm_set1_ps(scaling);

  464. 	while (unrolled--)

  465. 	{

  466. 		int i1 = *((int *) src);

  467. 		src+= src_skip;

  468. 		int i2 = *((int *) src);

  469. 		src+= src_skip;

  470. 		int i3 = *((int *) src);

  471. 		src+= src_skip;

  472. 		int i4 = *((int *) src);

  473. 		src+= src_skip;

  474. 

  475. 		__m128i src = _mm_set_epi32(i4, i3, i2, i1);

  476. 		/* sign extension - left shift will be reverted by scaling */

  477. 		__m128i shifted = _mm_slli_epi32(src, 8);

  478. 

  479. 		__m128 as_float = _mm_cvtepi32_ps(shifted);

  480. 		__m128 divided = _mm_mul_ps(as_float, factor);

  481. 

  482. 		_mm_storeu_ps(dst, divided);

  483. 

  484. 		dst += 4;

  485. 	}

  486. 	nsamples = nsamples & 3;

  487. #elif defined (__ARM_NEON__) || defined (__ARM_NEON)

  488. 	unsigned long unrolled = nsamples / 4;

  489. 	float32x4_t factor = vdupq_n_f32(scaling);

  490. 	while (unrolled--) {

  491. 		int32x4_t src128;

  492. 		switch(src_skip) {

  493. 			case 4:

  494. 				src128 = vld1q_s32((int32_t*)src);

  495. 				break;

  496. 			case 8:

  497. 				src128 = vld2q_s32((int32_t*)src).val[0];

  498. 				break;

  499. 			default:

  500. 				src128 = vld1q_lane_s32((int32_t*)src,              src128, 0);

  501. 				src128 = vld1q_lane_s32((int32_t*)(src+src_skip),   src128, 1);

  502. 				src128 = vld1q_lane_s32((int32_t*)(src+2*src_skip), src128, 2);

  503. 				src128 = vld1q_lane_s32((int32_t*)(src+3*src_skip), src128, 3);

  504. 				break;

  505. 		}

  506. 		/* sign extension  - left shift will be reverted by scaling */

  507. 		int32x4_t shifted = vshlq_n_s32(src128, 8);

  508. 		float32x4_t as_float = vcvtq_f32_s32(shifted);

  509. 		float32x4_t divided = vmulq_f32(as_float, factor);

  510. 		vst1q_f32(dst, divided);

  511. 

  512. 		src += 4*src_skip;

  513. 		dst += 4;

  514. 	}

  515. 	nsamples = nsamples & 3;

  516. #endif

  517. 

  518. 	/* ALERT: signed sign-extension portability !!! */

  519. 

  520. 	while (nsamples--) {

  521. 		/* sign extension  - left shift will be reverted by scaling */

  522. 		*dst = (*((int *) src) << 8) * scaling;

  523. 		dst++;

  524. 		src += src_skip;

  525. 	}

  526. }	

  527. 

  528. void sample_move_d24_sSs (char *dst, jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)

  529. {

  530. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  531. 	unsigned long unrolled = nsamples / 4;

  532. 	while (unrolled--) {

  533. 		int i;

  534. 		int32_t z[4];

  535. 		float32x4_t samples = vld1q_f32(src);

  536. 		int32x4_t converted = float_24_neon(samples);

  537. 		converted = vreinterpretq_s32_u8(vrev32q_u8(vreinterpretq_u8_s32(converted)));

  538. 		vst1q_s32(z, converted);

  539. 

  540. 		for (i = 0; i != 4; ++i) {

  541. 			memcpy (dst, ((char*)(z+i))+1, 3);

  542. 			dst += dst_skip;

  543. 		}

  544. 		src += 4;

  545. 	}

  546. 	nsamples = nsamples & 3;

  547. #endif

  548. 

  549. 	int32_t z;

  550. 

  551. 	while (nsamples--) {

  552. 		float_24 (*src, z);

  553. #if __BYTE_ORDER == __LITTLE_ENDIAN

  554. 		dst[0]=(char)(z>>16);

  555. 		dst[1]=(char)(z>>8);

  556. 		dst[2]=(char)(z);

  557. #elif __BYTE_ORDER == __BIG_ENDIAN

  558. 		dst[0]=(char)(z);

  559. 		dst[1]=(char)(z>>8);

  560. 		dst[2]=(char)(z>>16);

  561. #endif

  562. 		dst += dst_skip;

  563. 		src++;

  564. 	}

  565. }	

  566. 

  567. void sample_move_d24_sS (char *dst, jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)

  568. {

  569. #if defined (__SSE2__) && !defined (__sun__)

  570. 	_MM_SET_ROUNDING_MODE(_MM_ROUND_NEAREST);

  571. 	while (nsamples >= 4) {

  572. 		int i;

  573. 		int32_t z[4];

  574. 		__m128 samples = _mm_loadu_ps(src);

  575. 		__m128i converted = float_24_sse(samples);

  576. 

  577. #ifdef __SSE4_1__

  578. 		z[0] = _mm_extract_epi32(converted, 0);

  579. 		z[1] = _mm_extract_epi32(converted, 1);

  580. 		z[2] = _mm_extract_epi32(converted, 2);

  581. 		z[3] = _mm_extract_epi32(converted, 3);

  582. #else

  583. 		__m128i shuffled1 = _mm_shuffle_epi32(converted, _MM_SHUFFLE(0, 3, 2, 1));

  584. 		__m128i shuffled2 = _mm_shuffle_epi32(converted, _MM_SHUFFLE(1, 0, 3, 2));

  585. 		__m128i shuffled3 = _mm_shuffle_epi32(converted, _MM_SHUFFLE(2, 1, 0, 3));

  586. 

  587. 		_mm_store_ss((float*)z, (__m128)converted);

  588. 		_mm_store_ss((float*)z+1, (__m128)shuffled1);

  589. 		_mm_store_ss((float*)z+2, (__m128)shuffled2);

  590. 		_mm_store_ss((float*)z+3, (__m128)shuffled3);

  591. #endif

  592. 

  593. 		for (i = 0; i != 4; ++i) {

  594. 			memcpy (dst, z+i, 3);

  595. 			dst += dst_skip;

  596. 		}

  597. 

  598. 		nsamples -= 4;

  599. 		src += 4;

  600. 	}

  601. #elif defined (__ARM_NEON__) || defined (__ARM_NEON)

  602. 	unsigned long unrolled = nsamples / 4;

  603. 	while (unrolled--) {

  604. 		int i;

  605. 		int32_t z[4];

  606. 		float32x4_t samples = vld1q_f32(src);

  607. 		int32x4_t converted = float_24_neon(samples);

  608. 		vst1q_s32(z, converted);

  609. 

  610. 		for (i = 0; i != 4; ++i) {

  611. 			memcpy (dst, z+i, 3);

  612. 			dst += dst_skip;

  613. 		}

  614. 		src += 4;

  615. 	}

  616. 	nsamples = nsamples & 3;

  617. #endif

  618. 

  619.     int32_t z;

  620. 

  621. 	while (nsamples--) {

  622. 		float_24 (*src, z);

  623. #if __BYTE_ORDER == __LITTLE_ENDIAN

  624. 		memcpy (dst, &z, 3);

  625. #elif __BYTE_ORDER == __BIG_ENDIAN

  626. 		memcpy (dst, (char *)&z + 1, 3);

  627. #endif

  628. 		dst += dst_skip;

  629. 		src++;

  630. 	}

  631. }

  632. 

  633. void sample_move_dS_s24s (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip)

  634. {

  635. 	const jack_default_audio_sample_t scaling = 1.0/SAMPLE_24BIT_SCALING;

  636. 

  637. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  638. 	// we shift 8 to the right by dividing by 256.0 -> no sign extra handling

  639. 	const float32x4_t vscaling = vdupq_n_f32(scaling/256.0);

  640. 	int32_t x[4];

  641. 	memset(x, 0, sizeof(x));

  642. 	unsigned long unrolled = nsamples / 4;

  643. 	while (unrolled--) {

  644. #if __BYTE_ORDER == __BIG_ENDIAN	 /* ARM big endian?? */

  645. 		// right aligned / inverse sequence below -> *256

  646. 		memcpy(((char*)&x[0])+1, src, 3);

  647. 		memcpy(((char*)&x[1])+1, src+src_skip, 3);

  648. 		memcpy(((char*)&x[2])+1, src+2*src_skip, 3);

  649. 		memcpy(((char*)&x[3])+1, src+3*src_skip, 3);

  650. #else

  651. 		memcpy(&x[0], src, 3);

  652. 		memcpy(&x[1], src+src_skip, 3);

  653. 		memcpy(&x[2], src+2*src_skip, 3);

  654. 		memcpy(&x[3], src+3*src_skip, 3);

  655. #endif

  656. 		src += 4 * src_skip;

  657. 

  658. 		int32x4_t source = vld1q_s32(x);

  659. 		source = vreinterpretq_s32_u8(vrev32q_u8(vreinterpretq_u8_s32(source)));

  660. 		float32x4_t converted = vcvtq_f32_s32(source);

  661. 		float32x4_t scaled = vmulq_f32(converted, vscaling);

  662. 		vst1q_f32(dst, scaled);

  663. 		dst += 4;

  664. 	}

  665. 	nsamples = nsamples & 3;

  666. #endif

  667. 

  668. 	/* ALERT: signed sign-extension portability !!! */

  669. 

  670. 	while (nsamples--) {

  671. 		int x;

  672. #if __BYTE_ORDER == __LITTLE_ENDIAN

  673. 		x = (unsigned char)(src[0]);

  674. 		x <<= 8;

  675. 		x |= (unsigned char)(src[1]);

  676. 		x <<= 8;

  677. 		x |= (unsigned char)(src[2]);

  678. 		/* correct sign bit and the rest of the top byte */

  679. 		if (src[0] & 0x80) {

  680. 			x |= 0xff << 24;

  681. 		}

  682. #elif __BYTE_ORDER == __BIG_ENDIAN

  683. 		x = (unsigned char)(src[2]);

  684. 		x <<= 8;

  685. 		x |= (unsigned char)(src[1]);

  686. 		x <<= 8;

  687. 		x |= (unsigned char)(src[0]);

  688. 		/* correct sign bit and the rest of the top byte */

  689. 		if (src[2] & 0x80) {

  690. 			x |= 0xff << 24;

  691. 		}

  692. #endif

  693. 		*dst = x * scaling;

  694. 		dst++;

  695. 		src += src_skip;

  696. 	}

  697. }

  698. 

  699. void sample_move_dS_s24 (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip)

  700. {

  701. 	const jack_default_audio_sample_t scaling = 1.f/SAMPLE_24BIT_SCALING;

  702. 

  703. #if defined (__SSE2__) && !defined (__sun__)

  704. 	const __m128 scaling_block = _mm_set_ps1(scaling);

  705. 	while (nsamples >= 4) {

  706. 		int x0, x1, x2, x3;

  707. 

  708. 		memcpy((char*)&x0 + 1, src, 3);

  709. 		memcpy((char*)&x1 + 1, src+src_skip, 3);

  710. 		memcpy((char*)&x2 + 1, src+2*src_skip, 3);

  711. 		memcpy((char*)&x3 + 1, src+3*src_skip, 3);

  712. 		src += 4 * src_skip;

  713. 

  714. 		const __m128i block_i = _mm_set_epi32(x3, x2, x1, x0);

  715. 		const __m128i shifted = _mm_srai_epi32(block_i, 8);

  716. 		const __m128 converted = _mm_cvtepi32_ps (shifted);

  717. 		const __m128 scaled = _mm_mul_ps(converted, scaling_block);

  718. 		_mm_storeu_ps(dst, scaled);

  719. 		dst += 4;

  720. 		nsamples -= 4;

  721. 	}

  722. #elif defined (__ARM_NEON__) || defined (__ARM_NEON)

  723. 	// we shift 8 to the right by dividing by 256.0 -> no sign extra handling

  724. 	const float32x4_t vscaling = vdupq_n_f32(scaling/256.0);

  725. 	int32_t x[4];

  726. 	memset(x, 0, sizeof(x));

  727. 	unsigned long unrolled = nsamples / 4;

  728. 	while (unrolled--) {

  729. #if __BYTE_ORDER == __BIG_ENDIAN	/* ARM big endian?? */

  730. 		// left aligned -> *256

  731. 		memcpy(&x[0], src, 3);

  732. 		memcpy(&x[1], src+src_skip, 3);

  733. 		memcpy(&x[2], src+2*src_skip, 3);

  734. 		memcpy(&x[3], src+3*src_skip, 3);

  735. #else

  736. 		memcpy(((char*)&x[0])+1, src, 3);

  737. 		memcpy(((char*)&x[1])+1, src+src_skip, 3);

  738. 		memcpy(((char*)&x[2])+1, src+2*src_skip, 3);

  739. 		memcpy(((char*)&x[3])+1, src+3*src_skip, 3);

  740. #endif

  741. 		src += 4 * src_skip;

  742. 

  743. 		int32x4_t source = vld1q_s32(x);

  744. 		float32x4_t converted = vcvtq_f32_s32(source);

  745. 		float32x4_t scaled = vmulq_f32(converted, vscaling);

  746. 		vst1q_f32(dst, scaled);

  747. 		dst += 4;

  748. 	}

  749. 	nsamples = nsamples & 3;

  750. #endif

  751. 

  752. 	while (nsamples--) {

  753. 		int x;

  754. #if __BYTE_ORDER == __LITTLE_ENDIAN

  755. 		memcpy((char*)&x + 1, src, 3);

  756. #elif __BYTE_ORDER == __BIG_ENDIAN

  757. 		memcpy(&x, src, 3);

  758. #endif

  759. 		x >>= 8;

  760. 		*dst = x * scaling;

  761. 		dst++;

  762. 		src += src_skip;

  763. 	}

  764. }

  765. 

  766. 

  767. void sample_move_d16_sSs (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  768. {

  769. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  770. 	unsigned long unrolled = nsamples / 4;

  771. 	nsamples = nsamples & 3;

  772. 

  773. 	while (unrolled--) {

  774. 		float32x4_t samples = vld1q_f32(src);

  775. 		int16x4_t converted = float_16_neon(samples);

  776. 		converted = vreinterpret_s16_u8(vrev16_u8(vreinterpret_u8_s16(converted)));

  777. 

  778. 		switch(dst_skip) {

  779. 			case 2:

  780. 				vst1_s16((int16_t*)dst, converted);

  781. 				break;

  782. 			default:

  783. 				vst1_lane_s16((int16_t*)(dst),            converted, 0);

  784. 				vst1_lane_s16((int16_t*)(dst+dst_skip),   converted, 1);

  785. 				vst1_lane_s16((int16_t*)(dst+2*dst_skip), converted, 2);

  786. 				vst1_lane_s16((int16_t*)(dst+3*dst_skip), converted, 3);

  787. 				break;

  788. 		}

  789. 		dst += 4*dst_skip;

  790. 		src+= 4;

  791. 	}

  792. #endif

  793. 	int16_t tmp;

  794. 

  795. 	while (nsamples--) {

  796. 		// float_16 (*src, tmp);

  797. 

  798. 		if (*src <= NORMALIZED_FLOAT_MIN) {

  799. 			tmp = SAMPLE_16BIT_MIN;

  800. 		} else if (*src >= NORMALIZED_FLOAT_MAX) {

  801. 			tmp = SAMPLE_16BIT_MAX;

  802. 		} else {

  803. 			tmp = (int16_t) f_round (*src * SAMPLE_16BIT_SCALING);

  804. 		}

  805. 

  806. #if __BYTE_ORDER == __LITTLE_ENDIAN

  807. 		dst[0]=(char)(tmp>>8);

  808. 		dst[1]=(char)(tmp);

  809. #elif __BYTE_ORDER == __BIG_ENDIAN

  810. 		dst[0]=(char)(tmp);

  811. 		dst[1]=(char)(tmp>>8);

  812. #endif

  813. 		dst += dst_skip;

  814. 		src++;

  815. 	}

  816. }

  817. 

  818. void sample_move_d16_sS (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  819. {

  820. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  821. 	unsigned long unrolled = nsamples / 4;

  822. 	nsamples = nsamples & 3;

  823. 

  824. 	while (unrolled--) {

  825. 		float32x4_t samples = vld1q_f32(src);

  826. 		int16x4_t converted = float_16_neon(samples);

  827. 

  828. 		switch(dst_skip) {

  829. 			case 2:

  830. 				vst1_s16((int16_t*)dst, converted);

  831. 				break;

  832. 			default:

  833. 				vst1_lane_s16((int16_t*)(dst),            converted, 0);

  834. 				vst1_lane_s16((int16_t*)(dst+dst_skip),   converted, 1);

  835. 				vst1_lane_s16((int16_t*)(dst+2*dst_skip), converted, 2);

  836. 				vst1_lane_s16((int16_t*)(dst+3*dst_skip), converted, 3);

  837. 				break;

  838. 		}

  839. 		dst += 4*dst_skip;

  840. 		src+= 4;

  841. 	}

  842. #endif

  843. 	while (nsamples--) {

  844. 		float_16 (*src, *((int16_t*) dst));

  845. 		dst += dst_skip;

  846. 		src++;

  847. 	}

  848. }

  849. 

  850. void sample_move_dither_rect_d16_sSs (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  851. {

  852. 	jack_default_audio_sample_t val;

  853. 	int16_t      tmp;

  854. 

  855. 	while (nsamples--) {

  856. 		val = (*src * SAMPLE_16BIT_SCALING) + fast_rand() / (float) UINT_MAX - 0.5f;

  857. 		float_16_scaled (val, tmp);

  858. #if __BYTE_ORDER == __LITTLE_ENDIAN

  859. 		dst[0]=(char)(tmp>>8);

  860. 		dst[1]=(char)(tmp);

  861. #elif __BYTE_ORDER == __BIG_ENDIAN

  862. 		dst[0]=(char)(tmp);

  863. 		dst[1]=(char)(tmp>>8);

  864. #endif

  865. 		dst += dst_skip;

  866. 		src++;

  867. 	}

  868. }

  869. 

  870. void sample_move_dither_rect_d16_sS (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  871. {

  872. 	jack_default_audio_sample_t val;

  873. 

  874. 	while (nsamples--) {

  875. 		val = (*src * SAMPLE_16BIT_SCALING) + fast_rand() / (float)UINT_MAX - 0.5f;

  876. 		float_16_scaled (val, *((int16_t*) dst));

  877. 		dst += dst_skip;

  878. 		src++;

  879. 	}

  880. }

  881. 

  882. void sample_move_dither_tri_d16_sSs (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  883. {

  884. 	jack_default_audio_sample_t val;

  885. 	int16_t      tmp;

  886. 

  887. 	while (nsamples--) {

  888. 		val = (*src * SAMPLE_16BIT_SCALING) + ((float)fast_rand() + (float)fast_rand()) / (float)UINT_MAX - 1.0f;

  889. 		float_16_scaled (val, tmp);

  890. 

  891. #if __BYTE_ORDER == __LITTLE_ENDIAN

  892. 		dst[0]=(char)(tmp>>8);

  893. 		dst[1]=(char)(tmp);

  894. #elif __BYTE_ORDER == __BIG_ENDIAN

  895. 		dst[0]=(char)(tmp);

  896. 		dst[1]=(char)(tmp>>8);

  897. #endif

  898. 		dst += dst_skip;

  899. 		src++;

  900. 	}

  901. }

  902. 

  903. void sample_move_dither_tri_d16_sS (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  904. {

  905. 	jack_default_audio_sample_t val;

  906. 

  907. 	while (nsamples--) {

  908. 		val = (*src * SAMPLE_16BIT_SCALING) + ((float)fast_rand() + (float)fast_rand()) / (float)UINT_MAX - 1.0f;

  909. 		float_16_scaled (val, *((int16_t*) dst));

  910. 		dst += dst_skip;

  911. 		src++;

  912. 	}

  913. }

  914. 

  915. void sample_move_dither_shaped_d16_sSs (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  916. {

  917. 	jack_default_audio_sample_t     x;

  918. 	jack_default_audio_sample_t     xe; /* the innput sample - filtered error */

  919. 	jack_default_audio_sample_t     xp; /* x' */

  920. 	float        r;

  921. 	float        rm1 = state->rm1;

  922. 	unsigned int idx = state->idx;

  923. 	int16_t      tmp;

  924. 

  925. 	while (nsamples--) {

  926. 		x = *src * SAMPLE_16BIT_SCALING;

  927. 		r = ((float)fast_rand() + (float)fast_rand())  / (float)UINT_MAX - 1.0f;

  928. 		/* Filter the error with Lipshitz's minimally audible FIR:

  929. 		   [2.033 -2.165 1.959 -1.590 0.6149] */

  930. 		xe = x

  931. 		     - state->e[idx] * 2.033f

  932. 		     + state->e[(idx - 1) & DITHER_BUF_MASK] * 2.165f

  933. 		     - state->e[(idx - 2) & DITHER_BUF_MASK] * 1.959f

  934. 		     + state->e[(idx - 3) & DITHER_BUF_MASK] * 1.590f

  935. 		     - state->e[(idx - 4) & DITHER_BUF_MASK] * 0.6149f;

  936. 		xp = xe + r - rm1;

  937. 		rm1 = r;

  938. 

  939. 		float_16_scaled (xp, tmp);

  940. 

  941. 		/* Intrinsic z^-1 delay */

  942. 		idx = (idx + 1) & DITHER_BUF_MASK;

  943. 		state->e[idx] = xp - xe;

  944. 

  945. #if __BYTE_ORDER == __LITTLE_ENDIAN

  946. 		dst[0]=(char)(tmp>>8);

  947. 		dst[1]=(char)(tmp);

  948. #elif __BYTE_ORDER == __BIG_ENDIAN

  949. 		dst[0]=(char)(tmp);

  950. 		dst[1]=(char)(tmp>>8);

  951. #endif

  952. 		dst += dst_skip;

  953. 		src++;

  954. 	}

  955. 	state->rm1 = rm1;

  956. 	state->idx = idx;

  957. }

  958. 

  959. void sample_move_dither_shaped_d16_sS (char *dst,  jack_default_audio_sample_t *src, unsigned long nsamples, unsigned long dst_skip, dither_state_t *state)	

  960. {

  961. 	jack_default_audio_sample_t     x;

  962. 	jack_default_audio_sample_t     xe; /* the innput sample - filtered error */

  963. 	jack_default_audio_sample_t     xp; /* x' */

  964. 	float        r;

  965. 	float        rm1 = state->rm1;

  966. 	unsigned int idx = state->idx;

  967. 

  968. 	while (nsamples--) {

  969. 		x = *src * SAMPLE_16BIT_SCALING;

  970. 		r = ((float)fast_rand() + (float)fast_rand()) / (float)UINT_MAX - 1.0f;

  971. 		/* Filter the error with Lipshitz's minimally audible FIR:

  972. 		   [2.033 -2.165 1.959 -1.590 0.6149] */

  973. 		xe = x

  974. 		     - state->e[idx] * 2.033f

  975. 		     + state->e[(idx - 1) & DITHER_BUF_MASK] * 2.165f

  976. 		     - state->e[(idx - 2) & DITHER_BUF_MASK] * 1.959f

  977. 		     + state->e[(idx - 3) & DITHER_BUF_MASK] * 1.590f

  978. 		     - state->e[(idx - 4) & DITHER_BUF_MASK] * 0.6149f;

  979. 		xp = xe + r - rm1;

  980. 		rm1 = r;

  981. 

  982. 		float_16_scaled (xp, *((int16_t*) dst));

  983. 

  984. 		/* Intrinsic z^-1 delay */

  985. 		idx = (idx + 1) & DITHER_BUF_MASK;

  986. 		state->e[idx] = *((int16_t*) dst) - xe;

  987. 

  988. 		dst += dst_skip;

  989. 		src++;

  990. 	}

  991. 	state->rm1 = rm1;

  992. 	state->idx = idx;

  993. }

  994. 

  995. void sample_move_dS_s16s (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip) 	

  996. {

  997. 	short z;

  998. 	const jack_default_audio_sample_t scaling = 1.0/SAMPLE_16BIT_SCALING;

  999. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  1000. 	const float32x4_t vscaling = vdupq_n_f32(scaling);

  1001. 	unsigned long unrolled = nsamples / 4;

  1002. 	while (unrolled--) {

  1003. 		int16x4_t source16x4;

  1004. 		switch(src_skip) {

  1005. 			case 2:

  1006. 				source16x4 = vld1_s16((int16_t*)src);

  1007. 				break;

  1008. 			case 4:

  1009. 				source16x4 = vld2_s16((int16_t*)src).val[0];

  1010. 				break;

  1011. 			default:

  1012. 				source16x4 = vld1_lane_s16((int16_t*)src,              source16x4, 0);

  1013. 				source16x4 = vld1_lane_s16((int16_t*)(src+src_skip),   source16x4, 1);

  1014. 				source16x4 = vld1_lane_s16((int16_t*)(src+2*src_skip), source16x4, 2);

  1015. 				source16x4 = vld1_lane_s16((int16_t*)(src+3*src_skip), source16x4, 3);

  1016. 				break;

  1017. 		}

  1018. 		source16x4 = vreinterpret_s16_u8(vrev16_u8(vreinterpret_u8_s16(source16x4)));

  1019. 		int32x4_t source32x4 = vmovl_s16(source16x4);

  1020. 		src += 4 * src_skip;

  1021. 

  1022. 		float32x4_t converted = vcvtq_f32_s32(source32x4);

  1023. 		float32x4_t scaled = vmulq_f32(converted, vscaling);

  1024. 		vst1q_f32(dst, scaled);

  1025. 		dst += 4;

  1026. 	}

  1027. 	nsamples = nsamples & 3;

  1028. #endif

  1029. 

  1030. 	/* ALERT: signed sign-extension portability !!! */

  1031. 	while (nsamples--) {

  1032. #if __BYTE_ORDER == __LITTLE_ENDIAN

  1033. 		z = (unsigned char)(src[0]);

  1034. 		z <<= 8;

  1035. 		z |= (unsigned char)(src[1]);

  1036. #elif __BYTE_ORDER == __BIG_ENDIAN

  1037. 		z = (unsigned char)(src[1]);

  1038. 		z <<= 8;

  1039. 		z |= (unsigned char)(src[0]);

  1040. #endif

  1041. 		*dst = z * scaling;

  1042. 		dst++;

  1043. 		src += src_skip;

  1044. 	}

  1045. }	

  1046. 

  1047. void sample_move_dS_s16 (jack_default_audio_sample_t *dst, char *src, unsigned long nsamples, unsigned long src_skip) 

  1048. {

  1049. 	/* ALERT: signed sign-extension portability !!! */

  1050. 	const jack_default_audio_sample_t scaling = 1.0/SAMPLE_16BIT_SCALING;

  1051. #if defined (__ARM_NEON__) || defined (__ARM_NEON)

  1052. 	const float32x4_t vscaling = vdupq_n_f32(scaling);

  1053. 	unsigned long unrolled = nsamples / 4;

  1054. 	while (unrolled--) {

  1055. 		int16x4_t source16x4;

  1056. 		switch(src_skip) {

  1057. 			case 2:

  1058. 				source16x4 = vld1_s16((int16_t*)src);

  1059. 				break;

  1060. 			case 4:

  1061. 				source16x4 = vld2_s16((int16_t*)src).val[0];

  1062. 				break;

  1063. 			default:

  1064. 				source16x4 = vld1_lane_s16((int16_t*)src,              source16x4, 0);

  1065. 				source16x4 = vld1_lane_s16((int16_t*)(src+src_skip),   source16x4, 1);

  1066. 				source16x4 = vld1_lane_s16((int16_t*)(src+2*src_skip), source16x4, 2);

  1067. 				source16x4 = vld1_lane_s16((int16_t*)(src+3*src_skip), source16x4, 3);

  1068. 				break;

  1069. 		}

  1070. 		int32x4_t source32x4 = vmovl_s16(source16x4);

  1071. 		src += 4 * src_skip;

  1072. 

  1073. 		float32x4_t converted = vcvtq_f32_s32(source32x4);

  1074. 		float32x4_t scaled = vmulq_f32(converted, vscaling);

  1075. 		vst1q_f32(dst, scaled);

  1076. 		dst += 4;

  1077. 	}

  1078. 	nsamples = nsamples & 3;

  1079. #endif

  1080. 

  1081. 	while (nsamples--) {

  1082. 		*dst = (*((short *) src)) * scaling;

  1083. 		dst++;

  1084. 		src += src_skip;

  1085. 	}

  1086. }	

  1087. 

  1088. void memset_interleave (char *dst, char val, unsigned long bytes, 

  1089. 			unsigned long unit_bytes, 

  1090. 			unsigned long skip_bytes) 

  1091. {

  1092. 	switch (unit_bytes) {

  1093. 	case 1:

  1094. 		while (bytes--) {

  1095. 			*dst = val;

  1096. 			dst += skip_bytes;

  1097. 		}

  1098. 		break;

  1099. 	case 2:

  1100. 		while (bytes) {

  1101. 			*((short *) dst) = (short) val;

  1102. 			dst += skip_bytes;

  1103. 			bytes -= 2;

  1104. 		}

  1105. 		break;

  1106. 	case 4:		    

  1107. 		while (bytes) {

  1108. 			*((int *) dst) = (int) val;

  1109. 			dst += skip_bytes;

  1110. 			bytes -= 4;

  1111. 		}

  1112. 		break;

  1113. 	default:

  1114. 		while (bytes) {

  1115. 			memset(dst, val, unit_bytes);

  1116. 			dst += skip_bytes;

  1117. 			bytes -= unit_bytes;

  1118. 		}

  1119. 		break;

  1120. 	}

  1121. }

  1122. 

  1123. /* COPY FUNCTIONS: used to move data from an input channel to an

  1124.    output channel. Note that we assume that the skip distance

  1125.    is the same for both channels. This is completely fine

  1126.    unless the input and output were on different audio interfaces that

  1127.    were interleaved differently. We don't try to handle that.

  1128. */

  1129. 

  1130. void 

  1131. memcpy_fake (char *dst, char *src, unsigned long src_bytes, unsigned long foo, unsigned long bar)

  1132. {

  1133. 	memcpy (dst, src, src_bytes);

  1134. }

  1135. 

  1136. void 

  1137. memcpy_interleave_d16_s16 (char *dst, char *src, unsigned long src_bytes,

  1138. 			   unsigned long dst_skip_bytes, unsigned long src_skip_bytes)

  1139. {

  1140. 	while (src_bytes) {

  1141. 		*((short *) dst) = *((short *) src);

  1142. 		dst += dst_skip_bytes;

  1143. 		src += src_skip_bytes;

  1144. 		src_bytes -= 2;

  1145. 	}

  1146. }

  1147. 

  1148. void 

  1149. memcpy_interleave_d24_s24 (char *dst, char *src, unsigned long src_bytes,

  1150. 			   unsigned long dst_skip_bytes, unsigned long src_skip_bytes)

  1151. {

  1152. 	while (src_bytes) {

  1153. 		memcpy(dst, src, 3);

  1154. 		dst += dst_skip_bytes;

  1155. 		src += src_skip_bytes;

  1156. 		src_bytes -= 3;

  1157. 	}

  1158. }

  1159. 

  1160. void 

  1161. memcpy_interleave_d32_s32 (char *dst, char *src, unsigned long src_bytes,

  1162. 			   unsigned long dst_skip_bytes, unsigned long src_skip_bytes)

  1163. {

  1164. 	while (src_bytes) {

  1165. 		*((int *) dst) = *((int *) src);

  1166. 		dst += dst_skip_bytes;

  1167. 		src += src_skip_bytes;

  1168. 		src_bytes -= 4;

  1169. 	}

  1170. }