	.file	"570890.c"
	.intel_syntax noprefix
	.text
	.section	.rodata
.LC17:
	.string	"too many iterations in gaulag"
	.text
	.globl	gaulag
	.type	gaulag, @function
gaulag:
.LFB0:
	.cfi_startproc
	endbr64
	push	rbp
	.cfi_def_cfa_offset 16
	.cfi_offset 6, -16
	mov	rbp, rsp
	.cfi_def_cfa_register 6
	push	rbx
	sub	rsp, 104
	.cfi_offset 3, -24
	mov	QWORD PTR -88[rbp], rdi
	mov	QWORD PTR -96[rbp], rsi
	mov	DWORD PTR -100[rbp], edx
	movss	DWORD PTR -104[rbp], xmm0
	mov	DWORD PTR -80[rbp], 1
	jmp	.L2
.L14:
	cmp	DWORD PTR -80[rbp], 1
	jne	.L3
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -104[rbp]
	movsd	xmm0, QWORD PTR .LC0[rip]
	addsd	xmm1, xmm0
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -104[rbp]
	movsd	xmm0, QWORD PTR .LC1[rip]
	mulsd	xmm2, xmm0
	movsd	xmm0, QWORD PTR .LC2[rip]
	addsd	xmm0, xmm2
	mulsd	xmm0, xmm1
	pxor	xmm2, xmm2
	cvtsi2sd	xmm2, DWORD PTR -100[rbp]
	movsd	xmm1, QWORD PTR .LC3[rip]
	mulsd	xmm2, xmm1
	movsd	xmm1, QWORD PTR .LC0[rip]
	addsd	xmm2, xmm1
	pxor	xmm3, xmm3
	cvtss2sd	xmm3, DWORD PTR -104[rbp]
	movsd	xmm1, QWORD PTR .LC4[rip]
	mulsd	xmm1, xmm3
	addsd	xmm1, xmm2
	divsd	xmm0, xmm1
	movsd	QWORD PTR -40[rbp], xmm0
	jmp	.L4
.L3:
	cmp	DWORD PTR -80[rbp], 2
	jne	.L5
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -104[rbp]
	movsd	xmm0, QWORD PTR .LC5[rip]
	mulsd	xmm1, xmm0
	movsd	xmm0, QWORD PTR .LC6[rip]
	addsd	xmm0, xmm1
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -104[rbp]
	movsd	xmm1, QWORD PTR .LC7[rip]
	mulsd	xmm2, xmm1
	movsd	xmm1, QWORD PTR .LC0[rip]
	addsd	xmm2, xmm1
	pxor	xmm3, xmm3
	cvtsi2sd	xmm3, DWORD PTR -100[rbp]
	movsd	xmm1, QWORD PTR .LC8[rip]
	mulsd	xmm1, xmm3
	addsd	xmm1, xmm2
	divsd	xmm0, xmm1
	movsd	xmm1, QWORD PTR -40[rbp]
	addsd	xmm0, xmm1
	movsd	QWORD PTR -40[rbp], xmm0
	jmp	.L4
.L5:
	mov	eax, DWORD PTR -80[rbp]
	sub	eax, 2
	pxor	xmm0, xmm0
	cvtsi2ss	xmm0, eax
	movss	DWORD PTR -68[rbp], xmm0
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -68[rbp]
	movsd	xmm0, QWORD PTR .LC9[rip]
	mulsd	xmm1, xmm0
	movsd	xmm0, QWORD PTR .LC0[rip]
	addsd	xmm0, xmm1
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -68[rbp]
	movsd	xmm1, QWORD PTR .LC10[rip]
	mulsd	xmm2, xmm1
	movapd	xmm1, xmm0
	divsd	xmm1, xmm2
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -68[rbp]
	movsd	xmm0, QWORD PTR .LC11[rip]
	mulsd	xmm2, xmm0
	pxor	xmm0, xmm0
	cvtss2sd	xmm0, DWORD PTR -104[rbp]
	mulsd	xmm0, xmm2
	pxor	xmm3, xmm3
	cvtss2sd	xmm3, DWORD PTR -68[rbp]
	movsd	xmm2, QWORD PTR .LC12[rip]
	mulsd	xmm3, xmm2
	movsd	xmm2, QWORD PTR .LC0[rip]
	addsd	xmm2, xmm3
	divsd	xmm0, xmm2
	addsd	xmm1, xmm0
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	sal	rax, 2
	lea	rdx, -8[rax]
	mov	rax, QWORD PTR -88[rbp]
	add	rax, rdx
	movss	xmm0, DWORD PTR [rax]
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, xmm0
	movsd	xmm0, QWORD PTR -40[rbp]
	subsd	xmm0, xmm2
	mulsd	xmm0, xmm1
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -104[rbp]
	movsd	xmm1, QWORD PTR .LC13[rip]
	mulsd	xmm2, xmm1
	movsd	xmm1, QWORD PTR .LC0[rip]
	addsd	xmm1, xmm2
	divsd	xmm0, xmm1
	movsd	xmm1, QWORD PTR -40[rbp]
	addsd	xmm0, xmm1
	movsd	QWORD PTR -40[rbp], xmm0
