| -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 @@ */ #include "assign/assign_protein_type.h" +#include "assign/assign_blast_scores.h" #include "error/check_h5_error.h" +#include <petsc.h> #include <stdio.h> #include <signal.h> #define FILE "influenza.h5" int -main () +main (int argc, char **argv) { /* + * Initialize the PETSc database and MPI. + * + * http://www.mcs.anl.gov/petsc/petsc-2/snapshots/petsc-dev/docs/manualpages/Sys/PetscInitialize.html#PetscInitialize + */ + PetscInitialize (&argc, &argv, 0, 0); + + /* * Open the HDF5 file. */ hid_t file_id = H5Fopen (FILE, H5F_ACC_RDWR, H5P_DEFAULT); @@ -22,7 +31,14 @@ main () /* * Assign protein type values to the sequence records. */ + /* assign_protein_type (file_id); + */ + + /* + * Assign pairwise BLAST scores. + */ + assign_blast_scores (file_id); /* * Close the HDF5 file. @@ -31,5 +47,7 @@ main () if (status < 0) check_h5_error (__FILE__, __LINE__); + PetscFinalize (); + return 0; } |
