summaryrefslogtreecommitdiffstats
Unidiff
-rw-r--r--src/aggregator.c34
-rw-r--r--src/assign/assign_protein_type.c141
-rw-r--r--src/load/load_influenza_aa_dat.c53
-rw-r--r--src/load/load_influenza_aa_dat.h2
-rw-r--r--src/load/load_influenza_faa.c53
-rw-r--r--src/load/load_influenza_faa.h2
-rw-r--r--src/updator.c2
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>
910
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"
1114
12int15int
13main ()16main ()
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 }
2137
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);
2644
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);
3151
32 /*52 /*
33 * Close the HDF5 file.53 * Close the HDF5 file.

Valid XHTML 1.0 Strict

Copyright © 2009 Don Pellegrino All Rights Reserved.