1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
#include "parasail.h"
#include "parasail/matrices/blosum62.h"
#include "parasail/matrix_lookup.h"
int main(int argc, char **argv) {
const parasail_matrix_t *matrix = NULL;
parasail_matrix_t *user_matrix = NULL;
matrix = parasail_matrix_lookup("blosum62");
user_matrix = parasail_matrix_copy(matrix);
/* replace value at [2,4] location */
parasail_matrix_set_value(user_matrix, 2, 4, 200);
parasail_matrix_free(user_matrix);
user_matrix = parasail_matrix_create("ACGT", 2, -1);
/* replace value at [1,1] location, for a larger match */
parasail_matrix_set_value(user_matrix, 1, 1, 20);
parasail_matrix_free(user_matrix);
user_matrix = parasail_matrix_from_file("sample_matrix");
parasail_matrix_set_value(user_matrix, 3, 4, 50);
parasail_matrix_free(user_matrix);
return 0;
}
|
