// File: /workspace/github_cpp_makefile_star2_O1/kebwi[P]Trifid/src/FFT/fftlib.cpp
// Function: ffts1(double*, long, long, double*, short*)
void ffts1(Float_t *ioptr, long M, long Rows, Float_t *Utbl, short *BRLow){
/* Compute in-place complex fft on the rows of the input array	*/
/* INPUTS */
/* *ioptr = input data array	*/
/* M = log2 of fft size	(ex M=10 for 1024 point fft) */
/* Rows = number of rows in ioptr array (use Rows of 1 if ioptr is a 1 dimensional array)	*/
/* *Utbl = cosine table	*/
/* *BRLow = bit reversed counter table	*/
/* OUTPUTS */
/* *ioptr = output data array	*/

long 	StageCnt;
long 	NDiffU;

switch (M){
case 0:
	break;
case 1:
	for (;Rows>0;Rows--){
		fft2pt(ioptr);				/* a 2 pt fft */
		ioptr += 2*POW2(M);
	}
	break;
case 2:
	for (;Rows>0;Rows--){
		fft4pt(ioptr);				/* a 4 pt fft */
		ioptr += 2*POW2(M);
	}
	break;
case 3:
	for (;Rows>0;Rows--){
		fft8pt(ioptr);				/* an 8 pt fft */
		ioptr += 2*POW2(M);
	}
	break;
default:
	for (;Rows>0;Rows--){

		bitrevR2(ioptr, M, BRLow);						/* bit reverse and first radix 2 stage */
		
		StageCnt = (M-1) / 3;		// number of radix 8 stages
		NDiffU = 2;				// one radix 2 stage already complete

		if ( (M-1-(StageCnt * 3)) == 1 ){
			bfR2(ioptr, M, NDiffU);				/* 1 radix 2 stage */
			NDiffU *= 2;
		}

		if ( (M-1-(StageCnt * 3)) == 2 ){
			bfR4(ioptr, M, NDiffU);			/* 1 radix 4 stage */
			NDiffU *= 4;
		}

		if (M <= MCACHE)
			bfstages(ioptr, M, Utbl, 1, NDiffU, StageCnt);		/*  RADIX 8 Stages	*/

		else{
			fftrecurs(ioptr, M, Utbl, 1, NDiffU, StageCnt);		/*  RADIX 8 Stages	*/
		}

		ioptr += 2*POW2(M);
	}
}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x38, %rsp
movq %rdx, %r14
movq %rsi, %r15
movq %rdi, %rbx
cmpq $0x3, %rsi
ja 0x25047
leaq 0xac9d(%rip), %rax      # 0x2f8d0
movslq (%rax,%r15,4), %rcx
addq %rax, %rcx
jmpq *%rcx
testq %r14, %r14
jle 0x25280
movl $0x2, %eax
movl %r15d, %ecx
shlq %cl, %rax
addq $0x18, %rbx
shlq $0x3, %rax
incq %r14
movsd -0x18(%rbx), %xmm0
movsd -0x10(%rbx), %xmm1
movsd -0x8(%rbx), %xmm2
movsd (%rbx), %xmm3
movapd %xmm0, %xmm4
addsd %xmm2, %xmm4
movapd %xmm1, %xmm5
addsd %xmm3, %xmm5
subsd %xmm2, %xmm0
subsd %xmm3, %xmm1
movsd %xmm4, -0x18(%rbx)
movsd %xmm5, -0x10(%rbx)
movsd %xmm0, -0x8(%rbx)
movsd %xmm1, (%rbx)
addq %rax, %rbx
decq %r14
cmpq $0x1, %r14
jg 0x24c5b
jmp 0x25280
testq %r14, %r14
jle 0x25280
movl $0x2, %eax
movl %r15d, %ecx
shlq %cl, %rax
addq $0x38, %rbx
shlq $0x3, %rax
incq %r14
movsd -0x38(%rbx), %xmm1
movsd -0x30(%rbx), %xmm0
movsd -0x18(%rbx), %xmm4
movsd -0x10(%rbx), %xmm5
movsd -0x28(%rbx), %xmm6
movsd -0x20(%rbx), %xmm7
movsd -0x8(%rbx), %xmm8
movsd (%rbx), %xmm9
movapd %xmm1, %xmm2
addsd %xmm4, %xmm2
movapd %xmm0, %xmm3
addsd %xmm5, %xmm3
subsd %xmm4, %xmm1
subsd %xmm5, %xmm0
movapd %xmm6, %xmm4
subsd %xmm8, %xmm4
movapd %xmm7, %xmm5
subsd %xmm9, %xmm5
addsd %xmm8, %xmm6
addsd %xmm9, %xmm7
movapd %xmm2, %xmm8
addsd %xmm6, %xmm8
movapd %xmm3, %xmm9
addsd %xmm7, %xmm9
subsd %xmm6, %xmm2
subsd %xmm7, %xmm3
movapd %xmm1, %xmm6
subsd %xmm5, %xmm6
movapd %xmm0, %xmm7
addsd %xmm4, %xmm7
addsd %xmm1, %xmm5
subsd %xmm4, %xmm0
movsd %xmm8, -0x38(%rbx)
movsd %xmm9, -0x30(%rbx)
movsd %xmm5, -0x28(%rbx)
movsd %xmm0, -0x20(%rbx)
movsd %xmm2, -0x18(%rbx)
movsd %xmm3, -0x10(%rbx)
movsd %xmm6, -0x8(%rbx)
movsd %xmm7, (%rbx)
addq %rax, %rbx
decq %r14
cmpq $0x1, %r14
jg 0x24cc9
jmp 0x25280
testq %r14, %r14
jle 0x25280
movl $0x2, %eax
movl %r15d, %ecx
shlq %cl, %rax
addq $0x78, %rbx
shlq $0x3, %rax
incq %r14
movsd -0x78(%rbx), %xmm6
movsd -0x70(%rbx), %xmm3
movsd -0x38(%rbx), %xmm0
movsd -0x30(%rbx), %xmm4
movsd -0x58(%rbx), %xmm12
movsd -0x50(%rbx), %xmm5
movsd -0x18(%rbx), %xmm1
movsd -0x10(%rbx), %xmm2
movsd -0x68(%rbx), %xmm7
movsd %xmm7, 0x18(%rsp)
movsd -0x60(%rbx), %xmm7
movsd -0x28(%rbx), %xmm8
movsd -0x20(%rbx), %xmm13
movsd -0x48(%rbx), %xmm14
movsd -0x40(%rbx), %xmm11
movsd -0x8(%rbx), %xmm9
movapd %xmm6, %xmm10
addsd %xmm0, %xmm10
subsd %xmm0, %xmm6
movapd %xmm3, %xmm15
addsd %xmm4, %xmm15
subsd %xmm4, %xmm3
movapd %xmm12, %xmm0
subsd %xmm1, %xmm0
addsd %xmm1, %xmm12
movapd %xmm5, %xmm4
subsd %xmm2, %xmm4
addsd %xmm2, %xmm5
movapd %xmm10, %xmm1
addsd %xmm12, %xmm1
movsd %xmm1, 0x10(%rsp)
subsd %xmm12, %xmm10
movapd %xmm15, %xmm1
addsd %xmm5, %xmm1
movsd %xmm1, 0x20(%rsp)
subsd %xmm5, %xmm15
movapd %xmm6, %xmm5
subsd %xmm4, %xmm5
addsd %xmm6, %xmm4
movapd %xmm3, %xmm6
addsd %xmm0, %xmm6
subsd %xmm0, %xmm3
movsd 0x18(%rsp), %xmm0
movapd %xmm0, %xmm12
addsd %xmm8, %xmm12
subsd %xmm8, %xmm0
movapd %xmm0, %xmm2
movapd %xmm7, %xmm8
addsd %xmm13, %xmm8
subsd %xmm13, %xmm7
movapd %xmm14, %xmm1
subsd %xmm9, %xmm1
addsd %xmm9, %xmm14
movsd (%rbx), %xmm0
movapd %xmm11, %xmm9
subsd %xmm0, %xmm9
addsd %xmm0, %xmm11
movapd %xmm12, %xmm13
addsd %xmm14, %xmm13
subsd %xmm14, %xmm12
movapd %xmm8, %xmm14
addsd %xmm11, %xmm14
subsd %xmm11, %xmm8
movapd %xmm2, %xmm11
subsd %xmm9, %xmm11
addsd %xmm2, %xmm9
movapd %xmm7, %xmm2
addsd %xmm1, %xmm2
subsd %xmm1, %xmm7
movsd 0x10(%rsp), %xmm0
movapd %xmm0, %xmm1
subsd %xmm13, %xmm1
movsd %xmm1, 0x18(%rsp)
addsd %xmm0, %xmm13
movsd 0x20(%rsp), %xmm0
movapd %xmm0, %xmm1
subsd %xmm14, %xmm1
movsd %xmm1, 0x10(%rsp)
addsd %xmm0, %xmm14
movapd %xmm10, %xmm1
subsd %xmm8, %xmm1
addsd %xmm10, %xmm8
movapd %xmm15, %xmm10
addsd %xmm12, %xmm10
subsd %xmm12, %xmm15
xorpd 0x752e(%rip), %xmm7     # 0x2c490
movsd 0xa936(%rip), %xmm0     # 0x2f8a0
mulsd %xmm0, %xmm9
movapd %xmm4, %xmm12
subsd %xmm9, %xmm12
mulsd %xmm0, %xmm7
addsd %xmm7, %xmm12
addsd %xmm3, %xmm9
addsd %xmm7, %xmm9
mulsd %xmm0, %xmm11
movapd %xmm11, %xmm7
addsd %xmm5, %xmm7
mulsd %xmm0, %xmm2
subsd %xmm2, %xmm7
addsd %xmm6, %xmm11
addsd %xmm2, %xmm11
movsd %xmm13, -0x78(%rbx)
movsd %xmm14, -0x70(%rbx)
addsd %xmm4, %xmm4
subsd %xmm12, %xmm4
movsd %xmm4, -0x68(%rbx)
addsd %xmm3, %xmm3
subsd %xmm9, %xmm3
movsd %xmm3, -0x60(%rbx)
movsd %xmm8, -0x58(%rbx)
movsd %xmm15, -0x50(%rbx)
addsd %xmm5, %xmm5
subsd %xmm7, %xmm5
movsd %xmm5, -0x48(%rbx)
addsd %xmm6, %xmm6
subsd %xmm11, %xmm6
movsd %xmm6, -0x40(%rbx)
movsd 0x18(%rsp), %xmm0
movsd %xmm0, -0x38(%rbx)
movsd 0x10(%rsp), %xmm0
movsd %xmm0, -0x30(%rbx)
movsd %xmm12, -0x28(%rbx)
movsd %xmm9, -0x20(%rbx)
movsd %xmm1, -0x18(%rbx)
movsd %xmm10, -0x10(%rbx)
movsd %xmm7, -0x8(%rbx)
movsd %xmm11, (%rbx)
addq %rax, %rbx
decq %r14
cmpq $0x1, %r14
jg 0x24db7
jmp 0x25280
movq %rcx, 0x18(%rsp)
testq %r14, %r14
jle 0x25280
leaq -0x1(%r15), %rax
movabsq $0x5555555555555556, %rcx # imm = 0x5555555555555556
imulq %rcx
movq %rdx, %rax
shrq $0x3f, %rax
addq %rdx, %rax
movq %rax, 0x10(%rsp)
leaq (%rax,%rax,2), %rax
negq %rax
leaq (%r15,%rax), %rcx
decq %rcx
xorl %eax, %eax
movq %rcx, 0x20(%rsp)
cmpq $0x1, %rcx
sete %al
setne %dl
cmpq $0x3, %r15
setb %sil
movl $0x1, %edi
movl %r15d, %ecx
shlq %cl, %rdi
movl $0x2, %r9d
shlq %cl, %r9
movq %r9, 0x30(%rsp)
shrq $0x3, %rdi
movq %rdi, 0x28(%rsp)
orb %dl, %sil
movb %sil, 0xf(%rsp)
leaq 0x2(,%rax,2), %rbp
movq %rbx, %rdi
movq %r15, %rsi
movq %r8, %r13
movq %r8, %rdx
callq 0x270a1
movq %rbp, %r12
cmpb $0x0, 0xf(%rsp)
jne 0x25212
movq 0x28(%rsp), %rax
movq %rbx, %rcx
movsd 0x20(%rcx), %xmm0
movsd (%rcx), %xmm1
movsd 0x8(%rcx), %xmm2
movsd 0x28(%rcx), %xmm3
movsd 0x40(%rcx), %xmm4
movsd 0x60(%rcx), %xmm5
movsd 0x48(%rcx), %xmm6
movsd 0x68(%rcx), %xmm7
movapd %xmm1, %xmm8
addsd %xmm0, %xmm8
movapd %xmm2, %xmm9
addsd %xmm3, %xmm9
subsd %xmm0, %xmm1
subsd %xmm3, %xmm2
movapd %xmm4, %xmm0
addsd %xmm5, %xmm0
movapd %xmm6, %xmm3
addsd %xmm7, %xmm3
subsd %xmm5, %xmm4
subsd %xmm7, %xmm6
movsd %xmm8, (%rcx)
movsd %xmm9, 0x8(%rcx)
movsd %xmm1, 0x20(%rcx)
movsd %xmm2, 0x28(%rcx)
movsd %xmm0, 0x40(%rcx)
movsd %xmm3, 0x48(%rcx)
movsd %xmm4, 0x60(%rcx)
movsd %xmm6, 0x68(%rcx)
movsd 0x10(%rcx), %xmm0
movsd 0x38(%rcx), %xmm1
movsd 0x18(%rcx), %xmm2
movsd 0x30(%rcx), %xmm3
movsd 0x50(%rcx), %xmm4
movsd 0x78(%rcx), %xmm5
movsd 0x58(%rcx), %xmm6
movsd 0x70(%rcx), %xmm7
movapd %xmm0, %xmm8
addsd %xmm1, %xmm8
movapd %xmm2, %xmm9
subsd %xmm3, %xmm9
subsd %xmm1, %xmm0
addsd %xmm3, %xmm2
movapd %xmm4, %xmm1
addsd %xmm5, %xmm1
movapd %xmm6, %xmm3
subsd %xmm7, %xmm3
subsd %xmm5, %xmm4
addsd %xmm7, %xmm6
movsd %xmm8, 0x10(%rcx)
movsd %xmm9, 0x18(%rcx)
movsd %xmm0, 0x30(%rcx)
movsd %xmm2, 0x38(%rcx)
movsd %xmm1, 0x50(%rcx)
movsd %xmm3, 0x58(%rcx)
movsd %xmm4, 0x70(%rcx)
movsd %xmm6, 0x78(%rcx)
subq $-0x80, %rcx
decq %rax
jne 0x250f5
movl $0x4, %r12d
cmpq $0x2, 0x20(%rsp)
jne 0x2522c
movq %rbx, %rdi
movq %r15, %rsi
movq %r12, %rdx
callq 0x274fa
shll $0x2, %r12d
movl $0x1, %ecx
movq %rbx, %rdi
movq %r15, %rsi
cmpq $0xa, %r15
ja 0x25251
movq 0x18(%rsp), %rdx
movq %r12, %r8
movq 0x10(%rsp), %r9
callq 0x2622d
jmp 0x25263
movq 0x18(%rsp), %rdx
movq %r12, %r8
movq 0x10(%rsp), %r9
callq 0x24b5c
movq 0x30(%rsp), %rax
leaq (%rbx,%rax,8), %rbx
leaq -0x1(%r14), %rax
cmpq $0x1, %r14
movq %rax, %r14
movq %r13, %r8
jg 0x250ce
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq