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;
 }