// File: /workspace/github_cpp_makefile_star2_O1/ratulSharker[P]Havec/src/vector_with_rehash_8_hash_with_N_detection_5_bytes.cpp
// Function: storeKhmerInMemBlock5Byte(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char>>, khmerInfo5Byte*, unsigned long long, unsigned long long, char, unsigned int)
int storeKhmerInMemBlock5Byte(string sample,struct khmerInfo5Byte * memBlock,unsigned long long blockSize, unsigned long long khmer_h, char next_neucleotide, unsigned int whichHashFunc){

	unsigned long long index;
	unsigned long long qutioned; //this value must be under 12bit
	char neighbour_musk = 0xF0;
	char *p;

	//find the index value in the memBlock
	index = khmer_h % blockSize;
	
	qutioned = khmer_h / blockSize;//this value must be under 10bit

	//std::cout<<"data :: "<<(int)memBlock[index].data[1]<<std::endl;
	//cout<<"sample khmer "<<sample<<" kmer_h "<<khmer_h<<" hash func no "<<whichHashFunc<<endl;

	//now check the specified index is neighbourless or not (which is equivalent to -> is the space occupied)
	if((memBlock[index].data[4] & neighbour_musk) == 0){
		//no one is occupuying here put khmer_h here

		//putting the qutioned
		p = (char*) &qutioned;
		memBlock[index].data[0] = *p;
		p++;
		memBlock[index].data[1] = *p;
		p++;
		memBlock[index].data[2] = *p;
		p++;
		memBlock[index].data[3] = *p;
		p++;
		memBlock[index].data[4] = 0x01 & (*p);

		//putting "which hash value is stored"
		memBlock[index].data[4] = memBlock[index].data[4] | (whichHashFunc << 1);
		
		/*memBlock[index].data[0] = qutioned & 0xFF;
		memBlock[index].data[1] = (qutioned >> 8) & 0x0F;
		*/

		//std::cout<<"putting in the block array"<<std::endl;

		switch(toupper(next_neucleotide)){//putting the neighbour
		case 'A'://1000
			memBlock[index].data[4] = memBlock[index].data[4] | 0x80;
			break;
		case 'T'://0100
			memBlock[index].data[4] = memBlock[index].data[4] | 0x40;
			break;
		case 'C'://0010
			memBlock[index].data[4] = memBlock[index].data[4] | 0x20;
			break;
		case 'G'://0001
			memBlock[index].data[4] = memBlock[index].data[4] | 0x10;
			break;
		}
		//std::cout<<"index "<<index<<" quotioned :: "<< qutioned<< " data[0]=" <<(int)memBlock[index].data[0]<<" data[1]="<<(int)memBlock[index].data[1]<<std::endl;

		
		
		return MEM_BLOCK; //khmer putted in the mem_block so return mem_block
	}
	else if(((memBlock[index].data[4] >> 1) & 0x07) == whichHashFunc && //hash value matched
			(memBlock[index].data[0] == (qutioned & 0xFF)) &&
			(memBlock[index].data[1] == (qutioned & 0xFF00)>>8) &&
			(memBlock[index].data[2] == (qutioned & 0xFF0000)>>16) &&
			(memBlock[index].data[3] == (qutioned & 0xFF000000)>>24) &&
			((memBlock[index].data[4] & 0x01) == ((qutioned & 0x0100000000) >> 32))){
		//someone is sitting on the memBlock index & who is same as the khmer

		//now just update its neighbour
		//std::cout<<"now just update its neighbour"<<std::endl;
		switch(toupper(next_neucleotide)){//putting the neighbour
		case 'A'://1000
			memBlock[index].data[4] = memBlock[index].data[4] | 0x80;
			break;
		case 'T'://0100
			memBlock[index].data[4] = memBlock[index].data[4] | 0x40;
			break;
		case 'C'://0010
			memBlock[index].data[4] = memBlock[index].data[4] | 0x20;
			break;
		case 'G'://0001
			memBlock[index].data[4] = memBlock[index].data[4] | 0x10;
			break;
		}

		return MEM_BLOCK_FOUND; //khmer putted in the mem_block so return mem_block -- simply update the neighbour table
	}
	else if(whichHashFunc == LAST_HASH_FUNC){
		//only for the last time put in the alternative mem
		//neighbour field is not blank -- someone already sitting there
		//and who sitting there is not same as the given khmer
		//and the set pointer is null so it must be initialized
		//std::cout<<"else reached index "<<index<<" quotioned::"<<qutioned<< " data[0]=" <<(int)memBlock[index].data[0]<<" data[1]="<<(int)memBlock[index].data[1]<<std::endl;
		//expand the set pointer & place the khmer in there
		vector<struct vector_trial_5_byte>::iterator it;
		int index_vect=0;
		int index_v=index%num_of_vector;

		for(it=mapPointer5Byte[index_v].begin();it!=mapPointer5Byte[index_v].end();it++)
		{
			if(it->index_val==index)
			{
				//cout<<"existing value found for index "<<index<<" for khmer "<<khmer_h<<endl;
				break;
			}
			index_vect++;
		}


		if(it == mapPointer5Byte[index_v].end()){
			struct vector_trial_5_byte temp_vect;
			temp_vect.index_val=index;
			//temp_vect.vect_pnt=pp;
			mapPointer5Byte[index_v].push_back(temp_vect);
			//cout<<"new value inserted for index "<<index<<" for khmer "<<khmer_h<<" stored in index "<<index_v<<endl;
		}

		struct khmerInfo5Byte kin;//making a demo
		kin.data[0] = kin.data[1] = kin.data[2] = kin.data[3] = kin.data[4] = 0; //cleaning the demo info
		
		//cout<<"New value inserted in list for memblock index "<<index<<" for khmer "<<khmer_h<<" stored in list in index "<<index_v<<" hash func no "<<whichHashFunc<<endl;
		
		switch(toupper(next_neucleotide)){//putting the neighbour
		case 'A'://1000
			kin.data[4] = kin.data[4] | 0x80;
			break;
		case 'T'://0100
			kin.data[4] = kin.data[4] | 0x40;
			break;
		case 'C'://0010
			kin.data[4] = kin.data[4] | 0x20;
			break;
		case 'G'://s
			kin.data[4] = kin.data[4] | 0x10;
			break;
		}

		//putting the hash function data in the newly made khmerInfo5Byte obj
		
		/*****************************************************/
		//problem occur -- limon
		/*****************************************************/
		bool isFound = false;
		for(std::vector<khmerInfo5Byte>::iterator pi = mapPointer5Byte[index_v][index_vect].vect.begin();
				pi != mapPointer5Byte[index_v][index_vect].vect.end() ; pi++)
		{
			if(
				((pi->data[4] >> 1) & 0x07) == whichHashFunc && //hash value matched
				(pi->data[0] == (qutioned & 0xFF)) &&
				(pi->data[1] == (qutioned & 0xFF00)>>8) &&
				(pi->data[2] == (qutioned & 0xFF0000)>>16) &&
				(pi->data[3] == (qutioned & 0xFF000000)>>24) &&
				((pi->data[4] & 0x01) == ((qutioned & 0x0100000000) >> 32))){
			//(pi->data[0] == (qutioned & 0xFF)) && ((pi->data[1] & 0x0F) == ((qutioned & 0x0F00) >> 8))){
				pi->data[4] = pi->data[4] | kin.data[4];
				isFound = true;
				//cout<<"isFound true"<<endl;
				break;
			}
		}
		if(isFound == false){
		
			//putting the qutioned
			p = (char*) &qutioned;
			kin.data[0] = *p;
			p++;
			kin.data[1] = *p;
			p++;
			kin.data[2] = *p;
			p++;
			kin.data[3] = *p;
			p++;
			kin.data[4] = 0x01 & (*p);
			
			//putting "which hash value is stored"
			kin.data[4] = kin.data[4] | (whichHashFunc << 1); //--optimization		

			(mapPointer5Byte[index_v])[index_vect].vect.push_back(kin);
			
			return LIST;//new entry in the list
		}

		return LIST_FOUND; //khmer found in the list so return LIST_FOUND
	}

	//none of the if-else block is reached so it caused error
	//this means that particular khmer couldn't be placed inside the memblock
	//it will basically return for hash func 0,1,2 -> for hash func 3 if it fails
	//it will be stored in the list/vector/sqlite
	return ERROR;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x38, %rsp
movl %r9d, %ebx
movq %rcx, %rax
movq %rdx, %rcx
xorl %edx, %edx
divq %rcx
movq %rax, %r14
movl %r14d, %r12d
shrl $0x8, %r12d
movl %r14d, %edi
shrl $0x10, %edi
movl %r14d, %r9d
shrl $0x18, %r9d
movq %rax, %r15
shrq $0x20, %r15
leaq (%rdx,%rdx,4), %rax
leaq (%rsi,%rax), %rbp
movzbl 0x4(%rsi,%rax), %r13d
cmpl $0xf, %r13d
ja 0x3903
movb %r14b, (%rbp)
movb %r12b, 0x1(%rbp)
movb %dil, 0x2(%rbp)
movb %r9b, 0x3(%rbp)
andl $0x1, %r15d
leal (%r15,%rbx,2), %ebx
movb %bl, 0x4(%rbp)
movsbl %r8b, %edi
callq 0x3490
movl %eax, %ecx
movl $0x1, %eax
cmpl $0x46, %ecx
jg 0x3a43
cmpl $0x41, %ecx
je 0x3a6e
cmpl $0x43, %ecx
jne 0x3b41
orb $0x20, %bl
jmp 0x3a76
movq %rdx, %rcx
movl %r13d, %eax
shrl %eax
andl $0x7, %eax
cmpl %ebx, %eax
jne 0x396e
cmpb %r14b, (%rbp)
jne 0x396e
cmpb %r12b, 0x1(%rbp)
jne 0x396e
cmpb %dil, 0x2(%rbp)
jne 0x396e
cmpb %r9b, 0x3(%rbp)
jne 0x396e
movl %r13d, %eax
andl $0x1, %eax
movl %r15d, %edx
andl $0x1, %edx
cmpl %eax, %edx
jne 0x396e
movsbl %r8b, %edi
callq 0x3490
movl %eax, %ecx
movl $0x2, %eax
cmpl $0x46, %ecx
jg 0x3b60
cmpl $0x41, %ecx
je 0x3b70
cmpl $0x43, %ecx
jne 0x3b41
orb $0x20, %r13b
jmp 0x3b7a
leaq 0xa94f(%rip), %rdx      # 0xe2c4
movl $0x3, %eax
cmpl %ebx, (%rdx)
jne 0x3b41
leaq 0xa967(%rip), %rsi      # 0xe2f0
movq %rcx, %rax
xorl %edx, %edx
divq (%rsi)
xorl %ebp, %ebp
movq 0xa8f6(%rip), %rax      # 0xe290
movslq %edx, %rdx
leaq (%rdx,%rdx,2), %r13
movq (%rax,%r13,8), %rdx
movq 0x8(%rax,%r13,8), %rsi
cmpq %rsi, %rdx
je 0x39c1
cmpq %rcx, (%rdx)
je 0x39bf
incl %ebp
addq $0x20, %rdx
cmpq %rsi, %rdx
jne 0x39af
movl %ebp, %ebp
movl %r9d, 0xc(%rsp)
movl %edi, 0x8(%rsp)
cmpq %rsi, %rdx
jne 0x3a10
movl %r8d, 0x34(%rsp)
leaq (%rax,%r13,8), %rdi
xorps %xmm0, %xmm0
leaq 0x10(%rsp), %rsi
movups %xmm0, 0x8(%rsi)
movq $0x0, 0x18(%rsi)
movq %rcx, (%rsi)
callq 0x5e1e
movq 0x18(%rsp), %rdi
testq %rdi, %rdi
je 0x3a0b
movq 0x28(%rsp), %rsi
subq %rdi, %rsi
callq 0x3270
movl 0x34(%rsp), %r8d
movsbl %r8b, %edi
movb $0x0, 0x14(%rsp)
movl $0x0, 0x10(%rsp)
callq 0x3490
cmpl $0x46, %eax
jg 0x3a56
cmpl $0x41, %eax
movl 0x8(%rsp), %r8d
movl 0xc(%rsp), %r9d
je 0x3a7e
cmpl $0x43, %eax
jne 0x3a88
movb $0x20, %al
jmp 0x3a84
cmpl $0x47, %ecx
je 0x3a73
cmpl $0x54, %ecx
jne 0x3b41
orb $0x40, %bl
jmp 0x3a76
cmpl $0x47, %eax
movl 0x8(%rsp), %r8d
movl 0xc(%rsp), %r9d
je 0x3a82
cmpl $0x54, %eax
jne 0x3a88
movb $0x40, %al
jmp 0x3a84
orb $-0x80, %bl
jmp 0x3a76
orb $0x10, %bl
movb %bl, 0x4(%rbp)
jmp 0x3b41
movb $-0x80, %al
jmp 0x3a84
movb $0x10, %al
movb %al, 0x14(%rsp)
movq 0xa801(%rip), %rax      # 0xe290
movq (%rax,%r13,8), %rcx
shlq $0x5, %rbp
movq 0x8(%rcx,%rbp), %rax
movq 0x10(%rcx,%rbp), %rcx
cmpq %rcx, %rax
je 0x3ae4
movl %r15d, %edx
andl $0x1, %edx
movzbl 0x4(%rax), %esi
movl %esi, %edi
shrl %edi
andl $0x7, %edi
cmpl %ebx, %edi
jne 0x3adb
cmpb %r14b, (%rax)
jne 0x3adb
cmpb %r12b, 0x1(%rax)
jne 0x3adb
cmpb %r8b, 0x2(%rax)
jne 0x3adb
cmpb %r9b, 0x3(%rax)
jne 0x3adb
movl %esi, %edi
andl $0x1, %edi
cmpl %edi, %edx
je 0x3b50
addq $0x5, %rax
cmpq %rcx, %rax
jne 0x3aac
movb %r14b, 0x10(%rsp)
movb %r12b, 0x11(%rsp)
movb %r8b, 0x12(%rsp)
movb %r9b, 0x13(%rsp)
andl $0x1, %r15d
leal (%r15,%rbx,2), %eax
movb %al, 0x14(%rsp)
movq 0xa785(%rip), %rax      # 0xe290
movq (%rax,%r13,8), %rax
leaq (%rax,%rbp), %rdi
addq $0x8, %rdi
movq 0x8(%rdi), %rsi
cmpq 0x10(%rdi), %rsi
je 0x3b35
movb 0x14(%rsp), %al
movb %al, 0x4(%rsi)
movl 0x10(%rsp), %eax
movl %eax, (%rsi)
addq $0x5, 0x8(%rdi)
jmp 0x3b3f
leaq 0x10(%rsp), %rdx
callq 0x6244
xorl %eax, %eax
addq $0x38, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
orb 0x14(%rsp), %sil
movb %sil, 0x4(%rax)
movl $0x6, %eax
jmp 0x3b41
cmpl $0x47, %ecx
je 0x3b76
cmpl $0x54, %ecx
jne 0x3b41
orb $0x40, %r13b
jmp 0x3b7a
orb $-0x80, %r13b
jmp 0x3b7a
orb $0x10, %r13b
movb %r13b, 0x4(%rbp)
jmp 0x3b41
movq %rax, %rbx
movq 0x18(%rsp), %rdi
testq %rdi, %rdi
je 0x3b9a
movq 0x28(%rsp), %rsi
subq %rdi, %rsi
callq 0x3270
movq %rbx, %rdi
callq 0x3450