// File: /workspace/github_c_makefile_star2_O2/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:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x828, %rsp            # imm = 0x828
testl %ecx, %ecx
je 0x5a5c
movl %ecx, %r15d
movq %rdi, %rbp
leal -0x1(%r15), %eax
cmpb $0x0, (%rdi,%rax)
je 0x5a63
pushq $-0x1
popq %rbp
jmp 0x5d98
xorl %ebp, %ebp
jmp 0x5d98
movq %rdx, %r12
movq %rsi, %r13
movq %rax, (%rsp)
movl %r15d, %eax
movq %rax, 0x8(%rsp)
leaq (,%rax,4), %rdi
callq 0x76d0
movq %rax, %rbx
movl %r15d, %edi
callq 0x59fc
movq %rax, %r14
testq %rbx, %rbx
sete %al
testq %r14, %r14
sete %cl
orb %al, %cl
jne 0x5ab4
xorl %eax, %eax
cmpq $0x100, %rax            # imm = 0x100
je 0x5abc
andl $0x0, 0x24(%rsp,%rax,8)
incq %rax
jmp 0x5aa2
pushq $-0x2
popq %rbp
jmp 0x5d88
xorl %eax, %eax
cmpq %rax, 0x8(%rsp)
je 0x5ad3
movzbl (%rbp,%rax), %ecx
incl 0x24(%rsp,%rcx,8)
incq %rax
jmp 0x5abe
cmpl $0x1, 0x24(%rsp)
jne 0x5dac
xorl %eax, %eax
xorl %ecx, %ecx
cmpq $0x100, %rax            # imm = 0x100
je 0x5b03
movl 0x24(%rsp,%rax,8), %edx
testl %edx, %edx
je 0x5afe
movl %ecx, 0x20(%rsp,%rax,8)
movl %ecx, %esi
movl %edx, (%r14,%rsi,4)
addl %edx, %ecx
incq %rax
jmp 0x5ae2
movq (%rsp), %rax
cmpl %r15d, %eax
jae 0x5b38
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, (%rbx,%rcx,4)
movl %edx, (%r12,%rcx,4)
movl %eax, (%r13,%rdx,4)
decl %eax
jmp 0x5b07
movq (%rsp), %rax
testl %eax, %eax
je 0x5d0d
movl %eax, %ecx
movq %rcx, 0x18(%rsp)
movl (%r13,%rcx,4), %ecx
movl (%rbx,%rcx,4), %edx
andl $0x0, (%r14,%rdx,4)
movl %eax, %esi
cmpl %edx, %esi
jb 0x5b93
movl %esi, %ecx
movl (%r13,%rcx,4), %edi
leal -0x1(%rdi), %r8d
cmpl %r15d, %r8d
jae 0x5b8b
movl %r8d, %r9d
movl (%r12,%r9,4), %ecx
cmpl %eax, %ecx
jbe 0x5b7f
movl (%rbx,%r9,4), %r8d
jmp 0x5b69
cmpl %edx, %ecx
jb 0x5b8b
movl $0x1, (%r14,%rcx,4)
movl %r8d, (%rbx,%rdi,4)
decl %esi
jmp 0x5b5a
xorl %r8d, %r8d
movl %edx, %esi
cmpl %eax, %esi
ja 0x5bc9
movl %esi, %r9d
movl (%r13,%r9,4), %edi
movl %eax, (%r12,%rdi,4)
cmpl $0x0, (%r14,%r9,4)
jne 0x5bc5
leal (%r8,%rdx), %ecx
movl (%r13,%rcx,4), %r10d
movl %r10d, (%r13,%r9,4)
incl %r8d
movl %edi, (%r13,%rcx,4)
incl %esi
jmp 0x5b98
addl %edx, %r8d
xorl %esi, %esi
movl %edx, %edi
movl %r8d, %r9d
leal -0x1(%r9), %r10d
cmpl %edi, %r10d
jb 0x5c29
movl %r10d, %r11d
movl (%r13,%r11,4), %ecx
movl (%rbx,%rcx,4), %ebp
cmpl %r15d, %ebp
jae 0x5c02
cmpl %edx, (%r12,%rbp,4)
jb 0x5c12
movl (%rbx,%rbp,4), %r11d
movl %r11d, (%rbx,%rcx,4)
movl %ecx, (%rbx,%rbp,4)
decl %r10d
jmp 0x5bd7
movl %edi, %ecx
incl %edi
movl (%r13,%rcx,4), %ecx
movl %ecx, (%r13,%r11,4)
jmp 0x5bd7
decl %r8d
movl (%r13,%r8,4), %ecx
movl %ecx, (%r13,%r11,4)
movl %ebp, (%r13,%r8,4)
decl %r10d
jmp 0x5bd7
subl %r8d, %r9d
jbe 0x5c37
movl %r8d, %ecx
movl %r9d, (%r14,%rcx,4)
incl %esi
cmpl %r8d, %edi
jb 0x5bd0
jmp 0x5cda
movl %esi, 0x14(%rsp)
movl %edi, %r8d
movl (%r14,%r8,4), %ecx
leal (%rcx,%rdi), %r10d
decl %r10d
movl %ecx, %r9d
addl %edi, %r9d
cmpl %edi, %r10d
jb 0x5c9a
movl %r10d, %ecx
movl (%r13,%rcx,4), %ecx
movl (%rbx,%rcx,4), %r11d
movl (%r14,%r11,4), %ebp
leal -0x1(%rbp), %esi
movl %esi, (%r14,%r11,4)
leal -0x1(%rbp,%r11), %esi
movl (%r13,%rsi,4), %r11d
movl (%r12,%rcx,4), %ebp
movl %r11d, (%r13,%rbp,4)
movl %ebp, (%r12,%r11,4)
movl %ecx, (%r13,%rsi,4)
movl %esi, (%r12,%rcx,4)
decl %r10d
jmp 0x5c59
movl %r9d, %edi
movq %r8, %r10
cmpq %rdi, %r10
jae 0x5cbd
movl (%r13,%r10,4), %ecx
movl (%rbx,%rcx,4), %esi
movl (%r14,%rsi,4), %r11d
addl %esi, %r11d
movl %r11d, (%rbx,%rcx,4)
incq %r10
jmp 0x5ca0
movl 0x14(%rsp), %esi
cmpq %rdi, %r8
jae 0x5cd7
movl (%r13,%r8,4), %ecx
movl (%rbx,%rcx,4), %ecx
incl (%r14,%rcx,4)
incq %r8
jmp 0x5cc1
movl %r9d, %edi
subl $0x1, %esi
jae 0x5c41
jmp 0x5cf7
movl %eax, %ecx
movl (%r13,%rcx,4), %ecx
movl (%rbx,%rcx,4), %esi
movl (%rbx,%rsi,4), %esi
movl %esi, (%rbx,%rcx,4)
decl %eax
cmpl %edi, %eax
jae 0x5ce5
movq 0x18(%rsp), %rax
movl %edx, (%r13,%rax,4)
leal -0x1(%rdx), %eax
jmp 0x5b3c
movq (%rsp), %rax
movl %eax, (%r13)
xorl %eax, %eax
cmpq 0x8(%rsp), %rax
je 0x5d56
movl (%r13,%rax,4), %ecx
decl %ecx
cmpl %r15d, %ecx
jae 0x5d51
movl %ecx, %edx
movl (%r12,%rdx,4), %esi
cmpl %r15d, %esi
je 0x5d51
movl (%r13,%rsi,4), %edi
leal 0x1(%rdi), %r8d
movl %r8d, (%r13,%rsi,4)
movl %ecx, (%r13,%rdi,4)
movl %r15d, (%r12,%rdx,4)
movl (%rbx,%rdx,4), %ecx
jmp 0x5d25
incq %rax
jmp 0x5d17
xorl %eax, %eax
movq (%rsp), %rdx
cmpq %rax, 0x8(%rsp)
je 0x5d80
movl $0x1, (%r12,%rax,4)
incq %rax
jmp 0x5d5c
movl (%rbx,%rdx,4), %eax
cmpl %r15d, %eax
jae 0x5d80
movl (%r12,%rdx,4), %ecx
addl %ecx, (%r12,%rax,4)
subq $0x1, %rdx
jae 0x5d70
xorl %ebp, %ebp
movq %rbx, %rdi
callq 0x77a6
movq %r14, %rdi
callq 0x77a6
movl %ebp, %eax
addq $0x828, %rsp            # imm = 0x828
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq 0x6245(%rip), %rax      # 0xbff8
movq (%rax), %rcx
leaq 0x29f0(%rip), %rdi      # 0x87ad
pushq $0x15
popq %rsi
pushq $0x1
popq %rdx
callq 0x2280
pushq $-0x1
jmp 0x5ab6