| Index: third_party/openmax_dl/dl/sp/src/armSP_FFT_CToC_FC32_Radix2_unsafe_s.S
|
| diff --git a/third_party/openmax_dl/dl/sp/src/armSP_FFT_CToC_FC32_Radix2_unsafe_s.S b/third_party/openmax_dl/dl/sp/src/armSP_FFT_CToC_FC32_Radix2_unsafe_s.S
|
| new file mode 100644
|
| index 0000000000000000000000000000000000000000..73c1f4b82f3d7ffbb4e9f6e46628f576d92c0fff
|
| --- /dev/null
|
| +++ b/third_party/openmax_dl/dl/sp/src/armSP_FFT_CToC_FC32_Radix2_unsafe_s.S
|
| @@ -0,0 +1,191 @@
|
| +@//
|
| +@// Copyright (c) 2013 The WebRTC project authors. All Rights Reserved.
|
| +@//
|
| +@// Use of this source code is governed by a BSD-style license
|
| +@// that can be found in the LICENSE file in the root of the source
|
| +@// tree. An additional intellectual property rights grant can be found
|
| +@// in the file PATENTS. All contributing project authors may
|
| +@// be found in the AUTHORS file in the root of the source tree.
|
| +@//
|
| +@// This is a modification of armSP_FFT_CToC_SC32_Radix2_unsafe_s.s
|
| +@// to support float instead of SC32.
|
| +@//
|
| +
|
| +@// Description:
|
| +@// Compute a Radix 2 DIT in-order out-of-place FFT stage for an N point
|
| +@// complex signal. This handles the general stage, not the first or last
|
| +@// stage.
|
| +@//
|
| +@//
|
| +
|
| +
|
| +@// Include standard headers
|
| +
|
| +#include "dl/api/armCOMM_s.h"
|
| +#include "dl/api/omxtypes_s.h"
|
| +
|
| +
|
| +@// Import symbols required from other files
|
| +@// (For example tables)
|
| +
|
| +
|
| +
|
| +@// Set debugging level
|
| +@//DEBUG_ON SETL {TRUE}
|
| +
|
| +
|
| +
|
| +@// Guarding implementation by the processor name
|
| +
|
| +
|
| +
|
| +
|
| +@// Guarding implementation by the processor name
|
| +
|
| +
|
| +@//Input Registers
|
| +
|
| +#define pSrc r0
|
| +#define pDst r2
|
| +#define pTwiddle r1
|
| +#define subFFTNum r6
|
| +#define subFFTSize r7
|
| +
|
| +
|
| +@//Output Registers
|
| +
|
| +
|
| +@//Local Scratch Registers
|
| +
|
| +#define outPointStep r3
|
| +#define pointStep r4
|
| +#define grpCount r5
|
| +#define setCount r8
|
| +@//const RN 9
|
| +#define step r10
|
| +#define dstStep r11
|
| +#define pTable r9
|
| +#define pTmp r9
|
| +
|
| +@// Neon Registers
|
| +
|
| +#define dW D0.F32
|
| +#define dX0 D2.F32
|
| +#define dX1 D3.F32
|
| +#define dX2 D4.F32
|
| +#define dX3 D5.F32
|
| +#define dY0 D6.F32
|
| +#define dY1 D7.F32
|
| +#define dY2 D8.F32
|
| +#define dY3 D9.F32
|
| +#define qT0 D10.F32
|
| +#define qT1 D11.F32
|
| +
|
| +
|
| + .MACRO FFTSTAGE scaled, inverse, name
|
| +
|
| + @// Define stack arguments
|
| +
|
| +
|
| + @// Update grpCount and grpSize rightaway inorder to reuse pGrpCount
|
| + @// and pGrpSize regs
|
| +
|
| + LSR subFFTNum,subFFTNum,#1 @//grpSize
|
| + LSL grpCount,subFFTSize,#1
|
| +
|
| +
|
| + @// pT0+1 increments pT0 by 8 bytes
|
| + @// pT0+pointStep = increment of 8*pointStep bytes = 4*grpSize bytes
|
| + MOV pointStep,subFFTNum,LSL #2
|
| +
|
| + @// update subFFTSize for the next stage
|
| + MOV subFFTSize,grpCount
|
| +
|
| + @// pOut0+1 increments pOut0 by 8 bytes
|
| + @// pOut0+outPointStep == increment of 8*outPointStep bytes =
|
| + @// 4*size bytes
|
| + SMULBB outPointStep,grpCount,pointStep
|
| + LSL pointStep,pointStep,#1
|
| +
|
| +
|
| + RSB step,pointStep,#16
|
| + RSB dstStep,outPointStep,#16
|
| +
|
| + @// Loop on the groups
|
| +
|
| +radix2GrpLoop\name :
|
| + MOV setCount,pointStep,LSR #3
|
| + VLD1 dW,[pTwiddle],pointStep @//[wi | wr]
|
| +
|
| +
|
| + @// Loop on the sets
|
| +
|
| +
|
| +radix2SetLoop\name :
|
| +
|
| +
|
| + @// point0: dX0-real part dX1-img part
|
| + VLD2 {dX0,dX1},[pSrc],pointStep
|
| + @// point1: dX2-real part dX3-img part
|
| + VLD2 {dX2,dX3},[pSrc],step
|
| +
|
| + SUBS setCount,setCount,#2
|
| +
|
| + .ifeqs "\inverse", "TRUE"
|
| + VMUL qT0,dX2,dW[0]
|
| + VMLA qT0,dX3,dW[1] @// real part
|
| + VMUL qT1,dX3,dW[0]
|
| + VMLS qT1,dX2,dW[1] @// imag part
|
| +
|
| + .else
|
| +
|
| + VMUL qT0,dX2,dW[0]
|
| + VMLS qT0,dX3,dW[1] @// real part
|
| + VMUL qT1,dX3,dW[0]
|
| + VMLA qT1,dX2,dW[1] @// imag part
|
| +
|
| + .endif
|
| +
|
| + VSUB dY0,dX0,qT0
|
| + VSUB dY1,dX1,qT1
|
| + VADD dY2,dX0,qT0
|
| + VADD dY3,dX1,qT1
|
| +
|
| + VST2 {dY0,dY1},[pDst],outPointStep
|
| + @// dstStep = -outPointStep + 16
|
| + VST2 {dY2,dY3},[pDst],dstStep
|
| +
|
| + BGT radix2SetLoop\name
|
| +
|
| + SUBS grpCount,grpCount,#2
|
| + ADD pSrc,pSrc,pointStep
|
| + BGT radix2GrpLoop\name
|
| +
|
| +
|
| + @// Reset and Swap pSrc and pDst for the next stage
|
| + MOV pTmp,pDst
|
| + @// pDst -= 4*size; pSrc -= 8*size bytes
|
| + SUB pDst,pSrc,outPointStep,LSL #1
|
| + SUB pSrc,pTmp,outPointStep
|
| +
|
| + @// Reset pTwiddle for the next stage
|
| + @// pTwiddle -= 4*size bytes
|
| + SUB pTwiddle,pTwiddle,outPointStep
|
| +
|
| +
|
| + .endm
|
| +
|
| +
|
| +
|
| + M_START armSP_FFTFwd_CToC_FC32_Radix2_OutOfPlace_unsafe,r4
|
| + FFTSTAGE "FALSE","FALSE",FWD
|
| + M_END
|
| +
|
| +
|
| +
|
| + M_START armSP_FFTInv_CToC_FC32_Radix2_OutOfPlace_unsafe,r4
|
| + FFTSTAGE "FALSE","TRUE",INV
|
| + M_END
|
| +
|
| +
|
| + .end
|
|
|