summaryrefslogtreecommitdiffstats
authorDon Pellegrino <don@drexel.edu>2010-01-25 22:44:52 (GMT)
committer Don Pellegrino <don@drexel.edu>2010-01-25 22:44:52 (GMT)
commit401233a57b8aa86e5b068710c0f77f7b8e205b3b (patch) (side-by-side diff)
tree0a830a608e8ef99d4704a0f19fefb0727a6fca2d
parent9184eadaef5518db7941d9999b4ce54dda088c73 (diff)
downloadexp007-401233a57b8aa86e5b068710c0f77f7b8e205b3b.zip
exp007-401233a57b8aa86e5b068710c0f77f7b8e205b3b.tar.gz
exp007-401233a57b8aa86e5b068710c0f77f7b8e205b3b.tar.bz2
Test of writing a sparse matrix to disk.
-rw-r--r--src/assign/assign_blast_scores.c10
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;

Valid XHTML 1.0 Strict

Copyright © 2009 Don Pellegrino All Rights Reserved.