	.file	"125086.c"
	.intel_syntax noprefix
	.text
	.section	.rodata
.LC2:
	.string	"Bad x in routine betai"
	.text
	.globl	betai
	.type	betai, @function
betai:
.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, 56
	.cfi_offset 3, -24
	movss	DWORD PTR -36[rbp], xmm0
	movss	DWORD PTR -40[rbp], xmm1
	movss	DWORD PTR -44[rbp], xmm2
	pxor	xmm0, xmm0
	comiss	xmm0, DWORD PTR -44[rbp]
	ja	.L2
	movss	xmm0, DWORD PTR -44[rbp]
	movss	xmm1, DWORD PTR .LC1[rip]
	comiss	xmm0, xmm1
	jbe	.L3
.L2:
	lea	rdi, .LC2[rip]
	call	nrerror@PLT
.L3:
	pxor	xmm0, xmm0
	ucomiss	xmm0, DWORD PTR -44[rbp]
	jp	.L13
	pxor	xmm0, xmm0
	ucomiss	xmm0, DWORD PTR -44[rbp]
	je	.L5
.L13:
	movss	xmm0, DWORD PTR .LC1[rip]
	ucomiss	xmm0, DWORD PTR -44[rbp]
	jp	.L7
	movss	xmm0, DWORD PTR .LC1[rip]
	ucomiss	xmm0, DWORD PTR -44[rbp]
	jne	.L7
.L5:
	pxor	xmm0, xmm0
	movss	DWORD PTR -20[rbp], xmm0
	jmp	.L9
.L7:
	movss	xmm0, DWORD PTR -36[rbp]
	addss	xmm0, DWORD PTR -40[rbp]
	movd	eax, xmm0
	movd	xmm0, eax
	call	gammln@PLT
	movd	ebx, xmm0
	mov	eax, DWORD PTR -36[rbp]
	movd	xmm0, eax
	call	gammln@PLT
	movd	xmm4, ebx
	subss	xmm4, xmm0
	movss	DWORD PTR -56[rbp], xmm4
	mov	eax, DWORD PTR -40[rbp]
	movd	xmm0, eax
	call	gammln@PLT
	movss	xmm1, DWORD PTR -56[rbp]
	subss	xmm1, xmm0
	pxor	xmm3, xmm3
	cvtss2sd	xmm3, xmm1
	movsd	QWORD PTR -56[rbp], xmm3
	pxor	xmm5, xmm5
	cvtss2sd	xmm5, DWORD PTR -36[rbp]
	movsd	QWORD PTR -64[rbp], xmm5
	pxor	xmm4, xmm4
	cvtss2sd	xmm4, DWORD PTR -44[rbp]
	movq	rax, xmm4
	movq	xmm0, rax
	call	log@PLT
	mulsd	xmm0, QWORD PTR -64[rbp]
	addsd	xmm0, QWORD PTR -56[rbp]
	movsd	QWORD PTR -56[rbp], xmm0
	pxor	xmm6, xmm6
	cvtss2sd	xmm6, DWORD PTR -40[rbp]
	movsd	QWORD PTR -64[rbp], xmm6
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -44[rbp]
	movsd	xmm0, QWORD PTR .LC3[rip]
	movapd	xmm6, xmm0
	subsd	xmm6, xmm1
	movq	rax, xmm6
	movq	xmm0, rax
	call	log@PLT
	mulsd	xmm0, QWORD PTR -64[rbp]
	movapd	xmm7, xmm0
	addsd	xmm7, QWORD PTR -56[rbp]
	movq	rax, xmm7
	movq	xmm0, rax
	call	exp@PLT
	cvtsd2ss	xmm0, xmm0
	movss	DWORD PTR -20[rbp], xmm0
.L9:
	pxor	xmm2, xmm2
	cvtss2sd	xmm2, DWORD PTR -44[rbp]
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, DWORD PTR -36[rbp]
	movsd	xmm0, QWORD PTR .LC3[rip]
	addsd	xmm0, xmm1
	movss	xmm1, DWORD PTR -36[rbp]
	addss	xmm1, DWORD PTR -40[rbp]
	pxor	xmm3, xmm3
	cvtss2sd	xmm3, xmm1
	movsd	xmm1, QWORD PTR .LC4[rip]
	addsd	xmm1, xmm3
	divsd	xmm0, xmm1
	comisd	xmm0, xmm2
	jbe	.L15
	movss	xmm1, DWORD PTR -44[rbp]
	movss	xmm0, DWORD PTR -40[rbp]
	mov	eax, DWORD PTR -36[rbp]
	movaps	xmm2, xmm1
	movaps	xmm1, xmm0
	movd	xmm0, eax
	call	betacf@PLT
	movd	eax, xmm0
	movd	xmm0, eax
	mulss	xmm0, DWORD PTR -20[rbp]
	divss	xmm0, DWORD PTR -36[rbp]
	jmp	.L12
.L15:
	movss	xmm0, DWORD PTR .LC1[rip]
	movaps	xmm1, xmm0
	subss	xmm1, DWORD PTR -44[rbp]
	movss	xmm0, DWORD PTR -36[rbp]
	mov	eax, DWORD PTR -40[rbp]
	movaps	xmm2, xmm1
	movaps	xmm1, xmm0
	movd	xmm0, eax
	call	betacf@PLT
	movd	eax, xmm0
	movd	xmm0, eax
	mulss	xmm0, DWORD PTR -20[rbp]
	divss	xmm0, DWORD PTR -40[rbp]
	pxor	xmm1, xmm1
	cvtss2sd	xmm1, xmm0
	movsd	xmm0, QWORD PTR .LC3[rip]
	subsd	xmm0, xmm1
	cvtsd2ss	xmm0, xmm0
