	.file	"108029.c"
	.intel_syntax noprefix
	.text
	.globl	get_tau
	.type	get_tau, @function
get_tau:
.LFB6:
	.cfi_startproc
	endbr64
	push	rbp
	.cfi_def_cfa_offset 16
	.cfi_offset 6, -16
	mov	rbp, rsp
	.cfi_def_cfa_register 6
	sub	rsp, 208
	mov	QWORD PTR -120[rbp], rdi
	mov	QWORD PTR -128[rbp], rsi
	movsd	QWORD PTR -136[rbp], xmm0
	movsd	QWORD PTR -144[rbp], xmm1
	mov	DWORD PTR -148[rbp], edx
	movsd	QWORD PTR -160[rbp], xmm2
	movsd	QWORD PTR -168[rbp], xmm3
	mov	QWORD PTR -176[rbp], rcx
	mov	QWORD PTR -184[rbp], r8
	mov	QWORD PTR -192[rbp], r9
	mov	rax, QWORD PTR fs:40
	mov	QWORD PTR -8[rbp], rax
	xor	eax, eax
	movss	xmm0, DWORD PTR .LC0[rip]
	movss	DWORD PTR -76[rbp], xmm0
	mov	rax, QWORD PTR -168[rbp]
	movq	xmm0, rax
	call	nint@PLT
	mov	DWORD PTR -72[rbp], eax
	pxor	xmm0, xmm0
	cvtsi2sd	xmm0, DWORD PTR -72[rbp]
	subsd	xmm0, QWORD PTR -168[rbp]
	movq	xmm1, QWORD PTR .LC1[rip]
	andpd	xmm0, xmm1
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -76[rbp]
	comisd	xmm0, xmm1
	jbe	.L22
	mov	eax, 1
	jmp	.L16
.L22:
	lea	r8, -104[rbp]
	lea	rdi, -92[rbp]
	lea	rdx, -100[rbp]
	mov	esi, DWORD PTR -148[rbp]
	mov	rax, QWORD PTR -120[rbp]
	sub	rsp, 8
	lea	rcx, -104[rbp]
	push	rcx
	mov	r9, r8
	mov	r8, rdi
	mov	rcx, rdx
	mov	edx, 2
	mov	rdi, rax
	call	stats@PLT
	add	rsp, 16
	mov	DWORD PTR -80[rbp], 0
	jmp	.L5
