From 401233a57b8aa86e5b068710c0f77f7b8e205b3b Mon Sep 17 00:00:00 2001
From: Don Pellegrino <don@drexel.edu>
Date: Mon, 25 Jan 2010 22:44:52 +0000
Subject: Test of writing a sparse matrix to disk.

---
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;
--
cgit v0.8.3.1-22-g547a