-rw-r--r-- | src/aggregator.c | 34 | ||||
-rw-r--r-- | src/assign/assign_protein_type.c | 141 | ||||
-rw-r--r-- | src/load/load_influenza_aa_dat.c | 53 | ||||
-rw-r--r-- | src/load/load_influenza_aa_dat.h | 2 | ||||
-rw-r--r-- | src/load/load_influenza_faa.c | 53 | ||||
-rw-r--r-- | src/load/load_influenza_faa.h | 2 | ||||
-rw-r--r-- | src/updator.c | 2 |
7 files changed, 201 insertions, 86 deletions
diff --git a/src/aggregator.c b/src/aggregator.c index c00d912..c9a03b5 100644 --- a/src/aggregator.c +++ b/src/aggregator.c | |||
@@ -6,28 +6,48 @@ | |||
6 | #include "error/check_h5_error.h" | 6 | #include "error/check_h5_error.h" |
7 | #include "load/load_influenza_aa_dat.h" | 7 | #include "load/load_influenza_aa_dat.h" |
8 | #include "load/load_influenza_faa.h" | 8 | #include "load/load_influenza_faa.h" |
9 | #include <stdio.h> | ||
9 | 10 | ||
10 | #define FILE "influenza.h5" | 11 | #define H5FILE "influenza.h5" |
12 | #define INFLUENZA_AA_DAT "/home/don/exp004/genomes/INFLUENZA/influenza_aa.dat" | ||
13 | #define INFLUENZA_FAA "/home/don/exp004/genomes/INFLUENZA/influenza.faa" | ||
11 | 14 | ||
12 | int | 15 | int |
13 | main () | 16 | main () |
14 | { | 17 | { |
15 | /* | 18 | /* |
16 | * Create the HDF5 file. | 19 | * Create a new HDF5 file if it does not already exist. If an |
20 | * existing file is found then open it. | ||
17 | */ | 21 | */ |
18 | hid_t file_id = H5Fcreate (FILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT); | 22 | hid_t file_id = 0; |
19 | if (file_id < 0) | 23 | FILE *f = fopen (H5FILE, "r+"); |
20 | check_h5_error (file_id, __FILE__, __LINE__); | 24 | if (f == NULL) |
25 | { | ||
26 | file_id = H5Fcreate (H5FILE, H5F_ACC_EXCL, H5P_DEFAULT, H5P_DEFAULT); | ||
27 | if (file_id < 0) | ||
28 | check_h5_error (file_id, __FILE__, __LINE__); | ||
29 | } | ||
30 | else | ||
31 | { | ||
32 | fclose (f); | ||
33 | file_id = H5Fopen (H5FILE, H5F_ACC_RDWR, H5P_DEFAULT); | ||
34 | if (file_id < 0) | ||
35 | check_h5_error (file_id, __FILE__, __LINE__); | ||
36 | } | ||
21 | 37 | ||
22 | /* | 38 | /* |
23 | * Load the supplementary protein data file. | 39 | * Load the supplementary protein data file. |
24 | */ | 40 | */ |
25 | load_influenza_aa_dat (file_id); | 41 | printf ("Loading \"influenza_aa.dat\" with contents of %s.\n", |
42 | INFLUENZA_AA_DAT); | ||
43 | load_influenza_aa_dat (file_id, INFLUENZA_AA_DAT); | ||
26 | 44 | ||
27 | /* | 45 | /* |
28 | * Load the FASTA protein sequence data file. | 46 | * Load the FASTA protein sequence data file. |
29 | */ | 47 | */ |
30 | load_influenza_faa (file_id); | 48 | printf ("Loading \"influenza.faa\" with contents of %s.\n", |
49 | INFLUENZA_FAA); | ||
50 | load_influenza_faa (file_id, INFLUENZA_FAA); | ||
31 | 51 | ||
32 | /* | 52 | /* |
33 | * Close the HDF5 file. | 53 | * Close the HDF5 file. |