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