.L4:
	mov	DWORD PTR -76[rbp], 1
	jmp	.L6
.L12:
	movsd	xmm0, QWORD PTR .LC0[rip]
	movsd	QWORD PTR -64[rbp], xmm0
	pxor	xmm0, xmm0
	movsd	QWORD PTR -56[rbp], xmm0
	mov	DWORD PTR -72[rbp], 1
	jmp	.L7
.L8:
	movsd	xmm0, QWORD PTR -56[rbp]
	movsd	QWORD PTR -24[rbp], xmm0
	movsd	xmm0, QWORD PTR -64[rbp]
	movsd	QWORD PTR -56[rbp], xmm0
	mov	eax, DWORD PTR -72[rbp]
	add	eax, eax
	sub	eax, 1
	pxor	xmm0, xmm0
	cvtsi2ss	xmm0, eax
	addss	xmm0, DWORD PTR -104[rbp]
	cvtss2sd	xmm0, xmm0
	subsd	xmm0, QWORD PTR -40[rbp]
	mulsd	xmm0, QWORD PTR -56[rbp]
	mov	eax, DWORD PTR -72[rbp]
	sub	eax, 1
	pxor	xmm1, xmm1
	cvtsi2ss	xmm1, eax
	addss	xmm1, DWORD PTR -104[rbp]
	cvtss2sd	xmm1, xmm1
	mulsd	xmm1, QWORD PTR -24[rbp]
	subsd	xmm0, xmm1
	pxor	xmm1, xmm1
	cvtsi2sd	xmm1, DWORD PTR -72[rbp]
	divsd	xmm0, xmm1
	movsd	QWORD PTR -64[rbp], xmm0
	add	DWORD PTR -72[rbp], 1
.L7:
	mov	eax, DWORD PTR -72[rbp]
	cmp	eax, DWORD PTR -100[rbp]
	jle	.L8
	pxor	xmm0, xmm0
	cvtsi2sd	xmm0, DWORD PTR -100[rbp]
	mulsd	xmm0, QWORD PTR -64[rbp]
	pxor	xmm1, xmm1
	cvtsi2ss	xmm1, DWORD PTR -100[rbp]
	addss	xmm1, DWORD PTR -104[rbp]
	cvtss2sd	xmm1, xmm1
	mulsd	xmm1, QWORD PTR -56[rbp]
	subsd	xmm0, xmm1
	divsd	xmm0, QWORD PTR -40[rbp]
	movsd	QWORD PTR -48[rbp], xmm0
	movsd	xmm0, QWORD PTR -40[rbp]
	movsd	QWORD PTR -32[rbp], xmm0
	movsd	xmm0, QWORD PTR -64[rbp]
	movapd	xmm1, xmm0
	divsd	xmm1, QWORD PTR -48[rbp]
	movsd	xmm0, QWORD PTR -32[rbp]
	subsd	xmm0, xmm1
	movsd	QWORD PTR -40[rbp], xmm0
	movsd	xmm0, QWORD PTR -40[rbp]
	subsd	xmm0, QWORD PTR -32[rbp]
	movq	xmm1, QWORD PTR .LC15[rip]
	andpd	xmm1, xmm0
	movsd	xmm0, QWORD PTR .LC16[rip]
	comisd	xmm0, xmm1
	jnb	.L15
	add	DWORD PTR -76[rbp], 1