.L6:
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -120[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	xmm1, DWORD PTR -100[rbp]
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -120[rbp]
	add	rax, rdx
	subss	xmm0, xmm1
	movss	DWORD PTR [rax], xmm0
	add	DWORD PTR -80[rbp], 1
.L5:
	mov	eax, DWORD PTR -80[rbp]
	cmp	eax, DWORD PTR -148[rbp]
	jl	.L6
	lea	r8, -104[rbp]
	lea	rdi, -88[rbp]
	lea	rdx, -96[rbp]
	mov	esi, DWORD PTR -148[rbp]
	mov	rax, QWORD PTR -128[rbp]
	sub	rsp, 8
	lea	rcx, -104[rbp]
	push	rcx
	mov	r9, r8
	mov	r8, rdi
	mov	rcx, rdx
	mov	edx, 2
	mov	rdi, rax
	call	stats@PLT
	add	rsp, 16
	mov	DWORD PTR -80[rbp], 0
	jmp	.L7
.L8:
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -128[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	xmm1, DWORD PTR -96[rbp]
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -128[rbp]
	add	rax, rdx
	subss	xmm0, xmm1
	movss	DWORD PTR [rax], xmm0
	add	DWORD PTR -80[rbp], 1
.L7:
	mov	eax, DWORD PTR -80[rbp]
	cmp	eax, DWORD PTR -148[rbp]
	jl	.L8
	movsd	xmm0, QWORD PTR -160[rbp]
	mulsd	xmm0, QWORD PTR -168[rbp]
	movq	rax, xmm0
	movq	xmm0, rax
	call	nint@PLT
	mov	DWORD PTR -68[rbp], eax
	mov	eax, DWORD PTR -68[rbp]
	cmp	eax, DWORD PTR -148[rbp]
	jle	.L9
	mov	eax, 2
	jmp	.L16
.L9:
	mov	eax, DWORD PTR -68[rbp]
	add	eax, eax
	add	eax, 1
	mov	DWORD PTR -64[rbp], eax
	mov	eax, DWORD PTR -64[rbp]
	cdqe
	sal	rax, 2
	mov	rdi, rax
	call	malloc@PLT
	mov	QWORD PTR -16[rbp], rax
	movsd	xmm0, QWORD PTR .LC2[rip]
	divsd	xmm0, QWORD PTR -168[rbp]
	cvtsd2ss	xmm0, xmm0
	movss	DWORD PTR -60[rbp], xmm0
	lea	r9, -84[rbp]
	mov	r8d, DWORD PTR -68[rbp]
	mov	ecx, DWORD PTR -148[rbp]
	mov	rdx, QWORD PTR -16[rbp]
	mov	rsi, QWORD PTR -128[rbp]
	mov	rax, QWORD PTR -120[rbp]
	sub	rsp, 8
	lea	rdi, -108[rbp]
	push	rdi
	mov	rdi, rax
	call	correl
	add	rsp, 16
	mov	DWORD PTR -56[rbp], eax
	movss	xmm0, DWORD PTR -84[rbp]
	movss	xmm2, DWORD PTR -92[rbp]
	movss	xmm1, DWORD PTR -88[rbp]
	mulss	xmm1, xmm2
	divss	xmm0, xmm1
	mov	rax, QWORD PTR -184[rbp]
	movss	DWORD PTR [rax], xmm0
	movss	xmm0, DWORD PTR -84[rbp]
	movss	xmm2, DWORD PTR -92[rbp]
	movss	xmm1, DWORD PTR -92[rbp]
	mulss	xmm1, xmm2
	divss	xmm0, xmm1
	mov	rax, QWORD PTR -192[rbp]
	movss	DWORD PTR [rax], xmm0
	mov	rax, QWORD PTR -184[rbp]
	movss	xmm0, DWORD PTR [rax]
	movss	xmm1, DWORD PTR .LC3[rip]
	comiss	xmm0, xmm1
	jbe	.L10
	mov	rax, QWORD PTR -184[rbp]
	movss	xmm0, DWORD PTR .LC3[rip]
	movss	DWORD PTR [rax], xmm0
.L10:
	mov	rax, QWORD PTR -184[rbp]
	movss	xmm1, DWORD PTR [rax]
	movss	xmm0, DWORD PTR .LC4[rip]
	comiss	xmm0, xmm1
	jbe	.L12
	mov	rax, QWORD PTR -184[rbp]
	movss	xmm0, DWORD PTR .LC4[rip]
	movss	DWORD PTR [rax], xmm0
.L12:
	mov	eax, DWORD PTR -108[rbp]
	pxor	xmm0, xmm0
	cvtsi2ss	xmm0, eax
	movss	xmm1, DWORD PTR -60[rbp]
	mulss	xmm0, xmm1
	movss	DWORD PTR -52[rbp], xmm0
	mov	edx, DWORD PTR -108[rbp]
	mov	eax, DWORD PTR -68[rbp]
	add	eax, edx
	cdqe
	sal	rax, 2
	lea	rdx, -4[rax]
	mov	rax, QWORD PTR -16[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	DWORD PTR -48[rbp], xmm0
	mov	edx, DWORD PTR -108[rbp]
	mov	eax, DWORD PTR -68[rbp]
	add	eax, edx
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -16[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	DWORD PTR -44[rbp], xmm0
	mov	edx, DWORD PTR -108[rbp]
	mov	eax, DWORD PTR -68[rbp]
	add	eax, edx
	cdqe
	add	rax, 1
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -16[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	DWORD PTR -40[rbp], xmm0
	movss	xmm0, DWORD PTR -48[rbp]
	movss	xmm2, DWORD PTR .LC5[rip]
	movaps	xmm1, xmm0
	divss	xmm1, xmm2
	movss	xmm0, DWORD PTR -40[rbp]
	movss	xmm2, DWORD PTR .LC5[rip]
	divss	xmm0, xmm2
	addss	xmm0, xmm1
	subss	xmm0, DWORD PTR -44[rbp]
	pxor	xmm4, xmm4
	cvtss2sd	xmm4, xmm0
	movsd	QWORD PTR -200[rbp], xmm4
	pxor	xmm5, xmm5
	cvtss2sd	xmm5, DWORD PTR -60[rbp]
	movq	rax, xmm5
	movsd	xmm0, QWORD PTR .LC6[rip]
	movapd	xmm1, xmm0
	movq	xmm0, rax
	call	pow@PLT
	movsd	xmm1, QWORD PTR -200[rbp]
	divsd	xmm1, xmm0
	pxor	xmm0, xmm0
	cvtsd2ss	xmm0, xmm1
	movss	DWORD PTR -36[rbp], xmm0
	movss	xmm0, DWORD PTR -48[rbp]
	movss	xmm1, DWORD PTR .LC7[rip]
	xorps	xmm0, xmm1
	movss	xmm2, DWORD PTR .LC5[rip]
	movaps	xmm1, xmm0
	divss	xmm1, xmm2
	movss	xmm0, DWORD PTR -40[rbp]
	movss	xmm2, DWORD PTR .LC5[rip]
	divss	xmm0, xmm2
	addss	xmm0, xmm1
	divss	xmm0, DWORD PTR -60[rbp]
	movss	DWORD PTR -32[rbp], xmm0
	movss	xmm0, DWORD PTR -44[rbp]
	movss	DWORD PTR -28[rbp], xmm0
	pxor	xmm0, xmm0
	ucomiss	xmm0, DWORD PTR -36[rbp]
	jp	.L14
	pxor	xmm0, xmm0
	ucomiss	xmm0, DWORD PTR -36[rbp]
	jne	.L14
	mov	rax, QWORD PTR -16[rbp]
	mov	rdi, rax
	call	free@PLT
	mov	eax, 3
	jmp	.L16
.L14:
	movss	xmm0, DWORD PTR -32[rbp]
	movss	xmm1, DWORD PTR .LC7[rip]
	xorps	xmm1, xmm0
	movss	xmm0, DWORD PTR -36[rbp]
	movaps	xmm2, xmm0
	addss	xmm2, xmm0
	divss	xmm1, xmm2
	movaps	xmm0, xmm1
	movss	DWORD PTR -24[rbp], xmm0
	movss	xmm0, DWORD PTR -52[rbp]
	addss	xmm0, DWORD PTR -24[rbp]
	movss	DWORD PTR -20[rbp], xmm0
	movsd	xmm0, QWORD PTR -144[rbp]
	movapd	xmm1, xmm0
	subsd	xmm1, QWORD PTR -136[rbp]
	pxor	xmm0, xmm0
	cvtss2sd	xmm0, DWORD PTR -20[rbp]
	addsd	xmm0, xmm1
	mov	rax, QWORD PTR -176[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR -16[rbp]
	mov	rdi, rax
	call	free@PLT
	mov	eax, 0
.L16:
	mov	rcx, QWORD PTR -8[rbp]
	sub	rcx, QWORD PTR fs:40
	je	.L17
	call	__stack_chk_fail@PLT
.L17:
	leave
	.cfi_def_cfa 7, 8
	ret
	.cfi_endproc
.LFE6:
	.size	get_tau, .-get_tau
	.globl	correl
	.type	correl, @function
correl:
.LFB7:
	.cfi_startproc
	endbr64
	push	rbp
	.cfi_def_cfa_offset 16
	.cfi_offset 6, -16
	mov	rbp, rsp
	.cfi_def_cfa_register 6
	mov	QWORD PTR -24[rbp], rdi
	mov	QWORD PTR -32[rbp], rsi
	mov	QWORD PTR -40[rbp], rdx
	mov	DWORD PTR -44[rbp], ecx
	mov	DWORD PTR -48[rbp], r8d
	mov	QWORD PTR -56[rbp], r9
	mov	rax, QWORD PTR -56[rbp]
	pxor	xmm0, xmm0
	movss	DWORD PTR [rax], xmm0
	mov	eax, DWORD PTR -48[rbp]
	add	eax, eax
	add	eax, 1
	mov	DWORD PTR -4[rbp], eax
	mov	DWORD PTR -12[rbp], 0
	jmp	.L24
.L30:
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	pxor	xmm0, xmm0
	movss	DWORD PTR [rax], xmm0
	mov	DWORD PTR -8[rbp], 0
	jmp	.L25
.L27:
	mov	edx, DWORD PTR -8[rbp]
	mov	eax, DWORD PTR -12[rbp]
	add	eax, edx
	sub	eax, DWORD PTR -48[rbp]
	test	eax, eax
	js	.L26
	mov	edx, DWORD PTR -8[rbp]
	mov	eax, DWORD PTR -12[rbp]
	add	eax, edx
	sub	eax, DWORD PTR -48[rbp]
	cmp	DWORD PTR -44[rbp], eax
	jle	.L26
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movss	xmm1, DWORD PTR [rax]
	mov	eax, DWORD PTR -8[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -24[rbp]
	add	rax, rdx
	movss	xmm2, DWORD PTR [rax]
	mov	edx, DWORD PTR -8[rbp]
	mov	eax, DWORD PTR -12[rbp]
	add	eax, edx
	sub	eax, DWORD PTR -48[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -32[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	mulss	xmm0, xmm2
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	addss	xmm0, xmm1
	movss	DWORD PTR [rax], xmm0
.L26:
	add	DWORD PTR -8[rbp], 1
.L25:
	mov	eax, DWORD PTR -8[rbp]
	cmp	eax, DWORD PTR -44[rbp]
	jl	.L27
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	pxor	xmm1, xmm1
	cvtsi2ss	xmm1, DWORD PTR -44[rbp]
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	divss	xmm0, xmm1
	movss	DWORD PTR [rax], xmm0
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	movss	xmm1, DWORD PTR .LC9[rip]
	andps	xmm0, xmm1
	mov	rax, QWORD PTR -56[rbp]
	movss	xmm1, DWORD PTR [rax]
	movss	xmm2, DWORD PTR .LC9[rip]
	andps	xmm1, xmm2
	comiss	xmm0, xmm1
	jbe	.L28
	mov	eax, DWORD PTR -12[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	mov	rax, QWORD PTR -56[rbp]
	movss	DWORD PTR [rax], xmm0
	mov	eax, DWORD PTR -12[rbp]
	sub	eax, DWORD PTR -48[rbp]
	mov	edx, eax
	mov	rax, QWORD PTR 16[rbp]
	mov	DWORD PTR [rax], edx
.L28:
	add	DWORD PTR -12[rbp], 1
.L24:
	mov	eax, DWORD PTR -12[rbp]
	cmp	eax, DWORD PTR -4[rbp]
	jl	.L30
	mov	eax, 0
	pop	rbp
	.cfi_def_cfa 7, 8
	ret
	.cfi_endproc
.LFE7:
	.size	correl, .-correl
	.section	.rodata
	.align 4
.LC0:
	.long	1008981770
	.align 16
.LC1:
	.long	-1
	.long	2147483647
	.long	0
	.long	0
	.align 8
.LC2:
	.long	0
	.long	1072693248
	.align 4
.LC3:
	.long	1065353216
	.align 4
.LC4:
	.long	-1082130432
	.align 4
.LC5:
	.long	1073741824
	.align 8
.LC6:
	.long	0
	.long	1073741824
	.align 16
.LC7:
	.long	-2147483648
	.long	0
	.long	0
	.long	0
	.align 16
.LC9:
	.long	2147483647
	.long	0
	.long	0
	.long	0
	.ident	"GCC: (Ubuntu 10.5.0-1ubuntu1~20.04) 10.5.0"
	.section	.note.GNU-stack,"",@progbits
	.section	.note.gnu.property,"a"
	.align 8
	.long	 1f - 0f
	.long	 4f - 1f
	.long	 5
0:
	.string	 "GNU"
1:
	.align 8
	.long	 0xc0000002
	.long	 3f - 2f
2:
	.long	 0x3
3:
	.align 8
4:

	.arch armv5t
	.eabi_attribute 20, 1
	.eabi_attribute 21, 1
	.eabi_attribute 23, 3
	.eabi_attribute 24, 1
	.eabi_attribute 25, 1
	.eabi_attribute 26, 2
	.eabi_attribute 30, 6
	.eabi_attribute 34, 0
	.eabi_attribute 18, 4
	.file	"108029.c"
	.text
	.global	__aeabi_i2d
	.global	__aeabi_dsub
	.global	__aeabi_f2d
	.global	__aeabi_dcmpgt
	.global	__aeabi_fsub
	.global	__aeabi_dmul
	.global	__aeabi_ddiv
	.global	__aeabi_d2f
	.global	__aeabi_fmul
	.global	__aeabi_fdiv
	.global	__aeabi_fcmpgt
	.global	__aeabi_fcmplt
	.global	__aeabi_i2f
	.global	__aeabi_fadd
	.global	__aeabi_fcmpeq
	.global	__aeabi_dadd
	.section	.rodata
	.align	2
.LC0:
	.word	__stack_chk_guard
	.text
	.align	2
	.global	get_tau
	.syntax unified
	.arm
	.fpu softvfp
	.type	get_tau, %function
get_tau:
	@ args = 44, pretend = 0, frame = 136
	@ frame_needed = 1, uses_anonymous_args = 0
	push	{r4, r5, fp, lr}
	add	fp, sp, #12
	sub	sp, sp, #152
	str	r0, [fp, #-120]
	str	r1, [fp, #-124]
	str	r2, [fp, #-132]
	str	r3, [fp, #-128]
	ldr	r3, [fp, #36]
	str	r3, [fp, #-136]
	ldr	r3, [fp, #40]
	str	r3, [fp, #-140]
	ldr	r3, [fp, #44]
	str	r3, [fp, #-144]
	ldr	r3, .L24
	ldr	r3, [r3]
	str	r3, [fp, #-16]
	mov	r3,#0
	ldr	r3, .L24+4
	str	r3, [fp, #-80]	@ float
	add	r1, fp, #28
	ldmia	r1, {r0-r1}
	bl	nint
	str	r0, [fp, #-76]
	ldr	r0, [fp, #-76]
	bl	__aeabi_i2d
	add	r3, fp, #28
	ldmia	r3, {r2-r3}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	bic	r5, r3, #-2147483648
	ldr	r0, [fp, #-80]	@ float
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dcmpgt
	mov	r3, r0
	cmp	r3, #0
	beq	.L22
	mov	r3, #1
	b	.L16
.L22:
	sub	r2, fp, #104
	sub	r3, fp, #108
	str	r3, [sp, #8]
	sub	r3, fp, #108
	str	r3, [sp, #4]
	sub	r3, fp, #96
	str	r3, [sp]
	mov	r3, r2
	mov	r2, #2
	ldr	r1, [fp, #12]
	ldr	r0, [fp, #-120]
	bl	stats
	mov	r3, #0
	str	r3, [fp, #-84]
	b	.L5
.L6:
	ldr	r3, [fp, #-84]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-120]
	add	r3, r2, r3
	ldr	r0, [r3]	@ float
	ldr	r1, [fp, #-104]	@ float
	ldr	r3, [fp, #-84]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-120]
	add	r4, r2, r3
	bl	__aeabi_fsub
	mov	r3, r0
	str	r3, [r4]	@ float
	ldr	r3, [fp, #-84]
	add	r3, r3, #1
	str	r3, [fp, #-84]
.L5:
	ldr	r2, [fp, #-84]
	ldr	r3, [fp, #12]
	cmp	r2, r3
	blt	.L6
	sub	r2, fp, #100
	sub	r3, fp, #108
	str	r3, [sp, #8]
	sub	r3, fp, #108
	str	r3, [sp, #4]
	sub	r3, fp, #92
	str	r3, [sp]
	mov	r3, r2
	mov	r2, #2
	ldr	r1, [fp, #12]
	ldr	r0, [fp, #-124]
	bl	stats
	mov	r3, #0
	str	r3, [fp, #-84]
	b	.L7
.L8:
	ldr	r3, [fp, #-84]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-124]
	add	r3, r2, r3
	ldr	r0, [r3]	@ float
	ldr	r1, [fp, #-100]	@ float
	ldr	r3, [fp, #-84]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-124]
	add	r4, r2, r3
	bl	__aeabi_fsub
	mov	r3, r0
	str	r3, [r4]	@ float
	ldr	r3, [fp, #-84]
	add	r3, r3, #1
	str	r3, [fp, #-84]
.L7:
	ldr	r2, [fp, #-84]
	ldr	r3, [fp, #12]
	cmp	r2, r3
	blt	.L8
	add	r3, fp, #28
	ldmia	r3, {r2-r3}
	add	r1, fp, #20
	ldmia	r1, {r0-r1}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	nint
	str	r0, [fp, #-72]
	ldr	r2, [fp, #-72]
	ldr	r3, [fp, #12]
	cmp	r2, r3
	ble	.L9
	mov	r3, #2
	b	.L16
.L9:
	ldr	r3, [fp, #-72]
	lsl	r3, r3, #1
	add	r3, r3, #1
	str	r3, [fp, #-68]
	ldr	r3, [fp, #-68]
	lsl	r3, r3, #2
	mov	r0, r3
	bl	malloc
	mov	r3, r0
	str	r3, [fp, #-64]
	add	r3, fp, #28
	ldmia	r3, {r2-r3}
	mov	r0, #0
	ldr	r1, .L24+8
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	__aeabi_d2f
	mov	r3, r0
	str	r3, [fp, #-60]	@ float
	sub	r3, fp, #112
	str	r3, [sp, #8]
	sub	r3, fp, #88
	str	r3, [sp, #4]
	ldr	r3, [fp, #-72]
	str	r3, [sp]
	ldr	r3, [fp, #12]
	ldr	r2, [fp, #-64]
	ldr	r1, [fp, #-124]
	ldr	r0, [fp, #-120]
	bl	correl
	str	r0, [fp, #-56]
	ldr	r4, [fp, #-88]	@ float
	ldr	r3, [fp, #-96]	@ float
	ldr	r2, [fp, #-92]	@ float
	mov	r1, r2
	mov	r0, r3
	bl	__aeabi_fmul
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r2, r3
	ldr	r3, [fp, #-140]
	str	r2, [r3]	@ float
	ldr	r4, [fp, #-88]	@ float
	ldr	r3, [fp, #-96]	@ float
	ldr	r2, [fp, #-96]	@ float
	mov	r1, r2
	mov	r0, r3
	bl	__aeabi_fmul
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r2, r3
	ldr	r3, [fp, #-144]
	str	r2, [r3]	@ float
	ldr	r3, [fp, #-140]
	ldr	r3, [r3]	@ float
	mov	r1, #1065353216
	mov	r0, r3
	bl	__aeabi_fcmpgt
	mov	r3, r0
	cmp	r3, #0
	beq	.L10
	ldr	r3, [fp, #-140]
	mov	r2, #1065353216
	str	r2, [r3]	@ float
.L10:
	ldr	r3, [fp, #-140]
	ldr	r3, [r3]	@ float
	ldr	r1, .L24+12
	mov	r0, r3
	bl	__aeabi_fcmplt
	mov	r3, r0
	cmp	r3, #0
	beq	.L12
	ldr	r3, [fp, #-140]
	ldr	r2, .L24+12
	str	r2, [r3]	@ float
.L12:
	ldr	r3, [fp, #-112]
	mov	r0, r3
	bl	__aeabi_i2f
	mov	r3, r0
	mov	r1, r3
	ldr	r0, [fp, #-60]	@ float
	bl	__aeabi_fmul
	mov	r3, r0
	str	r3, [fp, #-52]	@ float
	ldr	r2, [fp, #-112]
	ldr	r3, [fp, #-72]
	add	r3, r2, r3
	sub	r3, r3, #-1073741823
	lsl	r3, r3, #2
	ldr	r2, [fp, #-64]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	str	r3, [fp, #-48]	@ float
	ldr	r2, [fp, #-112]
	ldr	r3, [fp, #-72]
	add	r3, r2, r3
	lsl	r3, r3, #2
	ldr	r2, [fp, #-64]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	str	r3, [fp, #-44]	@ float
	ldr	r2, [fp, #-112]
	ldr	r3, [fp, #-72]
	add	r3, r2, r3
	add	r3, r3, #1
	lsl	r3, r3, #2
	ldr	r2, [fp, #-64]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	str	r3, [fp, #-40]	@ float
	mov	r1, #1073741824
	ldr	r0, [fp, #-48]	@ float
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r4, r3
	mov	r1, #1073741824
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fadd
	mov	r3, r0
	ldr	r1, [fp, #-44]	@ float
	mov	r0, r3
	bl	__aeabi_fsub
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	mov	r4, r0
	mov	r5, r1
	ldr	r0, [fp, #-60]	@ float
	bl	__aeabi_f2d
	mov	r2, #0
	mov	r3, #1073741824
	bl	pow
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	__aeabi_d2f
	mov	r3, r0
	str	r3, [fp, #-36]	@ float
	ldr	r3, [fp, #-48]
	eor	r3, r3, #-2147483648
	mov	r1, #1073741824
	mov	r0, r3
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r4, r3
	mov	r1, #1073741824
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fadd
	mov	r3, r0
	ldr	r1, [fp, #-60]	@ float
	mov	r0, r3
	bl	__aeabi_fdiv
	mov	r3, r0
	str	r3, [fp, #-32]	@ float
	ldr	r3, [fp, #-44]	@ float
	str	r3, [fp, #-28]	@ float
	mov	r1, #0
	ldr	r0, [fp, #-36]	@ float
	bl	__aeabi_fcmpeq
	mov	r3, r0
	cmp	r3, #0
	beq	.L23
	ldr	r0, [fp, #-64]
	bl	free
	mov	r3, #3
	b	.L16
.L23:
	ldr	r3, [fp, #-32]
	eor	r4, r3, #-2147483648
	ldr	r3, [fp, #-36]	@ float
	mov	r1, r3
	mov	r0, r3
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fdiv
	mov	r3, r0
	str	r3, [fp, #-24]	@ float
	ldr	r1, [fp, #-24]	@ float
	ldr	r0, [fp, #-52]	@ float
	bl	__aeabi_fadd
	mov	r3, r0
	str	r3, [fp, #-20]	@ float
	sub	r3, fp, #132
	ldmia	r3, {r2-r3}
	ldmib	fp, {r0-r1}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-20]	@ float
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #-136]
	stm	r1, {r2-r3}
	ldr	r0, [fp, #-64]
	bl	free
	mov	r3, #0
.L16:
	ldr	r2, .L24
	ldr	r1, [r2]
	ldr	r2, [fp, #-16]
	eors	r1, r2, r1
	mov	r2, #0
	beq	.L17
	bl	__stack_chk_fail
.L17:
	mov	r0, r3
	sub	sp, fp, #12
	@ sp needed
	pop	{r4, r5, fp, pc}
.L25:
	.align	2
.L24:
	.word	.LC0
	.word	1008981770
	.word	1072693248
	.word	-1082130432
	.size	get_tau, .-get_tau
	.align	2
	.global	correl
	.syntax unified
	.arm
	.fpu softvfp
	.type	correl, %function
correl:
	@ args = 12, pretend = 0, frame = 32
	@ frame_needed = 1, uses_anonymous_args = 0
	push	{r4, r5, fp, lr}
	add	fp, sp, #12
	sub	sp, sp, #32
	str	r0, [fp, #-32]
	str	r1, [fp, #-36]
	str	r2, [fp, #-40]
	str	r3, [fp, #-44]
	ldr	r3, [fp, #8]
	mov	r2, #0
	str	r2, [r3]	@ float
	ldr	r3, [fp, #4]
	lsl	r3, r3, #1
	add	r3, r3, #1
	str	r3, [fp, #-16]
	mov	r3, #0
	str	r3, [fp, #-24]
	b	.L27
.L33:
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r3, r2, r3
	mov	r2, #0
	str	r2, [r3]	@ float
	mov	r3, #0
	str	r3, [fp, #-20]
	b	.L28
.L30:
	ldr	r2, [fp, #-20]
	ldr	r3, [fp, #-24]
	add	r2, r2, r3
	ldr	r3, [fp, #4]
	sub	r3, r2, r3
	cmp	r3, #0
	blt	.L29
	ldr	r2, [fp, #-20]
	ldr	r3, [fp, #-24]
	add	r2, r2, r3
	ldr	r3, [fp, #4]
	sub	r3, r2, r3
	ldr	r2, [fp, #-44]
	cmp	r2, r3
	ble	.L29
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r3, r2, r3
	ldr	r5, [r3]	@ float
	ldr	r3, [fp, #-20]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-32]
	add	r3, r2, r3
	ldr	r0, [r3]	@ float
	ldr	r2, [fp, #-20]
	ldr	r3, [fp, #-24]
	add	r2, r2, r3
	ldr	r3, [fp, #4]
	sub	r3, r2, r3
	lsl	r3, r3, #2
	ldr	r2, [fp, #-36]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	mov	r1, r3
	bl	__aeabi_fmul
	mov	r3, r0
	mov	r1, r3
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r4, r2, r3
	mov	r0, r5
	bl	__aeabi_fadd
	mov	r3, r0
	str	r3, [r4]	@ float
.L29:
	ldr	r3, [fp, #-20]
	add	r3, r3, #1
	str	r3, [fp, #-20]
.L28:
	ldr	r2, [fp, #-20]
	ldr	r3, [fp, #-44]
	cmp	r2, r3
	blt	.L30
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r3, r2, r3
	ldr	r5, [r3]	@ float
	ldr	r0, [fp, #-44]
	bl	__aeabi_i2f
	mov	r1, r0
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r4, r2, r3
	mov	r0, r5
	bl	__aeabi_fdiv
	mov	r3, r0
	str	r3, [r4]	@ float
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	bic	r2, r3, #-2147483648
	ldr	r3, [fp, #8]
	ldr	r3, [r3]	@ float
	bic	r3, r3, #-2147483648
	mov	r1, r3
	mov	r0, r2
	bl	__aeabi_fcmpgt
	mov	r3, r0
	cmp	r3, #0
	beq	.L31
	ldr	r3, [fp, #-24]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-40]
	add	r3, r2, r3
	ldr	r2, [r3]	@ float
	ldr	r3, [fp, #8]
	str	r2, [r3]	@ float
	ldr	r2, [fp, #-24]
	ldr	r3, [fp, #4]
	sub	r2, r2, r3
	ldr	r3, [fp, #12]
	str	r2, [r3]
.L31:
	ldr	r3, [fp, #-24]
	add	r3, r3, #1
	str	r3, [fp, #-24]
.L27:
	ldr	r2, [fp, #-24]
	ldr	r3, [fp, #-16]
	cmp	r2, r3
	blt	.L33
	mov	r3, #0
	mov	r0, r3
	sub	sp, fp, #12
	@ sp needed
	pop	{r4, r5, fp, pc}
	.size	correl, .-correl
	.ident	"GCC: (Ubuntu 9.4.0-1ubuntu1~20.04.2) 9.4.0"
	.section	.note.GNU-stack,"",%progbits
