// File: /workspace/github_c_makefile_star2_O3/felipelouza[P]lyndon-array/external/gsaca_cl/gsaca.c
// Function: gsaca_unsigned_cl_la
int gsaca_unsigned_cl_la(const unsigned char *S, unsigned int *SA, unsigned int *LA, unsigned int n) {
	unsigned int *ISA, *PREV;
	void *GSIZE;
	unsigned int i,j;
	unsigned int gstart, gend;
	unsigned int s, p, sr;
	unsigned int tmp, gstarttmp, gendtmp;
	unsigned int C[2*C_SIZE]; //counts and cumulative counts

	if (n == 0)		return 0;
	if (S[n-1] != '\0')	return -1;

	//set up needed structures
	ISA = LA;
	PREV = (unsigned int *)malloc( n * sizeof(unsigned int) );
	GSIZE = gsize_calloc( n );
	if (!PREV || !GSIZE)	{
		free(PREV); gsize_free(GSIZE);
		return -2;
	}

	//// PHASE 1: pre-sort suffixes ////
	//build initial group structure
	for (i = 0; i < C_SIZE; ++i)	C_CNT(i) = 0;
	for (i = 0; i < n; ++i) ++C_CNT(S[i]); //count characters

	if (C_CNT(0) != 1) { //more than 1 nullbyte
		fprintf(stderr, "more than 1 nullbyte\n");
		free(PREV); gsize_free(GSIZE);
		return -1;
	}
	j = 0;
	for (i = 0; i < C_SIZE; ++i) { //build cumulative counts and set up GSIZE
		if (C_CNT(i) > 0) {
			C_CUM(i) = j;
			gsize_set(GSIZE, j, C_CNT(i));
			j += C_CNT(i);
		}
	}
	for (i = n-1; i < n; --i) { //set up ISA, GLINK and SA
		gstart = C_CUM(S[i]);
		sr = gstart + --C_CNT(S[i]);
		GLINK[i] = gstart;
		ISA[i] = sr;
		SA[sr] = i;
	}

	//process groups from highest to lowest
	for (gend = n-1; gend > 0; gend = gstarttmp-1) {
		gstart = GLINK[SA[gend]];
		gstarttmp = gstart;	gendtmp = gend;

		//clear GSIZE group size for marking
		gsize_clear(GSIZE, gstart);

		//compute prev - pointers and mark suffixes of own group that
		//have a prev-pointer of own group pointing to them
		for (i = gend; i >= gstart; --i) {
			s = SA[i]; //use prev - pointers from already used groups
			for (p = s-1; p < n; p = PREV[p]) {
				if (ISA[p] <= gend) {
					if (ISA[p] >= gstart) {
						gsize_set(GSIZE, ISA[p], 1); //mark ISA[p]
					}
					break;
				}
			}
			PREV[s] = p;
		}

		//set GENDLINK of all suffixes for phase 2 and move unmarked
		// suffixes to the front of the actual group
		j = 0;
		for (i = gstart; i <= gend; ++i) {
			s = SA[i];
			GENDLINK[s] = gend;
			if (gsize_get(GSIZE, i) == 0) { //i is not marked
				//swap SA[i] and SA[gstart+j] TODO: changed, now a swap is performed
				SA[i] = SA[gstart+j];
				SA[gstart+(j++)] = s;
			}
		}

		//order the suffixes according on how much suffixes of same
		//group are jumped by them
		gend = gstart + j; //exclusive bound by now
		j = 0;
		do {
			i = gend-1; sr = gend;
			while (i >= gstart) {
				s = SA[i];
				p = PREV[s];
				if (p < n) {
					if (ISA[p] < gstarttmp) { //p is in a lex. smaller group
						SA[i--] = SA[--gend];
						SA[gend] = p; //push prev to back
					} else { //p is in same group
						PREV[s] = PREV[p];
						PREV[p] = s; //set a reference from p to s, as s's prev pointer will be correct later TODO changed, could be ignored before
						--i;
					}
				} else { //prev points to nothing
					SA[i] = SA[gstart++]; //remove entry
				}
			}
			//write number of suffixes written to end on stack using GSIZE
			if (gend < sr) {
				gsize_set(GSIZE, gend, sr-gend);
				++j; //also, count number of splitted groups
			}
		} while (gstart < gend);

		//rearrange previous suffixes stored in other groups
		while (j--) {
			gend = gstart + gsize_get(GSIZE, gstart);
			//decrement group count of previous group suffixes, and move them to back
			for (i = gend-1; i >= gstart; --i) {
				p = SA[i];
				sr = GLINK[p];
				sr += gsize_dec_get(GSIZE, sr);
				//move p to back by exchanging it with last suffix s of group
				s = SA[sr];
				tmp = ISA[p];
				SA[tmp] = s;	ISA[s] = tmp;
				SA[sr] = p;	ISA[p] = sr;				
			}
			//set new GLINK for moved suffixes
			for (i = gstart; i < gend; ++i) {
				p = SA[i];
				sr = GLINK[p];
				sr += gsize_get(GSIZE, sr);
				GLINK[p] = sr;
			}
			//set up GSIZE for newly created groups
			for (i = gstart; i < gend; ++i) {
				p = SA[i];
				sr = GLINK[p];
				gsize_inc(GSIZE, sr);
			}
			gstart = gend;
		}
		
		//set remaining prev pointers of own group to correct value TODO: new, was ignored before
		for (i = gendtmp; i >= gstart; --i) {
			s = SA[i];
			PREV[s] = PREV[PREV[s]]; //use stored references to get correct pointers
		}

		//prepare current group for phase 2
		SA[gendtmp] = gstarttmp; //counter where to place next entry TODO: removed because not necessary
	}

	//// PHASE 2: sort suffixes finally ////
	SA[0] = n-1;
	for (i = 0; i < n; i++) {
		s = SA[i]-1;
		while (s < n) {
			sr = GENDLINK[s];
			if (sr == n) //suffix already placed to SA, stop
				break;
			sr = SA[sr]++; //get position where to write s
			SA[sr] = s;
			//mark that suffix is placed in SA already
			GENDLINK[s] = n; 
			/* NOTE: Using annother technique
			   to indicate if a suffix is placed to SA already,
			   ISA could be computed too. A possible way would be to
			   set PREV[s] = s, and check that in the above if. 
			*/
			s = PREV[s]; //process next suffix
		}
	}

	//// PHASE 3: compute Lyndon array ////
	//compute for each suffix S[i..n] the length of its group context
	for (i = 0; i < n; i++) { //set all context lengths to 1
		LA[i] = 1;
	}

	//use prev pointers to compute context lengths
	for (i = n-1; i-- > 0; ) { //run from back to front
		p = PREV[i];
		if (p < n) { //ensure previous suffix exists
			LA[p] += LA[i]; //add own suffix length to that of p
		}
	}
	

	free(PREV); gsize_free(GSIZE);

	return 0;	
}