.L12:
	mov	rbx, QWORD PTR -8[rbp]
	leave
	.cfi_def_cfa 7, 8
	ret
	.cfi_endproc
.LFE0:
	.size	betai, .-betai
	.section	.rodata
	.align 4
.LC1:
	.long	1065353216
	.align 8
.LC3:
	.long	0
	.long	1072693248
	.align 8
.LC4:
	.long	0
	.long	1073741824
	.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	"125086.c"
	.text
	.global	__aeabi_fcmplt
	.global	__aeabi_fcmpgt
	.section	.rodata
	.align	2
.LC0:
	.ascii	"Bad x in routine betai\000"
	.global	__aeabi_fcmpeq
	.global	__aeabi_fadd
	.global	__aeabi_fsub
	.global	__aeabi_f2d
	.global	__aeabi_dmul
	.global	__aeabi_dadd
	.global	__aeabi_dsub
	.global	__aeabi_d2f
	.global	__aeabi_ddiv
	.global	__aeabi_dcmplt
	.global	__aeabi_fmul
	.global	__aeabi_fdiv
	.text
	.align	2
	.global	betai
	.syntax unified
	.arm
	.fpu softvfp
	.type	betai, %function
betai:
	@ args = 0, pretend = 0, frame = 24
	@ frame_needed = 1, uses_anonymous_args = 0
	push	{r4, r5, r6, r7, fp, lr}
	add	fp, sp, #20
	sub	sp, sp, #24
	str	r0, [fp, #-32]	@ float
	str	r1, [fp, #-36]	@ float
	str	r2, [fp, #-40]	@ float
	mov	r1, #0
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fcmplt
	mov	r3, r0
	cmp	r3, #0
	bne	.L2
	mov	r1, #1065353216
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fcmpgt
	mov	r3, r0
	cmp	r3, #0
	beq	.L3
.L2:
	ldr	r0, .L15
	bl	nrerror
.L3:
	mov	r1, #0
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fcmpeq
	mov	r3, r0
	cmp	r3, #0
	bne	.L5
	mov	r1, #1065353216
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_fcmpeq
	mov	r3, r0
	cmp	r3, #0
	beq	.L13
.L5:
	mov	r3, #0
	str	r3, [fp, #-24]	@ float
	b	.L8
.L13:
	ldr	r1, [fp, #-36]	@ float
	ldr	r0, [fp, #-32]	@ float
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	gammln
	mov	r4, r0
	ldr	r0, [fp, #-32]	@ float
	bl	gammln
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fsub
	mov	r3, r0
	mov	r4, r3
	ldr	r0, [fp, #-36]	@ float
	bl	gammln
	mov	r3, r0
	mov	r1, r3
	mov	r0, r4
	bl	__aeabi_fsub
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	mov	r4, r0
	mov	r5, r1
	ldr	r0, [fp, #-32]	@ float
	bl	__aeabi_f2d
	mov	r6, r0
	mov	r7, r1
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	log
	mov	r2, r0
	mov	r3, r1
	mov	r0, r6
	mov	r1, r7
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r0, [fp, #-36]	@ float
	bl	__aeabi_f2d
	mov	r4, r0
	mov	r5, r1
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_f2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, #0
	ldr	r1, .L15+4
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	log
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	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, r2
	mov	r1, r3
	bl	exp
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	bl	__aeabi_d2f
	mov	r3, r0
	str	r3, [fp, #-24]	@ float
.L8:
	ldr	r0, [fp, #-40]	@ float
	bl	__aeabi_f2d
	mov	r4, r0
	mov	r5, r1
	ldr	r0, [fp, #-32]	@ float
	bl	__aeabi_f2d
	mov	r2, #0
	ldr	r3, .L15+4
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r1, [fp, #-36]	@ float
	ldr	r0, [fp, #-32]	@ float
	bl	__aeabi_fadd
	mov	r3, r0
	mov	r0, r3
	bl	__aeabi_f2d
	mov	r2, #0
	mov	r3, #1073741824
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	mov	r0, r6
	mov	r1, r7
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dcmplt
	mov	r3, r0
	cmp	r3, #0
	beq	.L14
	ldr	r2, [fp, #-40]	@ float
	ldr	r1, [fp, #-36]	@ float
	ldr	r0, [fp, #-32]	@ float
	bl	betacf
	mov	r3, r0
	ldr	r1, [fp, #-24]	@ float
	mov	r0, r3
	bl	__aeabi_fmul
	mov	r3, r0
	ldr	r1, [fp, #-32]	@ float
	mov	r0, r3
	bl	__aeabi_fdiv
	mov	r3, r0
	b	.L11
.L14:
	ldr	r1, [fp, #-40]	@ float
	mov	r0, #1065353216
	bl	__aeabi_fsub
	mov	r3, r0
	mov	r2, r3
	ldr	r1, [fp, #-32]	@ float
	ldr	r0, [fp, #-36]	@ float
	bl	betacf
	mov	r3, r0
	ldr	r1, [fp, #-24]	@ float
	mov	r0, r3
	bl	__aeabi_fmul
	mov	r3, r0
	ldr	r1, [fp, #-36]	@ float
	mov	r0, r3
	bl	__aeabi_fdiv
	mov	r3, r0
	mov	r1, r3
	mov	r0, #1065353216
	bl	__aeabi_fsub
	mov	r3, r0
.L11:
	mov	r0, r3
	sub	sp, fp, #20
	@ sp needed
	pop	{r4, r5, r6, r7, fp, pc}
.L16:
	.align	2
.L15:
	.word	.LC0
	.word	1072693248
	.size	betai, .-betai
	.ident	"GCC: (Ubuntu 9.4.0-1ubuntu1~20.04.2) 9.4.0"
	.section	.note.GNU-stack,"",%progbits
