1 files changed, 10 insertions, 0 deletions

diff --git a/src/assign/assign_blast_scores.c b/src/assign/assign_blast_scores.c
index 1cd491a..b81caf8 100644
--- a/src/assign/assign_blast_scores.c
+++ b/src/assign/assign_blast_scores.c
@@ -21,6 +21,16 @@ assign_blast_scores (hid_t file_id)
   Mat M;
   MatCreateSeqAIJ (PETSC_COMM_SELF, 138769, 138769, 50, NULL, &M);
 
+  MatSetValue (M, 12, 8, 12.8, INSERT_VALUES);
+
+  /*
+   * Write the contents of the matrix to a HDF5 file.
+   */
+  PetscViewer viewer;
+  PetscErrorCode err = PetscViewerHDF5Open (PETSC_COMM_WORLD, "sparse.h5", FILE_MODE_WRITE, &viewer);
+  MatView (M, viewer);
+  PetscViewerDestroy (viewer);
+
   MatDestroy (M);
 
   return;