// Assembly:
testl %ecx, %ecx
je 0x6533
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x828, %rsp            # imm = 0x828
movl %ecx, %ebx
movq %rdi, %rbp
leal -0x1(%rbx), %ecx
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
cmpb $0x0, (%rdi,%rcx)
jne 0x6521
movq %rdx, %r14
movq %rsi, %r13
movq %rcx, (%rsp)
movl %ebx, %r12d
leaq (,%r12,4), %rdi
callq 0x855a
movq %rax, %r15
movl $0x4, %esi
movq %r12, 0x10(%rsp)
movq %r12, %rdi
callq 0x86b3
movq %rax, %r12
testq %r15, %r15
sete %al
testq %r12, %r12
sete %cl
movl $0xfffffffe, %esi       # imm = 0xFFFFFFFE
orb %al, %cl
jne 0x650d
xorl %eax, %eax
movl $0x0, 0x24(%rsp,%rax,8)
incq %rax
cmpq $0x100, %rax            # imm = 0x100
jne 0x6118
xorl %eax, %eax
movq 0x10(%rsp), %r9
movzbl (%rbp,%rax), %ecx
incl 0x24(%rsp,%rcx,8)
incq %rax
cmpq %rax, %r9
jne 0x6132
cmpl $0x1, 0x24(%rsp)
jne 0x6536
xorl %eax, %eax
xorl %ecx, %ecx
movl 0x24(%rsp,%rax,8), %edx
testl %edx, %edx
je 0x6166
movl %ecx, 0x20(%rsp,%rax,8)
movl %ecx, %esi
movl %edx, (%r12,%rsi,4)
addl %edx, %ecx
incq %rax
cmpq $0x100, %rax            # imm = 0x100
jne 0x6152
movq (%rsp), %rax
movl %eax, %ecx
movzbl (%rbp,%rcx), %edx
movl 0x20(%rsp,%rdx,8), %esi
movl 0x24(%rsp,%rdx,8), %edi
leal -0x1(%rdi), %r8d
movl %r8d, 0x24(%rsp,%rdx,8)
leal -0x1(%rdi,%rsi), %edx
movl %esi, (%r15,%rcx,4)
movl %edx, (%r14,%rcx,4)
movl %eax, (%r13,%rdx,4)
decl %eax
cmpl %ebx, %eax
jb 0x6175
cmpl $0x0, (%rsp)
je 0x63ad
movq (%rsp), %rax
movl %eax, %ecx
movq %rcx, 0x18(%rsp)
movl (%r13,%rcx,4), %ecx
movl (%r15,%rcx,4), %edx
movl $0x0, (%r12,%rdx,4)
xorl %esi, %esi
movl %eax, %edi
movl $0x0, %r8d
cmpl %edx, %eax
jb 0x624c
movl %edi, %ecx
movl (%r13,%rcx,4), %r8d
leal -0x1(%r8), %r9d
cmpl %ebx, %r9d
jae 0x6205
movl %r9d, %r10d
movl (%r14,%r10,4), %ecx
cmpl %eax, %ecx
jbe 0x61f9
movl (%r15,%r10,4), %r9d
jmp 0x61e3
cmpl %edx, %ecx
jb 0x6205
movl $0x1, (%r12,%rcx,4)
movl %r9d, (%r15,%r8,4)
decl %edi
cmpl %edx, %edi
jae 0x61d8
movl $0x0, %r8d
cmpl %edx, %eax
jb 0x624c
xorl %r8d, %r8d
movl %edx, %edi
movl %edi, %edi
movl (%r13,%rdi,4), %ecx
movl %eax, (%r14,%rcx,4)
cmpl $0x0, (%r12,%rdi,4)
jne 0x6246
leal (%r8,%rdx), %r9d
movl (%r13,%r9,4), %r10d
movl %r10d, (%r13,%rdi,4)
incl %r8d
movl %ecx, (%r13,%r9,4)
incl %edi
cmpl %eax, %edi
jbe 0x621e
addl %edx, %r8d
movl %edx, %edi
leal -0x1(%r8), %r10d
cmpl %edi, %r10d
jae 0x625f
movl %r8d, %r9d
jmp 0x62be
movl %r8d, %r9d
movl %r10d, %r11d
movl (%r13,%r11,4), %ecx
movl (%r15,%rcx,4), %ebp
cmpl %ebx, %ebp
jae 0x628f
decl %r10d
cmpl %edx, (%r14,%rbp,4)
jae 0x629f
decl %r9d
movl (%r13,%r9,4), %ecx
movl %ecx, (%r13,%r11,4)
movl %ebp, (%r13,%r9,4)
jmp 0x62ab
movl %edi, %ecx
incl %edi
movl (%r13,%rcx,4), %ecx
movl %ecx, (%r13,%r11,4)
jmp 0x62ab
movl (%r15,%rbp,4), %r11d
movl %r11d, (%r15,%rcx,4)
movl %ecx, (%r15,%rbp,4)
cmpl %edi, %r10d
jae 0x6262
subl %r9d, %r8d
jbe 0x62be
movl %r9d, %ecx
movl %r8d, (%r12,%rcx,4)
incl %esi
movl %r9d, %r8d
cmpl %r9d, %edi
jb 0x6251
testl %esi, %esi
je 0x6377
movl %esi, 0xc(%rsp)
movl %edi, %r9d
movl %edi, %r8d
movl (%r12,%r8,4), %edi
leal (%rdi,%r9), %r10d
decl %r10d
cmpl %r9d, %r10d
jb 0x6327
movl %r10d, %ecx
movl (%r13,%rcx,4), %ecx
movl (%r15,%rcx,4), %r11d
movl (%r12,%r11,4), %ebp
leal -0x1(%rbp), %esi
movl %esi, (%r12,%r11,4)
leal -0x1(%rbp,%r11), %esi
movl (%r13,%rsi,4), %r11d
movl (%r14,%rcx,4), %ebp
movl %r11d, (%r13,%rbp,4)
movl %ebp, (%r14,%r11,4)
movl %ecx, (%r13,%rsi,4)
movl %esi, (%r14,%rcx,4)
decl %r10d
cmpl %r8d, %r10d
jae 0x62e8
addl %r9d, %edi
cmpl %edi, %r8d
jae 0x636b
movl %edi, %r9d
movq %r8, %r10
movl (%r13,%r10,4), %ecx
movl (%r15,%rcx,4), %esi
movl (%r12,%rsi,4), %r11d
addl %esi, %r11d
movl %r11d, (%r15,%rcx,4)
incq %r10
cmpq %r10, %r9
jne 0x6335
cmpl %edi, %r8d
jae 0x636b
movl (%r13,%r8,4), %ecx
movl (%r15,%rcx,4), %ecx
incl (%r12,%rcx,4)
incq %r8
cmpq %r8, %r9
jne 0x6356
movl 0xc(%rsp), %esi
decl %esi
jne 0x62ce
cmpl %edi, %eax
movq 0x10(%rsp), %r9
jb 0x6399
movl %eax, %ecx
movl (%r13,%rcx,4), %ecx
movl (%r15,%rcx,4), %esi
movl (%r15,%rsi,4), %esi
movl %esi, (%r15,%rcx,4)
decl %eax
cmpl %edi, %eax
jae 0x6380
movq 0x18(%rsp), %rax
movl %edx, (%r13,%rax,4)
decl %edx
movl %edx, %eax
jne 0x61b2
movq (%rsp), %rax
movl %eax, (%r13)
xorl %eax, %eax
movl (%r13,%rax,4), %ecx
decl %ecx
cmpl %ebx, %ecx
jae 0x63e9
movl %ecx, %edx
movl (%r14,%rdx,4), %esi
cmpl %ebx, %esi
je 0x63e9
movl (%r13,%rsi,4), %edi
leal 0x1(%rdi), %r8d
movl %r8d, (%r13,%rsi,4)
movl %ecx, (%r13,%rdi,4)
movl %ebx, (%r14,%rdx,4)
movl (%r15,%rdx,4), %ecx
jmp 0x63be
incq %rax
cmpq %r9, %rax
jne 0x63b7
leaq 0x3(%r9), %rax
andq $-0x4, %rax
leaq -0x1(%r9), %rcx
movq %rcx, %xmm0
pshufd $0x44, %xmm0, %xmm0     # xmm0 = xmm0[0,1,0,1]
movdqa 0x33a1(%rip), %xmm1     # 0x97b0
movdqa 0x33a9(%rip), %xmm2     # 0x97c0
xorl %ecx, %ecx
movdqa 0x33af(%rip), %xmm3     # 0x97d0
pxor %xmm3, %xmm0
pcmpeqd %xmm4, %xmm4
movdqa 0x33af(%rip), %xmm5     # 0x97e0
movdqa %xmm2, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
movd %xmm6, %edx
notl %edx
testb $0x1, %dl
je 0x646d
movl $0x1, (%r14,%rcx,4)
pxor %xmm4, %xmm6
pextrw $0x4, %xmm6, %edx
testb $0x1, %dl
je 0x6484
movl $0x1, 0x4(%r14,%rcx,4)
movdqa %xmm1, %xmm6
pxor %xmm3, %xmm6
movdqa %xmm6, %xmm7
pcmpgtd %xmm0, %xmm7
pcmpeqd %xmm0, %xmm6
pshufd $0xf5, %xmm6, %xmm8     # xmm8 = xmm6[1,1,3,3]
pand %xmm7, %xmm8
pshufd $0xf5, %xmm7, %xmm6     # xmm6 = xmm7[1,1,3,3]
por %xmm8, %xmm6
pxor %xmm4, %xmm6
pextrw $0x0, %xmm6, %edx
testb $0x1, %dl
je 0x64c4
movl $0x1, 0x8(%r14,%rcx,4)
pextrw $0x4, %xmm6, %edx
testb $0x1, %dl
je 0x64d7
movl $0x1, 0xc(%r14,%rcx,4)
addq $0x4, %rcx
paddq %xmm5, %xmm2
paddq %xmm5, %xmm1
cmpq %rcx, %rax
jne 0x6431
xorl %esi, %esi
cmpl $0x0, (%rsp)
je 0x650d
leal -0x2(%rbx), %eax
movl (%r15,%rax,4), %ecx
cmpl %ebx, %ecx
jae 0x6507
movl (%r14,%rax,4), %edx
addl %edx, (%r14,%rcx,4)
addq $-0x1, %rax
jb 0x64f7
movq %r15, %rdi
movl %esi, %ebx
callq 0x8611
movq %r12, %rdi
callq 0x8611
movl %ebx, %eax
addq $0x828, %rsp            # imm = 0x828
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
xorl %eax, %eax
retq
callq 0x2f26
movl $0xffffffff, %esi       # imm = 0xFFFFFFFF
jmp 0x650d