	.file	"407128.c"
	.intel_syntax noprefix
	.text
	.section	.rodata
	.align 8
.LC1:
	.string	"n must be at least 2 in moment"
	.align 8
.LC4:
	.string	"No skew/kurtosis when variance = 0 (in moment)"
	.text
	.globl	moment
	.type	moment, @function
moment:
.LFB0:
	.cfi_startproc
	endbr64
	push	rbp
	.cfi_def_cfa_offset 16
	.cfi_offset 6, -16
	mov	rbp, rsp
	.cfi_def_cfa_register 6
	sub	rsp, 80
	mov	QWORD PTR -40[rbp], rdi
	mov	DWORD PTR -44[rbp], esi
	mov	QWORD PTR -56[rbp], rdx
	mov	QWORD PTR -64[rbp], rcx
	mov	QWORD PTR -72[rbp], r8
	mov	QWORD PTR -80[rbp], r9
	pxor	xmm0, xmm0
	movsd	QWORD PTR -24[rbp], xmm0
	cmp	DWORD PTR -44[rbp], 1
	jg	.L2
	lea	rdi, .LC1[rip]
	call	nrerror@PLT
.L2:
	pxor	xmm0, xmm0
	movsd	QWORD PTR -16[rbp], xmm0
	mov	DWORD PTR -28[rbp], 1
	jmp	.L3
.L4:
	mov	eax, DWORD PTR -28[rbp]
	cdqe
	lea	rdx, 0[0+rax*8]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movsd	xmm0, QWORD PTR [rax]
	movsd	xmm1, QWORD PTR -16[rbp]
	addsd	xmm0, xmm1
	movsd	QWORD PTR -16[rbp], xmm0
	add	DWORD PTR -28[rbp], 1
