Commit 4e62ce87 authored by Fabian Bauer's avatar Fabian Bauer Committed by Manuel Jander
Browse files

Revert "introduce FIX_1326_SUBSTITUTE_DIV_SQRT_IOSQRT_2"

This reverts commit 18d10e60.
parent 01fbdaf8
Loading
Loading
Loading
Loading
+1 −24
Original line number Diff line number Diff line
@@ -29,7 +29,7 @@
    the United Nations Convention on Contracts on the International Sales of Goods.

*******************************************************************************************************/
#define FIX_1326_SUBSTITUTE_DIV_SQRT_IOSQRT_2

#include <stdint.h>
#include "options.h"
#include <assert.h>
@@ -4232,28 +4232,6 @@ static void formulate2x2MixingMatrix_fx(

    eig2x2_fx( tmpRe_fx[0][0], tmpRe_fx[1][1], q_temp, tmpRe_fx[1][0], tmpIm_fx[1][0], q_temp, Ure_fx, Uim_fx, &q_U, D_fx, &q_D );


#ifdef FIX_1326_SUBSTITUTE_DIV_SQRT_IOSQRT_2
    IF( D_fx[0] == 0 )
    {
#ifndef FIX_1009_REPLACE_DIV_SQRT_BY_ISQRT_LC
        temp = BASOP_Util_Divide3232_Scale_cadence( ONE_IN_Q30, 4611686, &exp ); // 4611686 = 1e-12 in Q62
        exp = sub( exp, sub( Q30, 62 ) );
#else
        temp = ONE_DIV_EPSILON_MANT; /* Result of 1.0/eps with full precision */
        move32();
        exp = ONE_DIV_EPSILON_EXP;
        move16();
#endif
        div_fx[0] = Sqrt32( temp, &exp ); // Q = 31 - exp
        move32();
    }
    ELSE
    {
        exp = sub( 31, q_D );
        div_fx[0] = ISqrt32( D_fx[0], &exp );
    }
#else
    IF( D_fx[0] == 0 )
    {
        temp = ONE_DIV_EPSILON_MANT; /* Result of 1.0/eps with full precision */
@@ -4269,7 +4247,6 @@ static void formulate2x2MixingMatrix_fx(
    }
    div_fx[0] = Sqrt32( temp, &exp ); // Q = 31 - exp
    move32();
#endif /*FIX_1326_SUBSTITUTE_DIV_SQRT_IOSQRT_2*/

#ifdef FIX_1326_SUBSTITUTE_DIV_SQRT_IOSQRT
    // Sqrt(1)