obidmscolumn_seq.c 3 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
/****************************************************************************
 * OBIDMS_column_seq functions                                              *
 ****************************************************************************/

/**
 * @file obidsmcolumn_seq.c
 * @author Celine Mercier
 * @date November 18th 2015
 * @brief Functions handling OBIColumns containing data in the form of indices referring to DNA sequences.
 */


#include <stdlib.h>
#include <stdio.h>

#include "obidmscolumn.h"
#include "obitypes.h"
Celine Mercier's avatar
Celine Mercier committed
18
#include "dna_seq_indexer.h"
19 20 21 22 23 24 25 26


/**********************************************************************
 *
 * D E F I N I T I O N   O F   T H E   P U B L I C   F U N C T I O N S
 *
 **********************************************************************/

27
int obi_column_set_obiseq_with_elt_idx(OBIDMS_column_p column, index_t line_nb, index_t element_idx, const char* value)
28
{
29
	index_t idx;
30
	char*   new_indexer_name;
31

32 33
	if (obi_column_prepare_to_set_value(column, line_nb) < 0)
		return -1;
34

35 36
	// TODO NA

37
	// Add the value in the indexer
Celine Mercier's avatar
Celine Mercier committed
38
	idx = obi_index_dna_seq(column->indexer, value);
39 40 41 42
	if (idx == -1)	// An error occurred
	{
		if (obi_errno == OBI_READ_ONLY_INDEXER_ERROR)
		{
43
			// TODO PUT IN A COLUMN FUNCTION
44 45 46 47 48
			// If the error is that the indexer is read-only, clone it
			new_indexer_name = obi_build_indexer_name((column->header)->name, (column->header)->version);
			if (new_indexer_name == NULL)
				return -1;
			column->indexer = obi_clone_indexer(column->indexer, new_indexer_name);	// TODO Need to lock this somehow?
49 50 51 52
			strcpy((column->header)->indexer_name, new_indexer_name);
			free(new_indexer_name);
			obi_set_errno(0);

53 54 55 56 57 58 59 60
			// Add the value in the new indexer
			idx = obi_index_dna_seq(column->indexer, value);
			if (idx == -1)
				return -1;
		}
		else
			return -1;
	}
61
	// Add the value's index in the column
Celine Mercier's avatar
Celine Mercier committed
62
	*(((index_t*) (column->data)) + (line_nb * ((column->header)->nb_elements_per_line)) + element_idx) = idx;
63 64 65 66 67

	return 0;
}


68
char* obi_column_get_obiseq_with_elt_idx(OBIDMS_column_p column, index_t line_nb, index_t element_idx)
69
{
Celine Mercier's avatar
Celine Mercier committed
70
	index_t idx;
71

72
	if (obi_column_prepare_to_get_value(column, line_nb) < 0)
Celine Mercier's avatar
Celine Mercier committed
73
		return OBISeq_NA;
74

Celine Mercier's avatar
Celine Mercier committed
75
	idx = *(((index_t*) (column->data)) + (line_nb * ((column->header)->nb_elements_per_line)) + element_idx);
76 77 78

	// Check NA
	if (idx == OBIIdx_NA)
Celine Mercier's avatar
Celine Mercier committed
79
		return OBISeq_NA;
80

Celine Mercier's avatar
Celine Mercier committed
81
	return obi_retrieve_dna_seq(column->indexer, idx);
82 83 84
}


85
int obi_column_set_obiseq_with_elt_name(OBIDMS_column_p column, index_t line_nb, const char* element_name, const char* value)
86
{
87
	index_t element_idx = obi_column_get_element_index_from_name(column, element_name);
88 89
	if (element_idx == OBIIdx_NA)
		return -1;
Celine Mercier's avatar
Celine Mercier committed
90

91
	return obi_column_set_obiseq_with_elt_idx(column, line_nb, element_idx, value);
92 93 94
}


95
char* obi_column_get_obiseq_with_elt_name(OBIDMS_column_p column, index_t line_nb, const char* element_name)
96
{
97
	index_t element_idx = obi_column_get_element_index_from_name(column, element_name);
98
	if (element_idx == OBIIdx_NA)
Celine Mercier's avatar
Celine Mercier committed
99 100
		return OBISeq_NA;

101
	return obi_column_get_obiseq_with_elt_idx(column, line_nb, element_idx);
Celine Mercier's avatar
Celine Mercier committed
102 103
}