-rw-r--r-- | configure.ac | 62 | ||||
-rw-r--r-- | doc/build-cobalt.txt | 24 | ||||
-rw-r--r-- | src/Makefile.am | 6 | ||||
-rw-r--r-- | src/assign/assign_blast_scores.c | 27 | ||||
-rw-r--r-- | src/assign/assign_blast_scores.h | 12 | ||||
-rw-r--r-- | src/updator.c | 20 |
6 files changed, 127 insertions, 24 deletions
diff --git a/src/updator.c b/src/updator.c index f93e205..0bade65 100644 --- a/src/updator.c +++ b/src/updator.c | |||
@@ -3,16 +3,25 @@ | |||
3 | */ | 3 | */ |
4 | 4 | ||
5 | #include "assign/assign_protein_type.h" | 5 | #include "assign/assign_protein_type.h" |
6 | #include "assign/assign_blast_scores.h" | ||
6 | #include "error/check_h5_error.h" | 7 | #include "error/check_h5_error.h" |
8 | #include <petsc.h> | ||
7 | #include <stdio.h> | 9 | #include <stdio.h> |
8 | #include <signal.h> | 10 | #include <signal.h> |
9 | 11 | ||
10 | #define FILE "influenza.h5" | 12 | #define FILE "influenza.h5" |
11 | 13 | ||
12 | int | 14 | int |
13 | main () | 15 | main (int argc, char **argv) |
14 | { | 16 | { |
15 | /* | 17 | /* |
18 | * Initialize the PETSc database and MPI. | ||
19 | * | ||
20 | * http://www.mcs.anl.gov/petsc/petsc-2/snapshots/petsc-dev/docs/manualpages/Sys/PetscInitialize.html#PetscInitialize | ||
21 | */ | ||
22 | PetscInitialize (&argc, &argv, 0, 0); | ||
23 | |||
24 | /* | ||
16 | * Open the HDF5 file. | 25 | * Open the HDF5 file. |
17 | */ | 26 | */ |
18 | hid_t file_id = H5Fopen (FILE, H5F_ACC_RDWR, H5P_DEFAULT); | 27 | hid_t file_id = H5Fopen (FILE, H5F_ACC_RDWR, H5P_DEFAULT); |
@@ -22,7 +31,14 @@ main () | |||
22 | /* | 31 | /* |
23 | * Assign protein type values to the sequence records. | 32 | * Assign protein type values to the sequence records. |
24 | */ | 33 | */ |
34 | /* | ||
25 | assign_protein_type (file_id); | 35 | assign_protein_type (file_id); |
36 | */ | ||
37 | |||
38 | /* | ||
39 | * Assign pairwise BLAST scores. | ||
40 | */ | ||
41 | assign_blast_scores (file_id); | ||
26 | 42 | ||
27 | /* | 43 | /* |
28 | * Close the HDF5 file. | 44 | * Close the HDF5 file. |
@@ -31,5 +47,7 @@ main () | |||
31 | if (status < 0) | 47 | if (status < 0) |
32 | check_h5_error (__FILE__, __LINE__); | 48 | check_h5_error (__FILE__, __LINE__); |
33 | 49 | ||
50 | PetscFinalize (); | ||
51 | |||
34 | return 0; | 52 | return 0; |
35 | } | 53 | } |