.L3:
	mov	eax, DWORD PTR -28[rbp]
	cmp	eax, DWORD PTR -44[rbp]
	jle	.L4
	pxor	xmm1, xmm1
	cvtsi2sd	xmm1, DWORD PTR -44[rbp]
	movsd	xmm0, QWORD PTR -16[rbp]
	divsd	xmm0, xmm1
	mov	rax, QWORD PTR -56[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR 24[rbp]
	pxor	xmm0, xmm0
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR 24[rbp]
	movsd	xmm0, QWORD PTR [rax]
	mov	rax, QWORD PTR 16[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR 16[rbp]
	movsd	xmm0, QWORD PTR [rax]
	mov	rax, QWORD PTR -80[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm0, QWORD PTR [rax]
	mov	rax, QWORD PTR -64[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	DWORD PTR -28[rbp], 1
	jmp	.L5
.L6:
	mov	eax, DWORD PTR -28[rbp]
	cdqe
	lea	rdx, 0[0+rax*8]
	mov	rax, QWORD PTR -40[rbp]
	add	rax, rdx
	movsd	xmm0, QWORD PTR [rax]
	mov	rax, QWORD PTR -56[rbp]
	movsd	xmm1, QWORD PTR [rax]
	subsd	xmm0, xmm1
	movsd	QWORD PTR -16[rbp], xmm0
	movsd	xmm0, QWORD PTR -16[rbp]
	movq	xmm1, QWORD PTR .LC2[rip]
	andpd	xmm1, xmm0
	mov	rax, QWORD PTR -64[rbp]
	movsd	xmm0, QWORD PTR [rax]
	addsd	xmm0, xmm1
	mov	rax, QWORD PTR -64[rbp]
	movsd	QWORD PTR [rax], xmm0
	movsd	xmm0, QWORD PTR -24[rbp]
	addsd	xmm0, QWORD PTR -16[rbp]
	movsd	QWORD PTR -24[rbp], xmm0
	movsd	xmm0, QWORD PTR -16[rbp]
	mulsd	xmm0, xmm0
	movsd	QWORD PTR -8[rbp], xmm0
	movsd	xmm1, QWORD PTR -8[rbp]
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm0, QWORD PTR [rax]
	addsd	xmm0, xmm1
	mov	rax, QWORD PTR -80[rbp]
	movsd	QWORD PTR [rax], xmm0
	movsd	xmm0, QWORD PTR -8[rbp]
	mulsd	xmm0, QWORD PTR -16[rbp]
	movsd	QWORD PTR -8[rbp], xmm0
	movsd	xmm1, QWORD PTR -8[rbp]
	mov	rax, QWORD PTR 16[rbp]
	movsd	xmm0, QWORD PTR [rax]
	addsd	xmm0, xmm1
	mov	rax, QWORD PTR 16[rbp]
	movsd	QWORD PTR [rax], xmm0
	movsd	xmm0, QWORD PTR -8[rbp]
	mulsd	xmm0, QWORD PTR -16[rbp]
	movsd	QWORD PTR -8[rbp], xmm0
	movsd	xmm1, QWORD PTR -8[rbp]
	mov	rax, QWORD PTR 24[rbp]
	movsd	xmm0, QWORD PTR [rax]
	addsd	xmm0, xmm1
	mov	rax, QWORD PTR 24[rbp]
	movsd	QWORD PTR [rax], xmm0
	add	DWORD PTR -28[rbp], 1
.L5:
	mov	eax, DWORD PTR -28[rbp]
	cmp	eax, DWORD PTR -44[rbp]
	jle	.L6
	mov	rax, QWORD PTR -64[rbp]
	movsd	xmm0, QWORD PTR [rax]
	pxor	xmm1, xmm1
	cvtsi2sd	xmm1, DWORD PTR -44[rbp]
	divsd	xmm0, xmm1
	mov	rax, QWORD PTR -64[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm1, QWORD PTR [rax]
	movsd	xmm0, QWORD PTR -24[rbp]
	mulsd	xmm0, xmm0
	pxor	xmm3, xmm3
	cvtsi2sd	xmm3, DWORD PTR -44[rbp]
	movapd	xmm2, xmm0
	divsd	xmm2, xmm3
	subsd	xmm1, xmm2
	movapd	xmm0, xmm1
	mov	eax, DWORD PTR -44[rbp]
	sub	eax, 1
	pxor	xmm1, xmm1
	cvtsi2sd	xmm1, eax
	divsd	xmm0, xmm1
	mov	rax, QWORD PTR -80[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR -80[rbp]
	mov	rax, QWORD PTR [rax]
	movq	xmm0, rax
	call	sqrt@PLT
	movq	rax, xmm0
	mov	rdx, QWORD PTR -72[rbp]
	mov	QWORD PTR [rdx], rax
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm0, QWORD PTR [rax]
	pxor	xmm1, xmm1
	ucomisd	xmm0, xmm1
	jp	.L10
	pxor	xmm1, xmm1
	comisd	xmm0, xmm1
	je	.L11
.L10:
	mov	rax, QWORD PTR 16[rbp]
	movsd	xmm0, QWORD PTR [rax]
	pxor	xmm2, xmm2
	cvtsi2sd	xmm2, DWORD PTR -44[rbp]
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm1, QWORD PTR [rax]
	mulsd	xmm2, xmm1
	mov	rax, QWORD PTR -72[rbp]
	movsd	xmm1, QWORD PTR [rax]
	mulsd	xmm1, xmm2
	divsd	xmm0, xmm1
	mov	rax, QWORD PTR 16[rbp]
	movsd	QWORD PTR [rax], xmm0
	mov	rax, QWORD PTR 24[rbp]
	movsd	xmm0, QWORD PTR [rax]
	pxor	xmm2, xmm2
	cvtsi2sd	xmm2, DWORD PTR -44[rbp]
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm1, QWORD PTR [rax]
	mulsd	xmm2, xmm1
	mov	rax, QWORD PTR -80[rbp]
	movsd	xmm1, QWORD PTR [rax]
	mulsd	xmm1, xmm2
	divsd	xmm0, xmm1
	movsd	xmm1, QWORD PTR .LC3[rip]
	subsd	xmm0, xmm1
	mov	rax, QWORD PTR 24[rbp]
	movsd	QWORD PTR [rax], xmm0
	jmp	.L12
.L11:
	lea	rdi, .LC4[rip]
	call	nrerror@PLT
.L12:
	nop
	leave
	.cfi_def_cfa 7, 8
	ret
	.cfi_endproc
.LFE0:
	.size	moment, .-moment
	.section	.rodata
	.align 16
.LC2:
	.long	-1
	.long	2147483647
	.long	0
	.long	0
	.align 8
.LC3:
	.long	0
	.long	1074266112
	.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	"407128.c"
	.text
	.section	.rodata
	.align	2
.LC0:
	.ascii	"n must be at least 2 in moment\000"
	.global	__aeabi_dadd
	.global	__aeabi_i2d
	.global	__aeabi_ddiv
	.global	__aeabi_dsub
	.global	__aeabi_dmul
	.global	__aeabi_dcmpeq
	.align	2
.LC1:
	.ascii	"No skew/kurtosis when variance = 0 (in moment)\000"
	.text
	.align	2
	.global	moment
	.syntax unified
	.arm
	.fpu softvfp
	.type	moment, %function
moment:
	@ args = 16, pretend = 0, frame = 48
	@ frame_needed = 1, uses_anonymous_args = 0
	push	{r4, r5, r6, r7, fp, lr}
	add	fp, sp, #20
	sub	sp, sp, #48
	str	r0, [fp, #-56]
	str	r1, [fp, #-60]
	str	r2, [fp, #-64]
	str	r3, [fp, #-68]
	mov	r2, #0
	mov	r3, #0
	str	r2, [fp, #-44]
	str	r3, [fp, #-40]
	ldr	r3, [fp, #-60]
	cmp	r3, #1
	bgt	.L2
	ldr	r0, .L13
	bl	nrerror
.L2:
	mov	r2, #0
	mov	r3, #0
	str	r2, [fp, #-36]
	str	r3, [fp, #-32]
	mov	r3, #1
	str	r3, [fp, #-48]
	b	.L3
.L4:
	ldr	r3, [fp, #-48]
	lsl	r3, r3, #3
	ldr	r2, [fp, #-56]
	add	r3, r2, r3
	ldmia	r3, {r2-r3}
	sub	r1, fp, #36
	ldmia	r1, {r0-r1}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-36]
	str	r3, [fp, #-32]
	ldr	r3, [fp, #-48]
	add	r3, r3, #1
	str	r3, [fp, #-48]
.L3:
	ldr	r2, [fp, #-48]
	ldr	r3, [fp, #-60]
	cmp	r2, r3
	ble	.L4
	ldr	r0, [fp, #-60]
	bl	__aeabi_i2d
	mov	r2, r0
	mov	r3, r1
	sub	r1, fp, #36
	ldmia	r1, {r0-r1}
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #-64]
	stm	r1, {r2-r3}
	ldr	r1, [fp, #16]
	mov	r2, #0
	mov	r3, #0
	stm	r1, {r2-r3}
	ldr	r3, [fp, #16]
	ldmia	r3, {r2-r3}
	ldr	r1, [fp, #12]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #12]
	ldmia	r3, {r2-r3}
	ldr	r1, [fp, #8]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	ldr	r1, [fp, #-68]
	stm	r1, {r2-r3}
	mov	r3, #1
	str	r3, [fp, #-48]
	b	.L5
.L6:
	ldr	r3, [fp, #-48]
	lsl	r3, r3, #3
	ldr	r2, [fp, #-56]
	add	r3, r2, r3
	ldmia	r3, {r0-r1}
	ldr	r3, [fp, #-64]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-36]
	str	r3, [fp, #-32]
	ldr	r4, [fp, #-36]
	ldr	r3, [fp, #-32]
	bic	r5, r3, #-2147483648
	ldr	r3, [fp, #-68]
	ldmia	r3, {r2-r3}
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #-68]
	stm	r1, {r2-r3}
	sub	r3, fp, #36
	ldmia	r3, {r2-r3}
	sub	r1, fp, #44
	ldmia	r1, {r0-r1}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-44]
	str	r3, [fp, #-40]
	sub	r3, fp, #36
	ldmia	r3, {r2-r3}
	sub	r1, fp, #36
	ldmia	r1, {r0-r1}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-28]
	str	r3, [fp, #-24]
	sub	r1, fp, #28
	ldmia	r1, {r0-r1}
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #8]
	stm	r1, {r2-r3}
	sub	r3, fp, #36
	ldmia	r3, {r2-r3}
	sub	r1, fp, #28
	ldmia	r1, {r0-r1}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-28]
	str	r3, [fp, #-24]
	sub	r1, fp, #28
	ldmia	r1, {r0-r1}
	ldr	r3, [fp, #12]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #12]
	stm	r1, {r2-r3}
	sub	r3, fp, #36
	ldmia	r3, {r2-r3}
	sub	r1, fp, #28
	ldmia	r1, {r0-r1}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	str	r2, [fp, #-28]
	str	r3, [fp, #-24]
	sub	r1, fp, #28
	ldmia	r1, {r0-r1}
	ldr	r3, [fp, #16]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dadd
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #16]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #-48]
	add	r3, r3, #1
	str	r3, [fp, #-48]
.L5:
	ldr	r2, [fp, #-48]
	ldr	r3, [fp, #-60]
	cmp	r2, r3
	ble	.L6
	ldr	r3, [fp, #-68]
	ldmia	r3, {r4-r5}
	ldr	r0, [fp, #-60]
	bl	__aeabi_i2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #-68]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #8]
	ldmia	r3, {r4-r5}
	sub	r3, fp, #44
	ldmia	r3, {r2-r3}
	sub	r1, fp, #44
	ldmia	r1, {r0-r1}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r6, r2
	mov	r7, r3
	ldr	r0, [fp, #-60]
	bl	__aeabi_i2d
	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_dsub
	mov	r2, r0
	mov	r3, r1
	mov	r4, r2
	mov	r5, r3
	ldr	r3, [fp, #-60]
	sub	r3, r3, #1
	mov	r0, r3
	bl	__aeabi_i2d
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #8]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	mov	r0, r2
	mov	r1, r3
	bl	sqrt
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #4]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #8]
	ldmia	r3, {r0-r1}
	mov	r2, #0
	mov	r3, #0
	bl	__aeabi_dcmpeq
	mov	r3, r0
	cmp	r3, #0
	bne	.L11
	ldr	r3, [fp, #12]
	ldmia	r3, {r4-r5}
	ldr	r0, [fp, #-60]
	bl	__aeabi_i2d
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	ldr	r3, [fp, #4]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r4
	mov	r1, r5
	bl	__aeabi_ddiv
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #12]
	stm	r1, {r2-r3}
	ldr	r3, [fp, #16]
	ldmia	r3, {r4-r5}
	ldr	r0, [fp, #-60]
	bl	__aeabi_i2d
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	bl	__aeabi_dmul
	mov	r2, r0
	mov	r3, r1
	mov	r0, r2
	mov	r1, r3
	ldr	r3, [fp, #8]
	ldmia	r3, {r2-r3}
	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	r0, r2
	mov	r1, r3
	mov	r2, #0
	ldr	r3, .L13+4
	bl	__aeabi_dsub
	mov	r2, r0
	mov	r3, r1
	ldr	r1, [fp, #16]
	stm	r1, {r2-r3}
	b	.L12
.L11:
	ldr	r0, .L13+8
	bl	nrerror
.L12:
	nop
	sub	sp, fp, #20
	@ sp needed
	pop	{r4, r5, r6, r7, fp, pc}
.L14:
	.align	2
.L13:
	.word	.LC0
	.word	1074266112
	.word	.LC1
	.size	moment, .-moment
	.ident	"GCC: (Ubuntu 9.4.0-1ubuntu1~20.04.2) 9.4.0"
	.section	.note.GNU-stack,"",%progbits
