musix 0.3.3

Music player library for esoteric audio formats (music from C64,Amiga etc)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110

/******************************************************************************\
* Authors:  Iconoclast                                                         *
* Release:  2013.11.26                                                         *
* License:  CC0 Public Domain Dedication                                       *
*                                                                              *
* To the extent possible under law, the author(s) have dedicated all copyright *
* and related and neighboring rights to this software to the public domain     *
* worldwide. This software is distributed without any warranty.                *
*                                                                              *
* You should have received a copy of the CC0 Public Domain Dedication along    *
* with this software.                                                          *
* If not, see <http://creativecommons.org/publicdomain/zero/1.0/>.             *
\******************************************************************************/
#include "vu.h"

INLINE static void do_cr(usf_state_t * state, short* VD, short* VS, short* VT)
{
    ALIGNED short ge[N], le[N], sn[N];
    ALIGNED short VC[N];
    ALIGNED short cmp[N];
    register int i;

#ifdef ARCH_MIN_ARM_NEON

	int16x8_t v_sn, v_cmp, v_vc;

	int16x8_t zero = vdupq_n_s16(0);
	
	int16x8_t vs = vld1q_s16((const int16_t *)VS);
    int16x8_t vt = vld1q_s16((const int16_t *)VT);

	v_vc = vt;
	v_sn = veorq_s16(vs,vt);
	v_sn = vshrq_n_s16(v_sn,15);

	v_cmp = vandq_s16(vs, v_sn);
	v_cmp = vmvnq_s16(v_cmp);
	uint16x8_t v_le = vcleq_s16(vt,v_cmp);
	int16x8_t v_le_ = vnegq_s16((int16x8_t)v_le);

	v_cmp = vorrq_s16(vs, v_sn);
	uint16x8_t v_ge = vcgeq_s16(v_cmp, vt);
	int16x8_t v_ge_ = vnegq_s16((int16x8_t)v_ge);

	v_vc = veorq_s16(v_vc,v_sn);

	vst1q_s16(VC, v_vc);
	vst1q_s16(le, v_le_);
	vst1q_s16(ge, v_ge_);

	merge(VACC_L, le, VC, VS);
    vector_copy(VD, VACC_L);

	vst1q_s16(state->clip, v_ge_);
	vst1q_s16(state->comp, v_le_);
	vst1q_s16(state->ne,zero);
	vst1q_s16(state->co,zero);
	vst1q_s16(state->vce,zero);
	
	return;
	
#else
	
	
    for (i = 0; i < N; i++)
        VC[i] = VT[i];
    for (i = 0; i < N; i++)
        sn[i] = (signed short)(VS[i] ^ VT[i]) >> 15;
#if (0)
    for (i = 0; i < N; i++)
        le[i] = sn[i] ? (VT[i] <= ~VS[i]) : (VT[i] <= ~0x0000);
    for (i = 0; i < N; i++)
        ge[i] = sn[i] ? (~0x0000 >= VT[i]) : (VS[i] >= VT[i]);
#else
    for (i = 0; i < N; i++)
        cmp[i] = ~(VS[i] & sn[i]);
    for (i = 0; i < N; i++)
        le[i] = (VT[i] <= cmp[i]);
    for (i = 0; i < N; i++)
        cmp[i] =  (VS[i] | sn[i]);
    for (i = 0; i < N; i++)
        ge[i] = (cmp[i] >= VT[i]);
#endif
    for (i = 0; i < N; i++)
        VC[i] ^= sn[i]; /* if (sn == ~0) {VT = ~VT;} else {VT =  VT;} */
    merge(VACC_L, le, VC, VS);
    vector_copy(VD, VACC_L);

    for (i = 0; i < N; i++)
        state->clip[i] = ge[i];
    for (i = 0; i < N; i++)
        state->comp[i] = le[i];
    for (i = 0; i < N; i++)
        state->ne[i] = 0;
    for (i = 0; i < N; i++)
        state->co[i] = 0;
    for (i = 0; i < N; i++)
        state->vce[i] = 0;
    return;
#endif
}

static void VCR(usf_state_t * state, int vd, int vs, int vt, int e)
{
    ALIGNED short ST[N];

    SHUFFLE_VECTOR(ST, state->VR[vt], e);
    do_cr(state, state->VR[vd], state->VR[vs], ST);
    return;
}