.L6:
	cmp	DWORD PTR -76[rbp], 10
	jle	.L12
	jmp	.L11
.L15:
	nop
.L11:
	cmp	DWORD PTR -76[rbp], 10
	jle	.L13
	lea	rdi, .LC17[rip]
	call	nrerror@PLT
.L13:
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -88[rbp]
	add	rax, rdx
	pxor	xmm0, xmm0
	cvtsd2ss	xmm0, QWORD PTR -40[rbp]
	movss	DWORD PTR [rax], xmm0
	pxor	xmm0, xmm0
	cvtsi2ss	xmm0, DWORD PTR -100[rbp]
	addss	xmm0, DWORD PTR -104[rbp]
	movd	eax, xmm0
	movd	xmm0, eax
	call	gammln@PLT
	movd	ebx, xmm0
	pxor	xmm5, xmm5
	cvtsi2ss	xmm5, DWORD PTR -100[rbp]
	movd	eax, xmm5
	movd	xmm0, eax
	call	gammln@PLT
	movd	xmm1, ebx
	subss	xmm1, xmm0
	pxor	xmm6, xmm6
	cvtss2sd	xmm6, xmm1
	movq	rax, xmm6
	movq	xmm0, rax
	call	exp@PLT
	movq	rax, xmm0
	movq	xmm0, QWORD PTR .LC18[rip]
	movq	xmm4, rax
	xorpd	xmm4, xmm0
	movapd	xmm0, xmm4
	pxor	xmm1, xmm1
	cvtsi2sd	xmm1, DWORD PTR -100[rbp]
	mulsd	xmm1, QWORD PTR -48[rbp]
	mulsd	xmm1, QWORD PTR -56[rbp]
	divsd	xmm0, xmm1
	mov	eax, DWORD PTR -80[rbp]
	cdqe
	lea	rdx, 0[0+rax*4]
	mov	rax, QWORD PTR -96[rbp]
	add	rax, rdx
	cvtsd2ss	xmm0, xmm0
	movss	DWORD PTR [rax], xmm0
	add	DWORD PTR -80[rbp], 1
.L2:
	mov	eax, DWORD PTR -80[rbp]
	cmp	eax, DWORD PTR -100[rbp]
	jle	.L14
	nop
	nop
	mov	rbx, QWORD PTR -8[rbp]
	leave
	.cfi_def_cfa 7, 8
	ret
	.cfi_endproc
.LFE0:
	.size	gaulag, .-gaulag
	.section	.rodata
	.align 8
.LC0:
	.long	0
	.long	1072693248
	.align 8
.LC1:
	.long	-687194767
	.long	1072525475
	.align 8
.LC2:
	.long	0
	.long	1074266112
	.align 8
.LC3:
	.long	858993459
	.long	1073951539
	.align 8
.LC4:
	.long	-858993459
	.long	1073532108
	.align 8
.LC5:
	.long	0
	.long	1075380224
	.align 8
.LC6:
	.long	0
	.long	1076756480
	.align 8
.LC7:
	.long	-858993459
	.long	1072483532
	.align 8
.LC8:
	.long	0
	.long	1074003968
	.align 8
.LC9:
	.long	1717986918
	.long	1074030182
	.align 8
.LC10:
	.long	1717986918
	.long	1073636966
	.align 8
.LC11:
	.long	-1030792151
	.long	1072965877
	.align 8
.LC12:
	.long	0
	.long	1074528256
	.align 8
.LC13:
	.long	858993459
	.long	1070805811
	.align 16
.LC15:
	.long	-1
	.long	2147483647
	.long	0
	.long	0
	.align 8
.LC16:
	.long	-1527193420
	.long	1025565556
	.align 16
.LC18:
	.long	0
	.long	-2147483648
	.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	"570890.c"
	.text
	.global	__aeabi_f2d
	.global	__aeabi_dadd
	.global	__aeabi_dmul
	.global	__aeabi_i2d
	.global	__aeabi_ddiv
	.global	__aeabi_i2f
	.global	__aeabi_dsub
	.global	__aeabi_fadd
	.global	__aeabi_dcmple
	.section	.rodata
	.align	2
