// File: /workspace/github_c_makefile_star2_O1/felipelouza[P]lyndon-array/external/gsaca_cl/gsaca.c
// Function: gsaca_unsigned_cl
int gsaca_unsigned_cl(const unsigned char *S, unsigned int *SA, 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 = (unsigned int *)malloc( n * sizeof(unsigned int) );
	PREV = (unsigned int *)malloc( n * sizeof(unsigned int) );
	GSIZE = gsize_calloc( n );
	if (!ISA || !PREV || !GSIZE)	{
		free(ISA); 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(ISA); 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;
		}
		
		//prepare current group for phase 2
		//SA[gendtmp] = gstarttmp; //counter where to place next entry TODO: removed because not necessary
		
		//set remaining prev pointers of own group to correct value TODO: new, was ignored before
		while (gendtmp >= gstart) {
			s = SA[gendtmp--];
			PREV[s] = PREV[PREV[s]]; //use stored references to get correct pointers
		}
	}

	//// FOR FELIPE ////
	//compute for each suffix S[i..n] the length of its group context
	unsigned int *CLEN = SA; //context length of each suffix
	for (i = 0; i < n; i++) { //set all context lengths to 1
		CLEN[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
			CLEN[p] += CLEN[i]; //add own suffix length to that of p
		}
	}
	

	free(ISA); free(PREV); gsize_free(GSIZE);

	return 0;	
}

