author | Don Pellegrino <don@drexel.edu> | 2010-01-25 22:44:52 (GMT) |
---|---|---|
committer | Don Pellegrino <don@drexel.edu> | 2010-01-25 22:44:52 (GMT) |
commit | 401233a57b8aa86e5b068710c0f77f7b8e205b3b (patch) (side-by-side diff) | |
tree | 0a830a608e8ef99d4704a0f19fefb0727a6fca2d | |
parent | 9184eadaef5518db7941d9999b4ce54dda088c73 (diff) | |
download | exp007-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.c | 10 |
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; |