.LC0:
	.ascii	"too many iterations in gaulag\000"
	.global	__aeabi_d2f
	.global	__aeabi_fsub
	.text
	.align	2
	.global	gaulag
	.syntax unified
	.arm
	.fpu softvfp
	.type	gaulag, %function
gaulag:
	@ args = 0, pretend = 0, frame = 88
	@ frame_needed = 1, uses_anonymous_args = 0
	push	{r4, r5, r6, r7, r8, r9, fp, lr}
	add	fp, sp, #28
	sub	sp, sp, #88
	str	r0, [fp, #-96]
	str	r1, [fp, #-100]
	str	r2, [fp, #-104]
	str	r3, [fp, #-108]	@ float
	mov	r3, #1
	str	r3, [fp, #-92]
	b	.L2
.L14:
	ldr	r3, [fp, #-92]
	cmp	r3, #1
	bne	.L3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+4
	ldr	r3, .L16+8
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16+12
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2d
	ldr	r2, .L16+16
	ldr	r3, .L16+20
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+24
	ldr	r3, .L16+28
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r6
	mov	r1, r7
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-52]
	str	r3, [fp, #-48]
	b	.L4
.L3:
	ldr	r3, [fp, #-92]
	cmp	r3, #2
	bne	.L5
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	mov	r2, #0
	ldr	r3, .L16+32
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16+36
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+24
	ldr	r3, .L16+40
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2d
	mov	r2, #0
	ldr	r3, .L16+44
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r6
	mov	r1, r7
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #52
	ldmia	r1, {r0-r1}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-52]
	str	r3, [fp, #-48]
	b	.L4
.L5:
	ldr	r3, [fp, #-92]
	sub	r3, r3, #2
	mov	r0, r3
	bl	__aeabi_i2f
	mov	r3, r0
	str	r3, [fp, #-80]	@ float
	ldr	r0, [fp, #-80]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+48
	ldr	r3, .L16+52
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-80]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+48
	ldr	r3, .L16+56
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r0, [fp, #-80]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+60
	ldr	r3, .L16+64
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-80]	@ float
	bl	__aeabi_f2d
	mov	r2, #0
	ldr	r3, .L16+68
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r0, r6
	mov	r1, r7
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r3, [fp, #-92]
	sub	r3, r3, #-1073741822
	lsl	r3, r3, #2
	ldr	r2, [fp, #-96]
	add	r3, r2, r3
	ldr	r3, [r3]	@ float
	mov	r0, r3
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #52
	ldmia	r1, {r0-r1}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-108]	@ float
	bl	__aeabi_f2d
	ldr	r2, .L16+16
	ldr	r3, .L16+72
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L16
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #52
	ldmia	r1, {r0-r1}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-52]
	str	r3, [fp, #-48]
.L4:
	mov	r3, #1
	str	r3, [fp, #-88]
	b	.L6
.L12:
	mov	r2, #0
	ldr	r3, .L16
	str	r2, [fp, #-76]
	str	r3, [fp, #-72]
	mov	r2, #0
	mov	r3, #0
	str	r2, [fp, #-68]
	str	r3, [fp, #-64]
	mov	r3, #1
	str	r3, [fp, #-84]
	b	.L7
.L8:
	sub	r3, fp, #68
	ldmia	r3, {r2-r3}
	str	r2, [fp, #-36]
	str	r3, [fp, #-32]
	sub	r3, fp, #76
	ldmia	r3, {r2-r3}
	str	r2, [fp, #-68]
	str	r3, [fp, #-64]
	ldr	r3, [fp, #-84]
	lsl	r3, r3, #1
	sub	r3, r3, #1
	mov	r0, r3
	bl	__aeabi_i2f
	mov	r3, r0
	ldr	r1, [fp, #-108]	@ float
	mov	r0, r3
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	sub	r3, fp, #52
	ldmia	r3, {r2-r3}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	sub	r3, fp, #68
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r3, [fp, #-84]
	sub	r3, r3, #1
	mov	r0, r3
	bl	__aeabi_i2f
	mov	r3, r0
	ldr	r1, [fp, #-108]	@ float
	mov	r0, r3
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	sub	r3, fp, #36
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-84]
	bl	__aeabi_i2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-76]
	str	r3, [fp, #-72]
	ldr	r3, [fp, #-84]
	add	r3, r3, #1
	str	r3, [fp, #-84]
.L7:
	ldr	r2, [fp, #-84]
	ldr	r3, [fp, #-104]
	cmp	r2, r3
	ble	.L8
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2d
	sub	r3, fp, #76
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2f
	mov	r3, r0
	ldr	r1, [fp, #-108]	@ float
	mov	r0, r3
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	sub	r3, fp, #68
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	sub	r3, fp, #52
	ldmia	r3, {r2-r3}
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-60]
	str	r3, [fp, #-56]
	sub	r3, fp, #52
	ldmia	r3, {r2-r3}
	str	r2, [fp, #-44]
	str	r3, [fp, #-40]
	sub	r3, fp, #60
	ldmia	r3, {r2-r3}
	sub	r1, fp, #76
	ldmia	r1, {r0-r1}
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #44
	ldmia	r1, {r0-r1}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-52]
	str	r3, [fp, #-48]
	sub	r3, fp, #44
	ldmia	r3, {r2-r3}
	sub	r1, fp, #52
	ldmia	r1, {r0-r1}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r8, r2
	bic	r9, r3, #-2147483648
	ldr	r2, .L16+76
	ldr	r3, .L16+80
	mov	r0, r8
	mov	r1, r9
	bl	__aeabi_dcmple
	mov	r3, r0
	cmp	r3, #0
	bne	.L15
	ldr	r3, [fp, #-88]
	add	r3, r3, #1
	str	r3, [fp, #-88]
.L6:
	ldr	r3, [fp, #-88]
	cmp	r3, #10
	ble	.L12
	b	.L11
.L15:
	nop
.L11:
	ldr	r3, [fp, #-88]
	cmp	r3, #10
	ble	.L13
	ldr	r0, .L16+84
	bl	nrerror
.L13:
	ldr	r3, [fp, #-92]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-96]
	add	r4, r2, r3
	sub	r1, fp, #52
	ldmia	r1, {r0-r1}
	bl	__aeabi_d2f
	mov	r3, r0
	str	r3, [r4]	@ float
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2f
	mov	r3, r0
	ldr	r1, [fp, #-108]	@ float
	mov	r0, r3
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	gammln
	mov	r4, r0
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2f
	mov	r3, r0
	mov	r0, r3
	bl	gammln
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fsub
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	exp
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-116]
	eor	r3, r3, #-2147483648
	str	r3, [fp, #-112]
	ldr	r0, [fp, #-104]
	bl	__aeabi_i2d
	sub	r3, fp, #60
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	sub	r3, fp, #68
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #116
	ldmia	r1, {r0-r1}
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	ldr	r3, [fp, #-92]
	lsl	r3, r3, #2
	ldr	r2, [fp, #-100]
	add	r4, r2, r3
	bl	__aeabi_d2f
	mov	r3, r0
	str	r3, [r4]	@ float
	ldr	r3, [fp, #-92]
	add	r3, r3, #1
	str	r3, [fp, #-92]
.L2:
	ldr	r2, [fp, #-92]
	ldr	r3, [fp, #-104]
	cmp	r2, r3
	ble	.L14
	nop
	nop
	sub	sp, fp, #28
	@ sp needed
	pop	{r4, r5, r6, r7, r8, r9, fp, pc}
.L17:
	.align	2
.L16:
	.word	1072693248
	.word	-687194767
	.word	1072525475
	.word	1074266112
	.word	858993459
	.word	1073951539
	.word	-858993459
	.word	1073532108
	.word	1075380224
	.word	1076756480
	.word	1072483532
	.word	1074003968
	.word	1717986918
	.word	1074030182
	.word	1073636966
	.word	-1030792151
	.word	1072965877
	.word	1074528256
	.word	1070805811
	.word	-1527193420
	.word	1025565556
	.word	.LC0
	.size	gaulag, .-gaulag
	.ident	"GCC: (Ubuntu 9.4.0-1ubuntu1~20.04.2) 9.4.0"
	.section	.note.GNU-stack,"",%progbits
