cog/Frameworks/lazyusf/lazyusf/rsp/vu/vmulf.h

56 lines
2.1 KiB
C

/******************************************************************************\
* 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"
#ifndef SEMIFRAC
/*
* acc = VS * VT;
* acc = acc + 0x8000; // rounding value
* acc = acc << 1; // partial value shifting
*
* Wrong: ACC(HI) = -((INT32)(acc) < 0)
* Right: ACC(HI) = -(SEMIFRAC < 0)
*/
#define SEMIFRAC (VS[i]*VT[i]*2/2 + 0x8000/2)
#endif
INLINE static void do_mulf(usf_state_t * state, short* VD, short* VS, short* VT)
{
register int i;
for (i = 0; i < N; i++)
VACC_L[i] = (SEMIFRAC << 1) >> 0;
for (i = 0; i < N; i++)
VACC_M[i] = (SEMIFRAC << 1) >> 16;
for (i = 0; i < N; i++)
VACC_H[i] = -((VACC_M[i] < 0) & (VS[i] != VT[i])); /* -32768 * -32768 */
#ifndef ARCH_MIN_SSE2
vector_copy(VD, VACC_M);
for (i = 0; i < N; i++)
VD[i] -= (VACC_M[i] < 0) & (VS[i] == VT[i]); /* ACC b 31 set, min*min */
#else
SIGNED_CLAMP_AM(state, VD);
#endif
return;
}
static void VMULF(usf_state_t * state, int vd, int vs, int vt, int e)
{
ALIGNED short ST[N];
SHUFFLE_VECTOR(ST, state->VR[vt], e);
do_mulf(state, state->VR[vd], state->VR[vs], ST);
return;
}