how would I split this code into two different files without harming functionality. Everytime I try I tend to get different results than I do when the code is on one file

global compute$\_$variance

extern printf

segment $.$data

output$\_$float db $"\%1\mathrm{.}2$f$",10,0$

debug db "Debug", $10,0$

segment $.$bss

align $64$

storedata resb $832$

segment $.$text

compute$\_$variance:

; Back up components

push rbp

mov rbp$,$ rsp

push rbx

push rcx

push rdx

push rsi

push rdi

push r$8$

push r$9$

push r$10$

push r$11$

push r$12$

push r$13$

push r$14$

push r$15$

pushf

; Save all the floating $-$ point numbers

mov rax, $7$

mov rdx$,0$

xsave $[$storedata$]$

xor r$13,$ r$13$ ; r$13$ keeps track of the index of the loop

mov r$14,$ rdi ; rdi contains the array

mov r$15,$ rsi ; rsi contains the size of the array

pxor xmm$15,$ xmm$15$ ; Clear xmm$15$ to accumulate the sum

pxor xmm$14,$ xmm$14$ ; Clear xmm$14$ to use for intermediate calculations

mean$\_$loop:

; If the array is full, end the loop

cmp r$13,$ r$15$

jge mean$\_$finished

; Add n elements from the array into xmm$15$

addsd xmm$15,[$r$14+$ r$13*8]$

; Increase the index

inc r$13$

; Repeat the loop

jmp mean$\_$loop

mean$\_$finished:

; Calculate the mean by dividing the sum by the count

cvtsi$2$sd xmm$1,$ r$15$ ; Convert count to double

divsd xmm$15,$ xmm$1$ ; Divide sum by count to get mean

movsd xmm$14,$ xmm$15$

pxor xmm$15,$ xmm$15$ ; Copy mean to xmm$14$ for variance calculation

; Reset index for variance calculation

xor r$13,$ r$13$

variance$\_$loop:

cmp r$13,$ r$15$

jge variance$\_$finished

; Compute squared difference from the mean for each element

movsd xmm$0,[$r$14+$ r$13*8]$ ; Load element

subsd xmm$0,$ xmm$14$ ; Subtract mean

mulsd xmm$0,$ xmm$0$ ; Square the difference

addsd xmm$15,$ xmm$0$ ; Add to the sum of squared differences

; Increase the index

inc r$13$

; Repeat the loop

jmp variance$\_$loop

variance$\_$finished:

; Calculate the variance

cvtsi$2$sd xmm$1,$ r$15$ ; Convert count to double

divsd xmm$15,$ xmm$1$ ; Divide sum of squared differences by count

; Save the variance on the stack

sub rsp$,8$

movsd $[$rsp$],$ xmm$15$

; Restore all the floating $-$ point numbers

mov rax, $7$

mov rdx$,0$

xrstor $[$storedata$]$

; Restore the variance and send to manage

movsd xmm$0,[$rsp$]$

add rsp$,8$

; Restore the original values to the GPRs

popf

pop r$15$

pop r$14$

pop r$13$

pop r$12$

pop r$11$

pop r$10$

pop r$9$

pop r$8$

pop rdi

pop rsi

pop rdx

pop rcx

pop rbx

pop rbp

; Clean up

ret