db4d3af31d
This includes the following commits: changeset: 22:14a375563e2b summary: chang "while" command to "type -P" in bash based configure changeset: 21:c4bf7c0f33d2 summary: updated codec paper for RS01 life cycle changeset: 20:dbcaad8128aa summary: replaced build count with mercurial changeset changeset: 19:fba258a6acfa summary: Added tag 0.79.10 for changeset 49950be5a2ef changeset: 18:49950be5a2ef summary: merged some debian patches changeset: 17:815be3929c41 summary: merged "easy" patches from Stéphane Lesimple's version changeset: 16:7d15f8a958cb summary: Made printf format strings 32/64bit safe as suggested by Stéphane; changeset: 15:1055a53b8d6d summary: reorganized code for --with-gui=[yes|no] option changeset: 14:fbe2ae12a32c summary: Added tag 0.79.9 for changeset f2fdd6d3a1f5 changeset: 13:f2fdd6d3a1f5 summary: updated TODO and CHANGELOG And other changes that were needed to resolve the (many) conflicts.
115 lines
3.0 KiB
C
115 lines
3.0 KiB
C
/* dvdisaster: Additional error correction for optical media.
|
|
* Copyright (C) 2004-2017 Carsten Gnoerlich.
|
|
* Copyright (C) 2019-2021 The dvdisaster development team.
|
|
*
|
|
* Email: support@dvdisaster.org
|
|
*
|
|
* This file is part of dvdisaster.
|
|
*
|
|
* dvdisaster is free software: you can redistribute it and/or modify
|
|
* it under the terms of the GNU General Public License as published by
|
|
* the Free Software Foundation, either version 3 of the License, or
|
|
* (at your option) any later version.
|
|
*
|
|
* dvdisaster is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
* GNU General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU General Public License
|
|
* along with dvdisaster. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
/*** src type: no GUI code ***/
|
|
|
|
#include "dvdisaster.h"
|
|
|
|
#ifdef HAVE_SSE2
|
|
#include <emmintrin.h>
|
|
|
|
#ifdef HAVE_CPUID
|
|
#include <cpuid.h>
|
|
#else
|
|
#include "compat/cpuid.h"
|
|
#endif
|
|
#endif
|
|
|
|
/***
|
|
*** Reed-Solomon encoding using SSE2 intrinsics
|
|
***/
|
|
|
|
/* SSE 2 version */
|
|
|
|
#ifdef HAVE_SSE2
|
|
int ProbeSSE2(void)
|
|
{ unsigned int eax, ebx, ecx, edx;
|
|
|
|
if(!__get_cpuid(1, &eax, &ebx, &ecx, &edx))
|
|
{ Verbose("[ProbeSSE2: get_cpuid() failed]\n");
|
|
return 0;
|
|
}
|
|
|
|
if(edx & bit_SSE2)
|
|
{ Verbose("[ProbeSSE2: SSE2 available]\n");
|
|
return 1;
|
|
}
|
|
else
|
|
{ Verbose("[ProbeSSE2: no SSE2]\n");
|
|
return 0;
|
|
}
|
|
}
|
|
|
|
void encode_next_layer_sse2(ReedSolomonTables *rt, unsigned char *data, unsigned char *parity, guint64 layer_size, int shift)
|
|
{ gint32 *gf_index_of = rt->gfTables->indexOf;
|
|
gint32 *enc_alpha_to = rt->gfTables->encAlphaTo;
|
|
gint32 *rs_gpoly = rt->gpoly;
|
|
int nroots = rt->nroots;
|
|
int nroots_aligned = (nroots+15)&~15;
|
|
int nroots_full = nroots_aligned>>4;
|
|
int i,j;
|
|
|
|
for(i=0; i<layer_size; i++)
|
|
{ int feedback = gf_index_of[data[i] ^ parity[shift]];
|
|
int offset = nroots-shift-1;
|
|
|
|
if(feedback != GF_ALPHA0) /* non-zero feedback term */
|
|
{ guint8 *par_idx = (guint8*)parity;
|
|
guint8 *e_lut = rt->bLut[feedback]+offset;
|
|
__m128i par, lut, out;
|
|
|
|
/* Process lut in 128 bit steps */
|
|
|
|
for(j=nroots_full; j; j--)
|
|
{
|
|
par = _mm_load_si128((__m128i*)par_idx);
|
|
lut = _mm_loadu_si128((__m128i*)e_lut);
|
|
out = _mm_xor_si128(par, lut);
|
|
_mm_store_si128((__m128i*)par_idx, out);
|
|
par_idx += 16;
|
|
e_lut += 16;
|
|
}
|
|
|
|
parity[shift] = enc_alpha_to[feedback + rs_gpoly[0]];
|
|
}
|
|
else /* zero feedback term */
|
|
parity[shift] = 0;
|
|
|
|
parity += nroots_aligned;
|
|
}
|
|
}
|
|
#else /* don't have SSE2 */
|
|
/* Stub functions to keep the linker happy.
|
|
* Should never be executed.
|
|
*/
|
|
|
|
int ProbeSSE2()
|
|
{ return 0;
|
|
}
|
|
|
|
void encode_next_layer_sse2(ReedSolomonTables *rt, unsigned char *data, unsigned char *parity, guint64 layer_size, int shift)
|
|
{
|
|
Stop("Mega borkage - EncodeNextLayerSSE2() stub called.\n");
|
|
}
|
|
#endif /* HAVE_SSE2 */
|
|
|