aboutsummaryrefslogtreecommitdiffstats
path: root/libArithCoding/src/ac_arith_coder.cpp
diff options
context:
space:
mode:
authorFraunhofer IIS FDK <audio-fdk@iis.fraunhofer.de>2018-02-26 20:17:00 +0100
committerJean-Michel Trivi <jmtrivi@google.com>2018-04-19 11:21:15 -0700
commit6cfabd35363c3ef5e3b209b867169a500b3ccc3c (patch)
tree01c0a19f2735e8b5d2407555fe992d4230d089eb /libArithCoding/src/ac_arith_coder.cpp
parent6288a1e34c4dede4c2806beb1736ece6580558c7 (diff)
downloadfdk-aac-6cfabd35363c3ef5e3b209b867169a500b3ccc3c.tar.gz
fdk-aac-6cfabd35363c3ef5e3b209b867169a500b3ccc3c.tar.bz2
fdk-aac-6cfabd35363c3ef5e3b209b867169a500b3ccc3c.zip
Upgrade to FDKv2
Bug: 71430241 Test: CTS DecoderTest and DecoderTestAacDrc original-Change-Id: Iaa20f749b8a04d553b20247cfe1a8930ebbabe30 Apply clang-format also on header files. original-Change-Id: I14de1ef16bbc79ec0283e745f98356a10efeb2e4 Fixes for MPEG-D DRC original-Change-Id: If1de2d74bbbac84b3f67de3b88b83f6a23b8a15c Catch unsupported tw_mdct at an early stage original-Change-Id: Ied9dd00d754162a0e3ca1ae3e6b854315d818afe Fixing PVC transition frames original-Change-Id: Ib75725abe39252806c32d71176308f2c03547a4e Move qmf bands sanity check original-Change-Id: Iab540c3013c174d9490d2ae100a4576f51d8dbc4 Initialize scaling variable original-Change-Id: I3c4087101b70e998c71c1689b122b0d7762e0f9e Add 16 qmf band configuration to getSlotNrgHQ() original-Change-Id: I49a5d30f703a1b126ff163df9656db2540df21f1 Always apply byte alignment at the end of the AudioMuxElement original-Change-Id: I42d560287506d65d4c3de8bfe3eb9a4ebeb4efc7 Setup SBR element only if no parse error exists original-Change-Id: I1915b73704bc80ab882b9173d6bec59cbd073676 Additional array index check in HCR original-Change-Id: I18cc6e501ea683b5009f1bbee26de8ddd04d8267 Fix fade-in index selection in concealment module original-Change-Id: Ibf802ed6ed8c05e9257e1f3b6d0ac1162e9b81c1 Enable explicit backward compatible parser for AAC_LD original-Change-Id: I27e9c678dcb5d40ed760a6d1e06609563d02482d Skip spatial specific config in explicit backward compatible ASC original-Change-Id: Iff7cc365561319e886090cedf30533f562ea4d6e Update flags description in decoder API original-Change-Id: I9a5b4f8da76bb652f5580cbd3ba9760425c43830 Add QMF domain reset function original-Change-Id: I4f89a8a2c0277d18103380134e4ed86996e9d8d6 DRC upgrade v2.1.0 original-Change-Id: I5731c0540139dab220094cd978ef42099fc45b74 Fix integer overflow in sqrtFixp_lookup() original-Change-Id: I429a6f0d19aa2cc957e0f181066f0ca73968c914 Fix integer overflow in invSqrtNorm2() original-Change-Id: I84de5cbf9fb3adeb611db203fe492fabf4eb6155 Fix integer overflow in GenerateRandomVector() original-Change-Id: I3118a641008bd9484d479e5b0b1ee2b5d7d44d74 Fix integer overflow in adjustTimeSlot_EldGrid() original-Change-Id: I29d503c247c5c8282349b79df940416a512fb9d5 Fix integer overflow in FDKsbrEnc_codeEnvelope() original-Change-Id: I6b34b61ebb9d525b0c651ed08de2befc1f801449 Follow-up on: Fix integer overflow in adjustTimeSlot_EldGrid() original-Change-Id: I6f8f578cc7089e5eb7c7b93e580b72ca35ad689a Fix integer overflow in get_pk_v2() original-Change-Id: I63375bed40d45867f6eeaa72b20b1f33e815938c Fix integer overflow in Syn_filt_zero() original-Change-Id: Ie0c02fdfbe03988f9d3b20d10cd9fe4c002d1279 Fix integer overflow in CFac_CalcFacSignal() original-Change-Id: Id2d767c40066c591b51768e978eb8af3b803f0c5 Fix integer overflow in FDKaacEnc_FDKaacEnc_calcPeNoAH() original-Change-Id: Idcbd0f4a51ae2550ed106aa6f3d678d1f9724841 Fix integer overflow in sbrDecoder_calculateGainVec() original-Change-Id: I7081bcbe29c5cede9821b38d93de07c7add2d507 Fix integer overflow in CLpc_SynthesisLattice() original-Change-Id: I4a95ddc18de150102352d4a1845f06094764c881 Fix integer overflow in Pred_Lt4() original-Change-Id: I4dbd012b2de7d07c3e70a47b92e3bfae8dbc750a Fix integer overflow in FDKsbrEnc_InitSbrFastTransientDetector() original-Change-Id: I788cbec1a4a00f44c2f3a72ad7a4afa219807d04 Fix unsigned integer overflow in FDKaacEnc_WriteBitstream() original-Change-Id: I68fc75166e7d2cd5cd45b18dbe3d8c2a92f1822a Fix unsigned integer overflow in FDK_MetadataEnc_Init() original-Change-Id: Ie8d025f9bcdb2442c704bd196e61065c03c10af4 Fix overflow in pseudo random number generators original-Change-Id: I3e2551ee01356297ca14e3788436ede80bd5513c Fix unsigned integer overflow in sbrDecoder_Parse() original-Change-Id: I3f231b2f437e9c37db4d5b964164686710eee971 Fix unsigned integer overflow in longsub() original-Change-Id: I73c2bc50415cac26f1f5a29e125bbe75f9180a6e Fix unsigned integer overflow in CAacDecoder_DecodeFrame() original-Change-Id: Ifce2db4b1454b46fa5f887e9d383f1cc43b291e4 Fix overflow at CLpdChannelStream_Read() original-Change-Id: Idb9d822ce3a4272e4794b643644f5434e2d4bf3f Fix unsigned integer overflow in Hcr_State_BODY_SIGN_ESC__ESC_WORD() original-Change-Id: I1ccf77c0015684b85534c5eb97162740a870b71c Fix unsigned integer overflow in UsacConfig_Parse() original-Change-Id: Ie6d27f84b6ae7eef092ecbff4447941c77864d9f Fix unsigned integer overflow in aacDecoder_drcParse() original-Change-Id: I713f28e883eea3d70b6fa56a7b8f8c22bcf66ca0 Fix unsigned integer overflow in aacDecoder_drcReadCompression() original-Change-Id: Ia34dfeb88c4705c558bce34314f584965cafcf7a Fix unsigned integer overflow in CDataStreamElement_Read() original-Change-Id: Iae896cc1d11f0a893d21be6aa90bd3e60a2c25f0 Fix unsigned integer overflow in transportDec_AdjustEndOfAccessUnit() original-Change-Id: I64cf29a153ee784bb4a16fdc088baabebc0007dc Fix unsigned integer overflow in transportDec_GetAuBitsRemaining() original-Change-Id: I975b3420faa9c16a041874ba0db82e92035962e4 Fix unsigned integer overflow in extractExtendedData() original-Change-Id: I2a59eb09e2053cfb58dfb75fcecfad6b85a80a8f Fix signed integer overflow in CAacDecoder_ExtPayloadParse() original-Change-Id: I4ad5ca4e3b83b5d964f1c2f8c5e7b17c477c7929 Fix unsigned integer overflow in CAacDecoder_DecodeFrame() original-Change-Id: I29a39df77d45c52a0c9c5c83c1ba81f8d0f25090 Follow-up on: Fix integer overflow in CLpc_SynthesisLattice() original-Change-Id: I8fb194ffc073a3432a380845be71036a272d388f Fix signed integer overflow in _interpolateDrcGain() original-Change-Id: I879ec9ab14005069a7c47faf80e8bc6e03d22e60 Fix unsigned integer overflow in FDKreadBits() original-Change-Id: I1f47a6a8037ff70375aa8844947d5681bb4287ad Fix unsigned integer overflow in FDKbyteAlign() original-Change-Id: Id5f3a11a0c9e50fc6f76ed6c572dbd4e9f2af766 Fix unsigned integer overflow in FDK_get32() original-Change-Id: I9d33b8e97e3d38cbb80629cb859266ca0acdce96 Fix unsigned integer overflow in FDK_pushBack() original-Change-Id: Ic87f899bc8c6acf7a377a8ca7f3ba74c3a1e1c19 Fix unsigned integer overflow in FDK_pushForward() original-Change-Id: I3b754382f6776a34be1602e66694ede8e0b8effc Fix unsigned integer overflow in ReadPsData() original-Change-Id: I25361664ba8139e32bbbef2ca8c106a606ce9c37 Fix signed integer overflow in E_UTIL_residu() original-Change-Id: I8c3abd1f437ee869caa8fb5903ce7d3d641b6aad REVERT: Follow-up on: Integer overflow in CLpc_SynthesisLattice(). original-Change-Id: I3d340099acb0414795c8dfbe6362bc0a8f045f9b Follow-up on: Fix integer overflow in CLpc_SynthesisLattice() original-Change-Id: I4aedb8b3a187064e9f4d985175aa55bb99cc7590 Follow-up on: Fix unsigned integer overflow in aacDecoder_drcParse() original-Change-Id: I2aa2e13916213bf52a67e8b0518e7bf7e57fb37d Fix integer overflow in acelp original-Change-Id: Ie6390c136d84055f8b728aefbe4ebef6e029dc77 Fix unsigned integer overflow in aacDecoder_UpdateBitStreamCounters() original-Change-Id: I391ffd97ddb0b2c184cba76139bfb356a3b4d2e2 Adjust concealment default settings original-Change-Id: I6a95db935a327c47df348030bcceafcb29f54b21 Saturate estimatedStartPos original-Change-Id: I27be2085e0ae83ec9501409f65e003f6bcba1ab6 Negative shift exponent in _interpolateDrcGain() original-Change-Id: I18edb26b26d002aafd5e633d4914960f7a359c29 Negative shift exponent in calculateICC() original-Change-Id: I3dcd2ae98d2eb70ee0d59750863cbb2a6f4f8aba Too large shift exponent in FDK_put() original-Change-Id: Ib7d9aaa434d2d8de4a13b720ca0464b31ca9b671 Too large shift exponent in CalcInvLdData() original-Change-Id: I43e6e78d4cd12daeb1dcd5d82d1798bdc2550262 Member access within null pointer of type SBR_CHANNEL original-Change-Id: Idc5e4ea8997810376d2f36bbdf628923b135b097 Member access within null pointer of type CpePersistentData original-Change-Id: Ib6c91cb0d37882768e5baf63324e429589de0d9d Member access within null pointer FDKaacEnc_psyMain() original-Change-Id: I7729b7f4479970531d9dc823abff63ca52e01997 Member access within null pointer FDKaacEnc_GetPnsParam() original-Change-Id: I9aa3b9f3456ae2e0f7483dbd5b3dde95fc62da39 Member access within null pointer FDKsbrEnc_EnvEncodeFrame() original-Change-Id: I67936f90ea714e90b3e81bc0dd1472cc713eb23a Add HCR sanity check original-Change-Id: I6c1d9732ebcf6af12f50b7641400752f74be39f7 Fix memory issue for HBE edge case with 8:3 SBR original-Change-Id: I11ea58a61e69fbe8bf75034b640baee3011e63e9 Additional SBR parametrization sanity check for ELD original-Change-Id: Ie26026fbfe174c2c7b3691f6218b5ce63e322140 Add MPEG-D DRC channel layout check original-Change-Id: Iea70a74f171b227cce636a9eac4ba662777a2f72 Additional out-of-bounds checks in MPEG-D DRC original-Change-Id: Ife4a8c3452c6fde8a0a09e941154a39a769777d4 Change-Id: Ic63cb2f628720f54fe9b572b0cb528e2599c624e
Diffstat (limited to 'libArithCoding/src/ac_arith_coder.cpp')
-rw-r--r--libArithCoding/src/ac_arith_coder.cpp777
1 files changed, 777 insertions, 0 deletions
diff --git a/libArithCoding/src/ac_arith_coder.cpp b/libArithCoding/src/ac_arith_coder.cpp
new file mode 100644
index 0000000..b791f39
--- /dev/null
+++ b/libArithCoding/src/ac_arith_coder.cpp
@@ -0,0 +1,777 @@
+/* -----------------------------------------------------------------------------
+Software License for The Fraunhofer FDK AAC Codec Library for Android
+
+© Copyright 1995 - 2018 Fraunhofer-Gesellschaft zur Förderung der angewandten
+Forschung e.V. All rights reserved.
+
+ 1. INTRODUCTION
+The Fraunhofer FDK AAC Codec Library for Android ("FDK AAC Codec") is software
+that implements the MPEG Advanced Audio Coding ("AAC") encoding and decoding
+scheme for digital audio. This FDK AAC Codec software is intended to be used on
+a wide variety of Android devices.
+
+AAC's HE-AAC and HE-AAC v2 versions are regarded as today's most efficient
+general perceptual audio codecs. AAC-ELD is considered the best-performing
+full-bandwidth communications codec by independent studies and is widely
+deployed. AAC has been standardized by ISO and IEC as part of the MPEG
+specifications.
+
+Patent licenses for necessary patent claims for the FDK AAC Codec (including
+those of Fraunhofer) may be obtained through Via Licensing
+(www.vialicensing.com) or through the respective patent owners individually for
+the purpose of encoding or decoding bit streams in products that are compliant
+with the ISO/IEC MPEG audio standards. Please note that most manufacturers of
+Android devices already license these patent claims through Via Licensing or
+directly from the patent owners, and therefore FDK AAC Codec software may
+already be covered under those patent licenses when it is used for those
+licensed purposes only.
+
+Commercially-licensed AAC software libraries, including floating-point versions
+with enhanced sound quality, are also available from Fraunhofer. Users are
+encouraged to check the Fraunhofer website for additional applications
+information and documentation.
+
+2. COPYRIGHT LICENSE
+
+Redistribution and use in source and binary forms, with or without modification,
+are permitted without payment of copyright license fees provided that you
+satisfy the following conditions:
+
+You must retain the complete text of this software license in redistributions of
+the FDK AAC Codec or your modifications thereto in source code form.
+
+You must retain the complete text of this software license in the documentation
+and/or other materials provided with redistributions of the FDK AAC Codec or
+your modifications thereto in binary form. You must make available free of
+charge copies of the complete source code of the FDK AAC Codec and your
+modifications thereto to recipients of copies in binary form.
+
+The name of Fraunhofer may not be used to endorse or promote products derived
+from this library without prior written permission.
+
+You may not charge copyright license fees for anyone to use, copy or distribute
+the FDK AAC Codec software or your modifications thereto.
+
+Your modified versions of the FDK AAC Codec must carry prominent notices stating
+that you changed the software and the date of any change. For modified versions
+of the FDK AAC Codec, the term "Fraunhofer FDK AAC Codec Library for Android"
+must be replaced by the term "Third-Party Modified Version of the Fraunhofer FDK
+AAC Codec Library for Android."
+
+3. NO PATENT LICENSE
+
+NO EXPRESS OR IMPLIED LICENSES TO ANY PATENT CLAIMS, including without
+limitation the patents of Fraunhofer, ARE GRANTED BY THIS SOFTWARE LICENSE.
+Fraunhofer provides no warranty of patent non-infringement with respect to this
+software.
+
+You may use this FDK AAC Codec software or modifications thereto only for
+purposes that are authorized by appropriate patent licenses.
+
+4. DISCLAIMER
+
+This FDK AAC Codec software is provided by Fraunhofer on behalf of the copyright
+holders and contributors "AS IS" and WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES,
+including but not limited to the implied warranties of merchantability and
+fitness for a particular purpose. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
+CONTRIBUTORS BE LIABLE for any direct, indirect, incidental, special, exemplary,
+or consequential damages, including but not limited to procurement of substitute
+goods or services; loss of use, data, or profits, or business interruption,
+however caused and on any theory of liability, whether in contract, strict
+liability, or tort (including negligence), arising in any way out of the use of
+this software, even if advised of the possibility of such damage.
+
+5. CONTACT INFORMATION
+
+Fraunhofer Institute for Integrated Circuits IIS
+Attention: Audio and Multimedia Departments - FDK AAC LL
+Am Wolfsmantel 33
+91058 Erlangen, Germany
+
+www.iis.fraunhofer.de/amm
+amm-info@iis.fraunhofer.de
+----------------------------------------------------------------------------- */
+
+/************************** Arithmetic coder library ***************************
+
+ Author(s): Youliy Ninov, Oliver Weiss
+
+ Description: Definition of Spectral Noiseless Coding Scheme based on an
+ Arithmetic Coder in Conjunction with an Adaptive Context
+
+*******************************************************************************/
+
+#include "ac_arith_coder.h"
+
+#define cbitsnew 16
+#define stat_bitsnew 14
+#define ari_q4new (((long)1 << cbitsnew) - 1) /* 0xFFFF */
+#define ari_q1new (ari_q4new / 4 + 1) /* 0x4000 */
+#define ari_q2new (2 * ari_q1new) /* 0x8000 */
+#define ari_q3new (3 * ari_q1new) /* 0xC000 */
+
+#define VAL_ESC 16
+
+/* Arithmetic coder library info */
+#define AC_LIB_VL0 2
+#define AC_LIB_VL1 0
+#define AC_LIB_VL2 0
+#define AC_LIB_TITLE "Arithmetic Coder Lib"
+#ifdef __ANDROID__
+#define AC_LIB_BUILD_DATE ""
+#define AC_LIB_BUILD_TIME ""
+#else
+#define AC_LIB_BUILD_DATE __DATE__
+#define AC_LIB_BUILD_TIME __TIME__
+#endif
+
+const SHORT ari_lsb2[3][4] = {
+ {12571, 10569, 3696, 0}, {12661, 5700, 3751, 0}, {10827, 6884, 2929, 0}};
+
+H_ALLOC_MEM(ArcoData, CArcoData)
+/*! The structure ArcoData contains 2-tuple context of previous frame. <br>
+ Dimension: 1 */
+C_ALLOC_MEM(ArcoData, CArcoData, 1)
+
+/*
+ This define triggers the use of the pre-known return values of function
+ get_pk_v2() for the cases, where parameter s is in range
+ 0x00000000..0x0000000F. Note: These 16 bytes have been moved into the first 4
+ entries of ari_merged_hash_ps that are no more referenced.
+*/
+
+static const ULONG ari_merged_hash_ps[742] = {
+ 0x00001044UL, 0x00003D0AUL, 0x00005350UL, 0x000074D6UL, 0x0000A49FUL,
+ 0x0000F96EUL, 0x00111000UL, 0x01111E83UL, 0x01113146UL, 0x01114036UL,
+ 0x01116863UL, 0x011194E9UL, 0x0111F7EEUL, 0x0112269BUL, 0x01124775UL,
+ 0x01126DA1UL, 0x0112D912UL, 0x01131AF0UL, 0x011336DDUL, 0x01135CF5UL,
+ 0x01139DF8UL, 0x01141A5BUL, 0x01144773UL, 0x01146CF5UL, 0x0114FDE9UL,
+ 0x01166CF5UL, 0x0116FDE4UL, 0x01174CF3UL, 0x011FFDCFUL, 0x01211CC2UL,
+ 0x01213B2DUL, 0x01214036UL, 0x01216863UL, 0x012194D2UL, 0x0122197FUL,
+ 0x01223AADUL, 0x01224036UL, 0x01226878UL, 0x0122A929UL, 0x0122F4ABUL,
+ 0x01232B2DUL, 0x012347B6UL, 0x01237DF8UL, 0x0123B929UL, 0x012417DDUL,
+ 0x01245D76UL, 0x01249DF8UL, 0x0124F912UL, 0x01255D75UL, 0x0125FDE9UL,
+ 0x01265D75UL, 0x012B8DF7UL, 0x01311E2AUL, 0x01313B5EUL, 0x0131687BUL,
+ 0x01321A6DUL, 0x013237BCUL, 0x01326863UL, 0x0132F4EEUL, 0x01332B5EUL,
+ 0x01335DA1UL, 0x01338E24UL, 0x01341A5EUL, 0x01343DB6UL, 0x01348DF8UL,
+ 0x01351935UL, 0x01355DB7UL, 0x0135FE12UL, 0x01376DF7UL, 0x013FFE29UL,
+ 0x01400024UL, 0x01423821UL, 0x014318F6UL, 0x01433821UL, 0x0143F8E5UL,
+ 0x01443DA1UL, 0x01486E38UL, 0x014FF929UL, 0x01543EE3UL, 0x015FF912UL,
+ 0x016F298CUL, 0x018A5A69UL, 0x021007F1UL, 0x02112C2CUL, 0x02114B48UL,
+ 0x02117353UL, 0x0211F4AEUL, 0x02122FEAUL, 0x02124B48UL, 0x02127850UL,
+ 0x0212F72EUL, 0x02133A9EUL, 0x02134036UL, 0x02138864UL, 0x021414ADUL,
+ 0x0214379EUL, 0x02145DB6UL, 0x0214FE1FUL, 0x02166DB7UL, 0x02200DC4UL,
+ 0x02212FEAUL, 0x022147A0UL, 0x02218369UL, 0x0221F7EEUL, 0x02222AADUL,
+ 0x02224788UL, 0x02226863UL, 0x02229929UL, 0x0222F4ABUL, 0x02232A9EUL,
+ 0x02234F08UL, 0x02237864UL, 0x0223A929UL, 0x022417DEUL, 0x02244F36UL,
+ 0x02248863UL, 0x02251A74UL, 0x02256DA1UL, 0x0225FE12UL, 0x02263DB6UL,
+ 0x02276DF7UL, 0x022FFE29UL, 0x0231186DUL, 0x023137BCUL, 0x02314020UL,
+ 0x02319ED6UL, 0x0232196DUL, 0x023237BCUL, 0x02325809UL, 0x02329429UL,
+ 0x023317EDUL, 0x02333F08UL, 0x02335809UL, 0x023378E4UL, 0x02341A7CUL,
+ 0x02344221UL, 0x0234A8D3UL, 0x023514BCUL, 0x02354221UL, 0x0235F8DFUL,
+ 0x02364861UL, 0x023FFE29UL, 0x02400024UL, 0x0241583BUL, 0x024214C8UL,
+ 0x02424809UL, 0x02427EE6UL, 0x02431708UL, 0x02434809UL, 0x02436ED0UL,
+ 0x02441A76UL, 0x02443821UL, 0x024458E3UL, 0x0244F91FUL, 0x02454863UL,
+ 0x0246190AUL, 0x02464863UL, 0x024FF929UL, 0x02525ED0UL, 0x025314E1UL,
+ 0x025348FBUL, 0x025419A1UL, 0x025458D0UL, 0x0254F4E5UL, 0x02552861UL,
+ 0x025FF912UL, 0x02665993UL, 0x027F5A69UL, 0x029F1481UL, 0x02CF28A4UL,
+ 0x03100AF0UL, 0x031120AAUL, 0x031147A0UL, 0x03118356UL, 0x031217ECUL,
+ 0x03123B5EUL, 0x03124008UL, 0x03127350UL, 0x031314AAUL, 0x0313201EUL,
+ 0x03134F08UL, 0x03136863UL, 0x03141A5EUL, 0x03143F3CUL, 0x03200847UL,
+ 0x03212AADUL, 0x03214F20UL, 0x03218ED6UL, 0x032218ADUL, 0x032237BCUL,
+ 0x03225809UL, 0x03229416UL, 0x032317EDUL, 0x03234F20UL, 0x03237350UL,
+ 0x0323FA6BUL, 0x03243F08UL, 0x03246863UL, 0x0324F925UL, 0x03254221UL,
+ 0x0325F8DFUL, 0x03264821UL, 0x032FFE29UL, 0x03311E47UL, 0x03313F08UL,
+ 0x0331580DUL, 0x033214DEUL, 0x03323F08UL, 0x03324020UL, 0x03326350UL,
+ 0x033294E9UL, 0x033317DEUL, 0x03333F08UL, 0x0333627BUL, 0x0333A9A9UL,
+ 0x033417FCUL, 0x03343220UL, 0x0334627BUL, 0x0334A9A9UL, 0x0335148AUL,
+ 0x03353220UL, 0x033588E4UL, 0x03361A4AUL, 0x03363821UL, 0x0336F8D2UL,
+ 0x03376863UL, 0x03411939UL, 0x0341583BUL, 0x034214C8UL, 0x03424809UL,
+ 0x03426ED0UL, 0x03431588UL, 0x03434809UL, 0x03436ED0UL, 0x03441A48UL,
+ 0x0344480DUL, 0x03446ED0UL, 0x03451A4AUL, 0x03453809UL, 0x03455EFBUL,
+ 0x034614CAUL, 0x03463849UL, 0x034F8924UL, 0x03500A69UL, 0x035252D0UL,
+ 0x035314E0UL, 0x0353324DUL, 0x03535ED0UL, 0x035414E0UL, 0x0354324DUL,
+ 0x03545ED0UL, 0x0354F4E8UL, 0x0355384DUL, 0x03555ED0UL, 0x0355F4DFUL,
+ 0x03564350UL, 0x035969A6UL, 0x035FFA52UL, 0x036649A6UL, 0x036FFA52UL,
+ 0x037F4F66UL, 0x039D7492UL, 0x03BF6892UL, 0x03DF8A1FUL, 0x04100B84UL,
+ 0x04112107UL, 0x0411520DUL, 0x041214EAUL, 0x04124F20UL, 0x04131EDEUL,
+ 0x04133F08UL, 0x04135809UL, 0x0413F42BUL, 0x04142F08UL, 0x04200847UL,
+ 0x042121FCUL, 0x04214209UL, 0x04221407UL, 0x0422203CUL, 0x04224209UL,
+ 0x04226350UL, 0x04231A7CUL, 0x04234209UL, 0x0423637BUL, 0x04241A7CUL,
+ 0x04243220UL, 0x0424627BUL, 0x042514C8UL, 0x04254809UL, 0x042FF8E9UL,
+ 0x04311E47UL, 0x04313220UL, 0x0431527BUL, 0x043214FCUL, 0x04323220UL,
+ 0x04326250UL, 0x043315BCUL, 0x04333220UL, 0x0433527BUL, 0x04338413UL,
+ 0x04341488UL, 0x04344809UL, 0x04346ED0UL, 0x04351F48UL, 0x0435527BUL,
+ 0x0435F9A5UL, 0x04363809UL, 0x04375EFBUL, 0x043FF929UL, 0x04412E79UL,
+ 0x0441427BUL, 0x044219B9UL, 0x04423809UL, 0x0442537BUL, 0x044314C8UL,
+ 0x04432020UL, 0x0443527BUL, 0x044414CAUL, 0x04443809UL, 0x0444537BUL,
+ 0x04448993UL, 0x0445148AUL, 0x04453809UL, 0x04455ED0UL, 0x0445F4E5UL,
+ 0x0446384DUL, 0x045009A6UL, 0x045272D3UL, 0x045314A0UL, 0x0453324DUL,
+ 0x04535ED0UL, 0x045415A0UL, 0x0454324DUL, 0x04545ED0UL, 0x04551F60UL,
+ 0x0455324DUL, 0x04562989UL, 0x04564350UL, 0x045FF4D2UL, 0x04665993UL,
+ 0x047FFF62UL, 0x048FF725UL, 0x049F44BDUL, 0x04BFB7E5UL, 0x04EF8A25UL,
+ 0x04FFFB98UL, 0x051131F9UL, 0x051212C7UL, 0x05134209UL, 0x05200247UL,
+ 0x05211007UL, 0x05213E60UL, 0x052212C7UL, 0x05224209UL, 0x052319BCUL,
+ 0x05233220UL, 0x0523527BUL, 0x052414C8UL, 0x05243820UL, 0x053112F9UL,
+ 0x05313E49UL, 0x05321439UL, 0x05323E49UL, 0x0532537BUL, 0x053314C8UL,
+ 0x0533480DUL, 0x05337413UL, 0x05341488UL, 0x0534527BUL, 0x0534F4EBUL,
+ 0x05353809UL, 0x05356ED0UL, 0x0535F4E5UL, 0x0536427BUL, 0x054119B9UL,
+ 0x054212F9UL, 0x05423249UL, 0x05426ED3UL, 0x05431739UL, 0x05433249UL,
+ 0x05435ED0UL, 0x0543F4EBUL, 0x05443809UL, 0x05445ED0UL, 0x0544F4E8UL,
+ 0x0545324DUL, 0x054FF992UL, 0x055362D3UL, 0x0553F5ABUL, 0x05544350UL,
+ 0x055514CAUL, 0x0555427BUL, 0x0555F4E5UL, 0x0556327BUL, 0x055FF4D2UL,
+ 0x05665993UL, 0x05774F53UL, 0x059FF728UL, 0x05CC37FDUL, 0x05EFBA28UL,
+ 0x05FFFB98UL, 0x061131F9UL, 0x06121407UL, 0x06133E60UL, 0x061A72E4UL,
+ 0x06211E47UL, 0x06214E4BUL, 0x062214C7UL, 0x06223E60UL, 0x062312F9UL,
+ 0x06233E60UL, 0x063112F9UL, 0x06313E4CUL, 0x063219B9UL, 0x06323E49UL,
+ 0x06331439UL, 0x06333809UL, 0x06336EE6UL, 0x0633F5ABUL, 0x06343809UL,
+ 0x0634F42BUL, 0x0635427BUL, 0x063FF992UL, 0x064342FBUL, 0x0643F4EBUL,
+ 0x0644427BUL, 0x064524C9UL, 0x06655993UL, 0x0666170AUL, 0x066652E6UL,
+ 0x067A6F56UL, 0x0698473DUL, 0x06CF67D2UL, 0x06EF3A26UL, 0x06FFFAD8UL,
+ 0x071131CCUL, 0x07211307UL, 0x07222E79UL, 0x072292DCUL, 0x07234E4BUL,
+ 0x073112F9UL, 0x07322339UL, 0x073632CBUL, 0x073FF992UL, 0x074432CBUL,
+ 0x075549A6UL, 0x0776FF68UL, 0x07774350UL, 0x0788473DUL, 0x07CF4516UL,
+ 0x07EF3A26UL, 0x07FFFAD8UL, 0x08222E79UL, 0x083112F9UL, 0x0834330BUL,
+ 0x0845338BUL, 0x08756F5CUL, 0x0887F725UL, 0x08884366UL, 0x08AF649CUL,
+ 0x08F00898UL, 0x08FFFAD8UL, 0x091111C7UL, 0x0932330BUL, 0x0945338BUL,
+ 0x09774F7DUL, 0x0998C725UL, 0x09996416UL, 0x09EF87E5UL, 0x09FFFAD8UL,
+ 0x0A34330BUL, 0x0A45338BUL, 0x0A77467DUL, 0x0AA9F52BUL, 0x0AAA6416UL,
+ 0x0ABD67DFUL, 0x0AFFFA18UL, 0x0B33330BUL, 0x0B4443A6UL, 0x0B76467DUL,
+ 0x0BB9751FUL, 0x0BBB59BDUL, 0x0BEF5892UL, 0x0BFFFAD8UL, 0x0C221339UL,
+ 0x0C53338EUL, 0x0C76367DUL, 0x0CCAF52EUL, 0x0CCC6996UL, 0x0CFFFA18UL,
+ 0x0D44438EUL, 0x0D64264EUL, 0x0DDCF52EUL, 0x0DDD5996UL, 0x0DFFFA18UL,
+ 0x0E43338EUL, 0x0E68465CUL, 0x0EEE651CUL, 0x0EFFFA18UL, 0x0F33238EUL,
+ 0x0F553659UL, 0x0F8F451CUL, 0x0FAFF8AEUL, 0x0FF00A2EUL, 0x0FFF1ACCUL,
+ 0x0FFF33BDUL, 0x0FFF7522UL, 0x0FFFFAD8UL, 0x10002C72UL, 0x1111103EUL,
+ 0x11121E83UL, 0x11131E9AUL, 0x1121115AUL, 0x11221170UL, 0x112316F0UL,
+ 0x1124175DUL, 0x11311CC2UL, 0x11321182UL, 0x11331D42UL, 0x11411D48UL,
+ 0x11421836UL, 0x11431876UL, 0x11441DF5UL, 0x1152287BUL, 0x12111903UL,
+ 0x1212115AUL, 0x121316F0UL, 0x12211B30UL, 0x12221B30UL, 0x12231B02UL,
+ 0x12311184UL, 0x12321D04UL, 0x12331784UL, 0x12411D39UL, 0x12412020UL,
+ 0x12422220UL, 0x12511D89UL, 0x1252227BUL, 0x1258184AUL, 0x12832992UL,
+ 0x1311171AUL, 0x13121B30UL, 0x1312202CUL, 0x131320AAUL, 0x132120AAUL,
+ 0x132220ADUL, 0x13232FEDUL, 0x13312107UL, 0x13322134UL, 0x13332134UL,
+ 0x13411D39UL, 0x13431E74UL, 0x13441834UL, 0x134812B4UL, 0x1352230BUL,
+ 0x13611E4BUL, 0x136522E4UL, 0x141113C2UL, 0x141211C4UL, 0x143121F9UL,
+ 0x143221F9UL, 0x143321CAUL, 0x14351D34UL, 0x14431E47UL, 0x14441E74UL,
+ 0x144612B4UL, 0x1452230EUL, 0x14551E74UL, 0x1471130EUL, 0x151113C2UL,
+ 0x152121F9UL, 0x153121F9UL, 0x153221F9UL, 0x15331007UL, 0x15522E4EUL,
+ 0x15551E74UL, 0x1571130EUL, 0x161113C7UL, 0x162121F9UL, 0x163121F9UL,
+ 0x16611E79UL, 0x16661334UL, 0x171113C7UL, 0x172121F9UL, 0x17451E47UL,
+ 0x1771130CUL, 0x181113C7UL, 0x18211E47UL, 0x18511E4CUL, 0x1882130CUL,
+ 0x191113C7UL, 0x19331E79UL, 0x1A111307UL, 0x1A311E79UL, 0x1F52230EUL,
+ 0x200003C1UL, 0x20001027UL, 0x20004467UL, 0x200079E7UL, 0x2000E5EFUL,
+ 0x21100BC0UL, 0x211129C0UL, 0x21114011UL, 0x211189E7UL, 0x2111F5EFUL,
+ 0x21124011UL, 0x21127455UL, 0x211325C0UL, 0x21134011UL, 0x21137455UL,
+ 0x211425C0UL, 0x21212440UL, 0x21213001UL, 0x2121F9EFUL, 0x21222540UL,
+ 0x21226455UL, 0x2122F5EFUL, 0x21233051UL, 0x2123F56FUL, 0x21244451UL,
+ 0x21312551UL, 0x21323451UL, 0x21332551UL, 0x21844555UL, 0x221125C0UL,
+ 0x22113011UL, 0x2211F9EFUL, 0x22123051UL, 0x2212F9EFUL, 0x221329D1UL,
+ 0x22212541UL, 0x22213011UL, 0x2221F9EFUL, 0x22223451UL, 0x2222F9EFUL,
+ 0x22232551UL, 0x2223F56FUL, 0x22312551UL, 0x223229D1UL, 0x2232F56FUL,
+ 0x22332551UL, 0x2233F56FUL, 0x22875555UL, 0x22DAB5D7UL, 0x23112BD1UL,
+ 0x23115467UL, 0x231225D1UL, 0x232129D1UL, 0x232229D1UL, 0x2322F9EFUL,
+ 0x23233451UL, 0x2323F9EFUL, 0x23312551UL, 0x233229D1UL, 0x2332F9EFUL,
+ 0x2333F56FUL, 0x237FF557UL, 0x238569D5UL, 0x23D955D7UL, 0x24100BE7UL,
+ 0x248789E7UL, 0x24E315D7UL, 0x24FFFBEFUL, 0x259869E7UL, 0x25DFF5EFUL,
+ 0x25FFFBEFUL, 0x268789E7UL, 0x26DFA5D7UL, 0x26FFFBEFUL, 0x279649E7UL,
+ 0x27E425D7UL, 0x27FFFBEFUL, 0x288879E7UL, 0x28EFF5EFUL, 0x28FFFBEFUL,
+ 0x298439E7UL, 0x29F115EFUL, 0x29FFFBEFUL, 0x2A7659E7UL, 0x2AEF75D7UL,
+ 0x2AFFFBEFUL, 0x2B7C89E7UL, 0x2BEF95D7UL, 0x2BFFFBEFUL, 0x2C6659E7UL,
+ 0x2CD555D7UL, 0x2CFFFBEFUL, 0x2D6329E7UL, 0x2DDD55E7UL, 0x2DFFFBEBUL,
+ 0x2E8479D7UL, 0x2EEE35E7UL, 0x2EFFFBEFUL, 0x2F5459E7UL, 0x2FCF85D7UL,
+ 0x2FFEFBEBUL, 0x2FFFA5EFUL, 0x2FFFEBEFUL, 0x30001AE7UL, 0x30002001UL,
+ 0x311129C0UL, 0x31221015UL, 0x31232000UL, 0x31332451UL, 0x32112540UL,
+ 0x32131027UL, 0x32212440UL, 0x33452455UL, 0x4000F9D7UL, 0x4122F9D7UL,
+ 0x43F65555UL, 0x43FFF5D7UL, 0x44F55567UL, 0x44FFF5D7UL, 0x45F00557UL,
+ 0x45FFF5D7UL, 0x46F659D7UL, 0x471005E7UL, 0x47F449E7UL, 0x481005E7UL,
+ 0x48EFA9D5UL, 0x48FFF5EFUL, 0x49F449E7UL, 0x49FFF5EFUL, 0x4AEA79E7UL,
+ 0x4AFFF5EFUL, 0x4BE9C9D5UL, 0x4BFFF5EFUL, 0x4CE549E7UL, 0x4CFFF5EFUL,
+ 0x4DE359E7UL, 0x4DFFF5D7UL, 0x4EE469E7UL, 0x4EFFF5D7UL, 0x4FEF39E7UL,
+ 0x4FFFF5EFUL, 0x6000F9E7UL, 0x69FFF557UL, 0x6FFFF9D7UL, 0x811009D7UL,
+ 0x8EFFF555UL, 0xFFFFF9E7UL};
+
+static const SHORT ari_pk[64][17] = {
+ {708, 706, 579, 569, 568, 567, 479, 469, 297, 138, 97, 91, 72, 52, 38, 34,
+ 0},
+ {7619, 6917, 6519, 6412, 5514, 5003, 4683, 4563, 3907, 3297, 3125, 3060,
+ 2904, 2718, 2631, 2590, 0},
+ {7263, 4888, 4810, 4803, 1889, 415, 335, 327, 195, 72, 52, 49, 36, 20, 15,
+ 14, 0},
+ {3626, 2197, 2188, 2187, 582, 57, 47, 46, 30, 12, 9, 8, 6, 4, 3, 2, 0},
+ {7806, 5541, 5451, 5441, 2720, 834, 691, 674, 487, 243, 179, 167, 139, 98,
+ 77, 70, 0},
+ {6684, 4101, 4058, 4055, 1748, 426, 368, 364, 322, 257, 235, 232, 228, 222,
+ 217, 215, 0},
+ {9162, 5964, 5831, 5819, 3269, 866, 658, 638, 535, 348, 258, 244, 234, 214,
+ 195, 186, 0},
+ {10638, 8491, 8365, 8351, 4418, 2067, 1859, 1834, 1190, 601, 495, 478, 356,
+ 217, 174, 164, 0},
+ {13389, 10514, 10032, 9961, 7166, 3488, 2655, 2524, 2015, 1140, 760, 672,
+ 585, 426, 325, 283, 0},
+ {14861, 12788, 12115, 11952, 9987, 6657, 5323, 4984, 4324, 3001, 2205, 1943,
+ 1764, 1394, 1115, 978, 0},
+ {12876, 10004, 9661, 9610, 7107, 3435, 2711, 2595, 2257, 1508, 1059, 952,
+ 893, 753, 609, 538, 0},
+ {15125, 13591, 13049, 12874, 11192, 8543, 7406, 7023, 6291, 4922, 4104,
+ 3769, 3465, 2890, 2486, 2275, 0},
+ {14574, 13106, 12731, 12638, 10453, 7947, 7233, 7037, 6031, 4618, 4081,
+ 3906, 3465, 2802, 2476, 2349, 0},
+ {15070, 13179, 12517, 12351, 10742, 7657, 6200, 5825, 5264, 3998, 3014,
+ 2662, 2510, 2153, 1799, 1564, 0},
+ {15542, 14466, 14007, 13844, 12489, 10409, 9481, 9132, 8305, 6940, 6193,
+ 5867, 5458, 4743, 4291, 4047, 0},
+ {15165, 14384, 14084, 13934, 12911, 11485, 10844, 10513, 10002, 8993, 8380,
+ 8051, 7711, 7036, 6514, 6233, 0},
+ {15642, 14279, 13625, 13393, 12348, 9971, 8405, 7858, 7335, 6119, 4918,
+ 4376, 4185, 3719, 3231, 2860, 0},
+ {13408, 13407, 11471, 11218, 11217, 11216, 9473, 9216, 6480, 3689, 2857,
+ 2690, 2256, 1732, 1405, 1302, 0},
+ {16098, 15584, 15191, 14931, 14514, 13578, 12703, 12103, 11830, 11172,
+ 10475, 9867, 9695, 9281, 8825, 8389, 0},
+ {15844, 14873, 14277, 13996, 13230, 11535, 10205, 9543, 9107, 8086, 7085,
+ 6419, 6214, 5713, 5195, 4731, 0},
+ {16131, 15720, 15443, 15276, 14848, 13971, 13314, 12910, 12591, 11874,
+ 11225, 10788, 10573, 10077, 9585, 9209, 0},
+ {16331, 16330, 12283, 11435, 11434, 11433, 8725, 8049, 6065, 4138, 3187,
+ 2842, 2529, 2171, 1907, 1745, 0},
+ {16011, 15292, 14782, 14528, 14008, 12767, 11556, 10921, 10591, 9759, 8813,
+ 8043, 7855, 7383, 6863, 6282, 0},
+ {16380, 16379, 15159, 14610, 14609, 14608, 12859, 12111, 11046, 9536, 8348,
+ 7713, 7216, 6533, 5964, 5546, 0},
+ {16367, 16333, 16294, 16253, 16222, 16143, 16048, 15947, 15915, 15832,
+ 15731, 15619, 15589, 15512, 15416, 15310, 0},
+ {15967, 15319, 14937, 14753, 14010, 12638, 11787, 11360, 10805, 9706, 8934,
+ 8515, 8166, 7456, 6911, 6575, 0},
+ {4906, 3005, 2985, 2984, 875, 102, 83, 81, 47, 17, 12, 11, 8, 5, 4, 3, 0},
+ {7217, 4346, 4269, 4264, 1924, 428, 340, 332, 280, 203, 179, 175, 171, 164,
+ 159, 157, 0},
+ {16010, 15415, 15032, 14805, 14228, 13043, 12168, 11634, 11265, 10419, 9645,
+ 9110, 8892, 8378, 7850, 7437, 0},
+ {8573, 5218, 5046, 5032, 2787, 771, 555, 533, 443, 286, 218, 205, 197, 181,
+ 168, 162, 0},
+ {11474, 8095, 7822, 7796, 4632, 1443, 1046, 1004, 748, 351, 218, 194, 167,
+ 121, 93, 83, 0},
+ {16152, 15764, 15463, 15264, 14925, 14189, 13536, 13070, 12846, 12314,
+ 11763, 11277, 11131, 10777, 10383, 10011, 0},
+ {14187, 11654, 11043, 10919, 8498, 4885, 3778, 3552, 2947, 1835, 1283, 1134,
+ 998, 749, 585, 514, 0},
+ {14162, 11527, 10759, 10557, 8601, 5417, 4105, 3753, 3286, 2353, 1708, 1473,
+ 1370, 1148, 959, 840, 0},
+ {16205, 15902, 15669, 15498, 15213, 14601, 14068, 13674, 13463, 12970,
+ 12471, 12061, 11916, 11564, 11183, 10841, 0},
+ {15043, 12972, 12092, 11792, 10265, 7446, 5934, 5379, 4883, 3825, 3036,
+ 2647, 2507, 2185, 1901, 1699, 0},
+ {15320, 13694, 12782, 12352, 11191, 8936, 7433, 6671, 6255, 5366, 4622,
+ 4158, 4020, 3712, 3420, 3198, 0},
+ {16255, 16020, 15768, 15600, 15416, 14963, 14440, 14006, 13875, 13534,
+ 13137, 12697, 12602, 12364, 12084, 11781, 0},
+ {15627, 14503, 13906, 13622, 12557, 10527, 9269, 8661, 8117, 6933, 5994,
+ 5474, 5222, 4664, 4166, 3841, 0},
+ {16366, 16365, 14547, 14160, 14159, 14158, 11969, 11473, 8735, 6147, 4911,
+ 4530, 3865, 3180, 2710, 2473, 0},
+ {16257, 16038, 15871, 15754, 15536, 15071, 14673, 14390, 14230, 13842,
+ 13452, 13136, 13021, 12745, 12434, 12154, 0},
+ {15855, 14971, 14338, 13939, 13239, 11782, 10585, 9805, 9444, 8623, 7846,
+ 7254, 7079, 6673, 6262, 5923, 0},
+ {9492, 6318, 6197, 6189, 3004, 652, 489, 477, 333, 143, 96, 90, 78, 60, 50,
+ 47, 0},
+ {16313, 16191, 16063, 15968, 15851, 15590, 15303, 15082, 14968, 14704,
+ 14427, 14177, 14095, 13899, 13674, 13457, 0},
+ {8485, 5473, 5389, 5383, 2411, 494, 386, 377, 278, 150, 117, 112, 103, 89,
+ 81, 78, 0},
+ {10497, 7154, 6959, 6943, 3788, 1004, 734, 709, 517, 238, 152, 138, 120, 90,
+ 72, 66, 0},
+ {16317, 16226, 16127, 16040, 15955, 15762, 15547, 15345, 15277, 15111,
+ 14922, 14723, 14671, 14546, 14396, 14239, 0},
+ {16382, 16381, 15858, 15540, 15539, 15538, 14704, 14168, 13768, 13092,
+ 12452, 11925, 11683, 11268, 10841, 10460, 0},
+ {5974, 3798, 3758, 3755, 1275, 205, 166, 162, 95, 35, 26, 24, 18, 11, 8, 7,
+ 0},
+ {3532, 2258, 2246, 2244, 731, 135, 118, 115, 87, 45, 36, 34, 29, 21, 17, 16,
+ 0},
+ {7466, 4882, 4821, 4811, 2476, 886, 788, 771, 688, 531, 469, 457, 437, 400,
+ 369, 361, 0},
+ {9580, 5772, 5291, 5216, 3444, 1496, 1025, 928, 806, 578, 433, 384, 366,
+ 331, 296, 273, 0},
+ {10692, 7730, 7543, 7521, 4679, 1746, 1391, 1346, 1128, 692, 495, 458, 424,
+ 353, 291, 268, 0},
+ {11040, 7132, 6549, 6452, 4377, 1875, 1253, 1130, 958, 631, 431, 370, 346,
+ 296, 253, 227, 0},
+ {12687, 9332, 8701, 8585, 6266, 3093, 2182, 2004, 1683, 1072, 712, 608, 559,
+ 458, 373, 323, 0},
+ {13429, 9853, 8860, 8584, 6806, 4039, 2862, 2478, 2239, 1764, 1409, 1224,
+ 1178, 1077, 979, 903, 0},
+ {14685, 12163, 11061, 10668, 9101, 6345, 4871, 4263, 3908, 3200, 2668, 2368,
+ 2285, 2106, 1942, 1819, 0},
+ {13295, 11302, 10999, 10945, 7947, 5036, 4490, 4385, 3391, 2185, 1836, 1757,
+ 1424, 998, 833, 785, 0},
+ {4992, 2993, 2972, 2970, 1269, 575, 552, 549, 530, 505, 497, 495, 493, 489,
+ 486, 485, 0},
+ {15419, 13862, 13104, 12819, 11429, 8753, 7220, 6651, 6020, 4667, 3663,
+ 3220, 2995, 2511, 2107, 1871, 0},
+ {12468, 9263, 8912, 8873, 5758, 2193, 1625, 1556, 1187, 589, 371, 330, 283,
+ 200, 149, 131, 0},
+ {15870, 15076, 14615, 14369, 13586, 12034, 10990, 10423, 9953, 8908, 8031,
+ 7488, 7233, 6648, 6101, 5712, 0},
+ {1693, 978, 976, 975, 194, 18, 16, 15, 11, 7, 6, 5, 4, 3, 2, 1, 0},
+ {7992, 5218, 5147, 5143, 2152, 366, 282, 276, 173, 59, 38, 35, 27, 16, 11,
+ 10, 0}};
+
+typedef struct {
+ int low;
+ int high;
+ int vobf;
+} Tastat;
+
+static inline INT mul_sbc_14bits(INT r, INT c) {
+ return (((INT)r) * ((INT)c)) >> stat_bitsnew;
+}
+
+static inline INT ari_decode_14bits(HANDLE_FDK_BITSTREAM hBs, Tastat *s,
+ const SHORT *RESTRICT c_freq, int cfl) {
+ INT symbol;
+ INT low, high, range, value;
+ INT c;
+ const SHORT *p;
+
+ low = s->low;
+ high = s->high;
+ value = s->vobf;
+
+ range = high - low + 1;
+ c = (((int)(value - low + 1)) << stat_bitsnew) - ((int)1);
+ p = (const SHORT *)(c_freq - 1);
+
+ if (cfl == (VAL_ESC + 1)) {
+ /* In 50% of all cases, the first entry is the right one, so we check it
+ * prior to all others */
+ if ((p[1] * range) > c) {
+ p += 1;
+ if ((p[8] * range) > c) {
+ p += 8;
+ }
+ if ((p[4] * range) > c) {
+ p += 4;
+ }
+ if ((p[2] * range) > c) {
+ p += 2;
+ }
+ if ((p[1] * range) > c) {
+ p += 1;
+ }
+ }
+ } else if (cfl == 4) {
+ if ((p[2] * range) > c) {
+ p += 2;
+ }
+ if ((p[1] * range) > c) {
+ p += 1;
+ }
+ } else if (cfl == 2) {
+ if ((p[1] * range) > c) {
+ p += 1;
+ }
+ } else if (cfl == 27) {
+ const SHORT *p_24 = p + 24;
+
+ if ((p[16] * range) > c) {
+ p += 16;
+ }
+ if ((p[8] * range) > c) {
+ p += 8;
+ }
+ if (p != p_24) {
+ if ((p[4] * range) > c) {
+ p += 4;
+ }
+ }
+ if ((p[2] * range) > c) {
+ p += 2;
+ }
+
+ if (p != &p_24[2]) {
+ if ((p[1] * range) > c) {
+ p += 1;
+ }
+ }
+ }
+
+ symbol = (INT)(p - (const SHORT *)(c_freq - 1));
+
+ if (symbol) {
+ high = low + mul_sbc_14bits(range, c_freq[symbol - 1]) - 1;
+ }
+
+ low += mul_sbc_14bits(range, c_freq[symbol]);
+
+ USHORT us_high = (USHORT)high;
+ USHORT us_low = (USHORT)low;
+ while (1) {
+ if (us_high & 0x8000) {
+ if (!(us_low & 0x8000)) {
+ if (us_low & 0x4000 && !(us_high & 0x4000)) {
+ us_low -= 0x4000;
+ us_high -= 0x4000;
+ value -= 0x4000;
+ } else
+ break;
+ }
+ }
+ us_low = us_low << 1;
+ us_high = (us_high << 1) | 1;
+ value = (value << 1) | FDKreadBit(hBs);
+ }
+ s->low = (int)us_low;
+ s->high = (int)us_high;
+ s->vobf = value & 0xFFFF;
+
+ return symbol;
+}
+
+static inline void copyTableAmrwbArith2(UCHAR tab[], int sizeIn, int sizeOut) {
+ int i;
+ int j;
+ int k = 2;
+
+ tab += 2;
+
+ if (sizeIn < sizeOut) {
+ tab[sizeOut + 0] = tab[sizeIn + 0];
+ tab[sizeOut + 1] = tab[sizeIn + 1];
+ if (sizeIn < (sizeOut >> 2)) {
+ k = 8;
+ } else if (sizeIn == (sizeOut >> 2)) {
+ k = 4;
+ }
+
+ i = sizeOut - 1;
+ j = sizeIn - 1;
+
+ for (; i >= 0; j--) {
+ UCHAR tq_data0 = tab[j];
+
+ for (int l = (k >> 1); l > 0; l--) {
+ tab[i--] = tq_data0;
+ tab[i--] = tq_data0;
+ }
+ }
+ } else {
+ if (sizeOut < (sizeIn >> 2)) {
+ k = 8;
+ } else if (sizeOut == (sizeIn >> 2)) {
+ k = 4;
+ }
+
+ for (i = 0, j = 0; i < sizeOut; j += k) {
+ UCHAR tq_data0 = tab[j];
+
+ tab[i++] = tq_data0;
+ }
+ tab[sizeOut + 0] = tab[sizeIn + 0];
+ tab[sizeOut + 1] = tab[sizeIn + 1];
+ }
+}
+
+static inline ULONG get_pk_v2(ULONG s) {
+ const ULONG *p = ari_merged_hash_ps;
+ ULONG s12 = (fMax((UINT)s, (UINT)1) << 12) - 1;
+ if (s12 > p[485]) {
+ p += 486; /* 742 - 256 = 486 */
+ } else {
+ if (s12 > p[255]) p += 256;
+ }
+
+ if (s12 > p[127]) {
+ p += 128;
+ }
+ if (s12 > p[63]) {
+ p += 64;
+ }
+ if (s12 > p[31]) {
+ p += 32;
+ }
+ if (s12 > p[15]) {
+ p += 16;
+ }
+ if (s12 > p[7]) {
+ p += 8;
+ }
+ if (s12 > p[3]) {
+ p += 4;
+ }
+ if (s12 > p[1]) {
+ p += 2;
+ }
+ ULONG j = p[0];
+ if (s12 > j) j = p[1];
+ if (s != (j >> 12)) j >>= 6;
+ return (j & 0x3F);
+}
+
+static void decode2(HANDLE_FDK_BITSTREAM bbuf, UCHAR *RESTRICT c_prev,
+ FIXP_DBL *RESTRICT pSpectralCoefficient, INT n, INT nt) {
+ Tastat as;
+ int i, l, r;
+ INT lev, esc_nb, pki;
+ USHORT state_inc;
+ UINT s;
+
+ int c_3 = 0; /* context of current frame 3 time steps ago */
+ int c_2 = 0; /* context of current frame 2 time steps ago */
+ int c_1 = 0; /* context of current frame 1 time steps ago */
+ int c_0 = 1; /* context of current frame to be calculated */
+
+ /* ari_start_decoding_14bits */
+ as.low = 0;
+ as.high = ari_q4new;
+ as.vobf = FDKreadBits(bbuf, cbitsnew);
+
+ /* arith_map_context */
+ state_inc = c_prev[0] << 12;
+
+ for (i = 0; i < n; i++) {
+ /* arith_get_context */
+ s = state_inc >> 8;
+ s = s + (c_prev[i + 1] << 8);
+ s = (s << 4) + c_1;
+
+ state_inc = s;
+
+ if (i > 3) {
+ /* Cumulative amplitude below 2 */
+ if ((c_1 + c_2 + c_3) < 5) {
+ s += 0x10000;
+ }
+ }
+
+ /* MSBs decoding */
+ for (lev = esc_nb = 0;;) {
+ pki = get_pk_v2(s + (esc_nb << (VAL_ESC + 1)));
+ r = ari_decode_14bits(bbuf, &as, ari_pk[pki], VAL_ESC + 1);
+ if (r < VAL_ESC) {
+ break;
+ }
+
+ lev++;
+
+ if (esc_nb < 7) {
+ esc_nb++;
+ }
+ }
+
+ /* Stop symbol */
+ if (r == 0) {
+ if (esc_nb > 0) {
+ break; /* Stop symbol */
+ }
+ c_0 = 1;
+ } else /* if (r==0) */
+ {
+ INT b = r >> 2;
+ INT a = r & 0x3;
+
+ /* LSBs decoding */
+ for (l = 0; l < lev; l++) {
+ {
+ int pidx = (a == 0) ? 1 : ((b == 0) ? 0 : 2);
+ r = ari_decode_14bits(bbuf, &as, ari_lsb2[pidx], 4);
+ }
+ a = (a << 1) | (r & 1);
+ b = (b << 1) | (r >> 1);
+ }
+
+ pSpectralCoefficient[2 * i] = (FIXP_DBL)a;
+ pSpectralCoefficient[2 * i + 1] = (FIXP_DBL)b;
+
+ c_0 = a + b + 1;
+ if (c_0 > 0xF) {
+ c_0 = 0xF;
+ }
+
+ } /* endif (r==0) */
+
+ /* arith_update_context */
+ c_3 = c_2;
+ c_2 = c_1;
+ c_1 = c_0;
+ c_prev[i] = (UCHAR)c_0;
+
+ } /* for (i=0; i<n; i++) */
+
+ FDKpushBack(bbuf, cbitsnew - 2);
+
+ /* We need to run only from 0 to i-1 since all other q[i][1].a,b will be
+ * cleared later */
+ int j = i;
+ for (i = 0; i < j; i++) {
+ int bits = 0;
+ if (pSpectralCoefficient[2 * i] != (FIXP_DBL)0) bits++;
+ if (pSpectralCoefficient[2 * i + 1] != (FIXP_DBL)0) bits++;
+
+ if (bits) {
+ r = FDKreadBits(bbuf, bits);
+ if (pSpectralCoefficient[2 * i] != (FIXP_DBL)0 && !(r >> (bits - 1))) {
+ pSpectralCoefficient[2 * i] = -pSpectralCoefficient[2 * i];
+ }
+ if (pSpectralCoefficient[2 * i + 1] != (FIXP_DBL)0 && !(r & 1)) {
+ pSpectralCoefficient[2 * i + 1] = -pSpectralCoefficient[2 * i + 1];
+ }
+ }
+ }
+
+ FDKmemset(&c_prev[i], 1, sizeof(c_prev[0]) * (nt - i));
+}
+
+CArcoData *CArco_Create(void) { return GetArcoData(); }
+
+void CArco_Destroy(CArcoData *pArcoData) { FreeArcoData(&pArcoData); }
+
+ARITH_CODING_ERROR CArco_DecodeArithData(CArcoData *pArcoData,
+ HANDLE_FDK_BITSTREAM hBs,
+ FIXP_DBL *RESTRICT mdctSpectrum,
+ int lg, int lg_max,
+ int arith_reset_flag) {
+ ARITH_CODING_ERROR ErrorStatus = ARITH_CODER_OK;
+
+ /* Check lg and lg_max consistency. */
+ if (lg_max < lg) {
+ return ARITH_CODER_ERROR;
+ }
+
+ FDKmemclear(mdctSpectrum, lg_max * sizeof(FIXP_DBL));
+
+ /* arith_map_context */
+ if (arith_reset_flag) {
+ FDKmemclear(pArcoData->c_prev,
+ sizeof(pArcoData->c_prev[0]) * ((lg_max / 2) + 4));
+ } else {
+ if (lg_max != pArcoData->m_numberLinesPrev) {
+ if (pArcoData->m_numberLinesPrev == 0) {
+ /* Cannot decode without a valid AC context */
+ return ARITH_CODER_ERROR;
+ }
+
+ /* short-to-long or long-to-short block transition */
+ /* Current length differs compared to previous - perform up/downmix of
+ * m_qbuf */
+ copyTableAmrwbArith2(pArcoData->c_prev, pArcoData->m_numberLinesPrev >> 1,
+ lg_max >> 1);
+ }
+ }
+
+ pArcoData->m_numberLinesPrev = lg_max;
+
+ if (lg > 0) {
+ decode2(hBs, pArcoData->c_prev + 2, mdctSpectrum, lg >> 1, lg_max >> 1);
+ } else {
+ FDKmemset(&pArcoData->c_prev[2], 1,
+ sizeof(pArcoData->c_prev[2]) * (lg_max >> 1));
+ }
+
+ if ((INT)FDKgetValidBits(hBs) < 0) {
+ return ARITH_CODER_ERROR;
+ }
+
+ return ErrorStatus;
+}