// Assembly:
testl %edx, %edx
je 0x5f53
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x818, %rsp            # imm = 0x818
movl %edx, %ebx
movq %rdi, %rbp
leal -0x1(%rbx), %ecx
movl $0xffffffff, %eax       # imm = 0xFFFFFFFF
cmpb $0x0, (%rdi,%rcx)
jne 0x5f41
movq %rsi, %r14
movq %rcx, (%rsp)
movl %ebx, %r13d
leaq (,%r13,4), %r12
movq %r12, %rdi
callq 0x7f0e
movq %rax, %r15
movq %r12, %rdi
callq 0x7f0e
movq %rax, %r12
movl $0x4, %esi
movq %r13, 0x8(%rsp)
movq %r13, %rdi
callq 0x8067
movq %rax, %r13
testq %r15, %r15
sete %al
testq %r12, %r12
sete %cl
orb %al, %cl
testq %r13, %r13
sete %al
movl $0xfffffffe, %esi       # imm = 0xFFFFFFFE
orb %cl, %al
jne 0x5f25
xorl %eax, %eax
movl $0x0, 0x14(%rsp,%rax,8)
incq %rax
cmpq $0x100, %rax            # imm = 0x100
jne 0x5c88
xorl %eax, %eax
movq 0x8(%rsp), %rdx
movzbl (%rbp,%rax), %ecx
incl 0x14(%rsp,%rcx,8)
incq %rax
cmpq %rax, %rdx
jne 0x5ca2
cmpl $0x1, 0x14(%rsp)
jne 0x5f56
xorl %eax, %eax
xorl %ecx, %ecx
movl 0x14(%rsp,%rax,8), %edx
testl %edx, %edx
je 0x5cd7
movl %ecx, 0x10(%rsp,%rax,8)
movl %ecx, %esi
movl %edx, (%r13,%rsi,4)
addl %edx, %ecx
incq %rax
cmpq $0x100, %rax            # imm = 0x100
jne 0x5cc2
movq (%rsp), %rax
movl %eax, %ecx
movzbl (%rbp,%rcx), %edx
movl 0x10(%rsp,%rdx,8), %esi
movl 0x14(%rsp,%rdx,8), %edi
leal -0x1(%rdi), %r8d
movl %r8d, 0x14(%rsp,%rdx,8)
leal -0x1(%rdi,%rsi), %edx
movl %esi, (%r12,%rcx,4)
movl %edx, (%r15,%rcx,4)
movl %eax, (%r14,%rdx,4)
decl %eax
cmpl %ebx, %eax
jb 0x5ce6
cmpl $0x0, (%rsp)
je 0x5eea
movq (%rsp), %rax
movl %eax, %ecx
movl (%r14,%rcx,4), %ecx
movl (%r12,%rcx,4), %ecx
movl $0x0, (%r13,%rcx,4)
movl %eax, %edx
cmpl %ecx, %eax
jb 0x5d72
movl %edx, %esi
movl (%r14,%rsi,4), %esi
leal -0x1(%rsi), %edi
cmpl %ebx, %edi
jae 0x5d68
movl %edi, %r9d
movl (%r15,%r9,4), %r8d
cmpl %eax, %r8d
jbe 0x5d5a
movl (%r12,%r9,4), %edi
jmp 0x5d44
cmpl %ecx, %r8d
jb 0x5d68
movl $0x1, (%r13,%r8,4)
movl %edi, (%r12,%rsi,4)
decl %edx
cmpl %ecx, %edx
jae 0x5d3b
xorl %edx, %edx
movl $0x0, %edi
cmpl %eax, %ecx
ja 0x5dab
xorl %edi, %edi
movl %ecx, %esi
movl %esi, %esi
movl (%r14,%rsi,4), %r8d
movl %eax, (%r15,%r8,4)
cmpl $0x0, (%r13,%rsi,4)
jne 0x5da5
leal (%rdi,%rcx), %r9d
movl (%r14,%r9,4), %r10d
movl %r10d, (%r14,%rsi,4)
incl %edi
movl %r8d, (%r14,%r9,4)
incl %esi
cmpl %eax, %esi
jbe 0x5d81
addl %ecx, %edi
movl %ecx, %esi
leal -0x1(%rdi), %r9d
movl %edi, %r8d
cmpl %esi, %r9d
jb 0x5e02
movl %r9d, %r10d
movl (%r14,%r10,4), %ebp
movl (%r12,%rbp,4), %r11d
cmpl %ebx, %r11d
jae 0x5de5
decl %r9d
cmpl %ecx, (%r15,%r11,4)
jae 0x5df4
decl %r8d
movl (%r14,%r8,4), %ebp
movl %ebp, (%r14,%r10,4)
movl %r11d, (%r14,%r8,4)
jmp 0x5db6
movl %esi, %r11d
incl %esi
movl (%r14,%r11,4), %r11d
movl %r11d, (%r14,%r10,4)
jmp 0x5db6
movl (%r12,%r11,4), %r10d
movl %r10d, (%r12,%rbp,4)
movl %ebp, (%r12,%r11,4)
jmp 0x5db6
subl %r8d, %edi
jbe 0x5e11
movl %r8d, %r9d
movl %edi, (%r13,%r9,4)
incl %edx
movl %r8d, %edi
cmpl %r8d, %esi
jb 0x5daf
testl %edx, %edx
je 0x5edc
movl %edx, 0x8(%rsp)
movl %esi, %r8d
movl %esi, %edi
movl (%r13,%rdi,4), %esi
leal (%rsi,%r8), %r9d
decl %r9d
cmpl %r8d, %r9d
jb 0x5e78
movl %r9d, %r10d
movl (%r14,%r10,4), %r10d
movl (%r12,%r10,4), %r11d
movl (%r13,%r11,4), %ebp
leal -0x1(%rbp), %edx
movl %edx, (%r13,%r11,4)
leal -0x1(%rbp,%r11), %edx
movl (%r14,%rdx,4), %r11d
movl (%r15,%r10,4), %ebp
movl %r11d, (%r14,%rbp,4)
movl %ebp, (%r15,%r11,4)
movl %r10d, (%r14,%rdx,4)
movl %edx, (%r15,%r10,4)
decl %r9d
cmpl %edi, %r9d
jae 0x5e3b
addl %r8d, %esi
movl %esi, %r8d
cmpl %esi, %edi
jae 0x5ea1
movq %rdi, %r9
movl (%r14,%r9,4), %edx
movl (%r12,%rdx,4), %r10d
movl (%r13,%r10,4), %r11d
addl %r10d, %r11d
movl %r11d, (%r12,%rdx,4)
incq %r9
cmpq %r9, %r8
jne 0x5e85
cmpl %esi, %edi
jae 0x5eba
movl (%r14,%rdi,4), %edx
movl (%r12,%rdx,4), %edx
incl (%r13,%rdx,4)
incq %rdi
cmpq %rdi, %r8
jne 0x5ea5
movl 0x8(%rsp), %edx
decl %edx
jne 0x5e21
jmp 0x5edc
movl %eax, %edx
decl %eax
movl (%r14,%rdx,4), %edx
movl (%r12,%rdx,4), %edi
movl (%r12,%rdi,4), %edi
movl %edi, (%r12,%rdx,4)
cmpl %esi, %eax
jae 0x5ec8
decl %ecx
movl %ecx, %eax
jne 0x5d22
cmpl $0x1, %ebx
movl %ebx, %eax
adcl $0x0, %eax
xorl %ecx, %ecx
movl $0x1, (%r14,%rcx,4)
incq %rcx
cmpq %rcx, %rax
jne 0x5ef4
xorl %esi, %esi
cmpl $0x0, (%rsp)
je 0x5f25
leal -0x2(%rbx), %eax
movl (%r12,%rax,4), %ecx
cmpl %ebx, %ecx
jae 0x5f1f
movl (%r14,%rax,4), %edx
addl %edx, (%r14,%rcx,4)
addq $-0x1, %rax
jb 0x5f0f
movq %r15, %rdi
movl %esi, %ebx
callq 0x7fc5
movq %r12, %rdi
callq 0x7fc5
movq %r13, %rdi
callq 0x7fc5
movl %ebx, %eax
addq $0x818, %rsp            # imm = 0x818
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
xorl %eax, %eax
retq
callq 0x2ea3
movl $0xffffffff, %esi       # imm = 0xFFFFFFFF
jmp 0x5f25