<div dir="ltr"><div><div><div>Satish et al.,<br><br></div>I noticed a failure when testing with SeqAIJPERM and have made a small fix that is on branch rmills/fix-seqaijperm_seqmaij. If such a small fix needs to be tested in next, then please merge to next for testing. Otherwise, please merge to master.<br><br></div>Because it is so small, for convenience I have appended the git commit.<br><br></div>--Richard<br><br><span style="font-family:monospace,monospace">rmills@Richards-MBP:tests(rmills/fix-seqaijperm_seqmaij=)arch-osx-clang-gfortran7-debug$ git show<br>commit 1fde5765381347da5003c8b2045311701bd352b4 (HEAD -> rmills/fix-seqaijperm_seqmaij, origin/rmills/fix-seqaijperm_seqmaij)<br>Author: Richard Tran Mills <<a href="mailto:rmills@rmills.org">rmills@rmills.org</a>><br>Date:   Wed Aug 16 15:32:10 2017 -0700<br><br>    Associated function for MatPtAP_seqaijperm_seqmaij_C.<br><br>diff --git a/src/mat/impls/maij/maij.c b/src/mat/impls/maij/maij.c<br>index 83a4c181fe..b2c1711bd5 100644<br>--- a/src/mat/impls/maij/maij.c<br>+++ b/src/mat/impls/maij/maij.c<br>@@ -96,6 +96,7 @@ PetscErrorCode MatDestroy_SeqMAIJ(Mat A)<br>   ierr = PetscFree(A->data);CHKERRQ(ierr);<br>   ierr = PetscObjectComposeFunction((PetscObject)A,"MatConvert_seqmaij_seqaij_C",NULL);CHKERRQ(ierr);<br>   ierr = PetscObjectComposeFunction((PetscObject)A,"MatPtAP_seqaij_seqmaij_C",NULL);CHKERRQ(ierr);<br>+  ierr = PetscObjectComposeFunction((PetscObject)A,"MatPtAP_seqaijperm_seqmaij_C",NULL);CHKERRQ(ierr);<br>   PetscFunctionReturn(0);<br> }<br><br>@@ -3377,6 +3378,7 @@ PetscErrorCode  MatCreateMAIJ(Mat A,PetscInt dof,Mat *maij)<br>       }<br>       ierr = PetscObjectComposeFunction((PetscObject)B,"MatConvert_seqmaij_seqaij_C",MatConvert_SeqMAIJ_SeqAIJ);CHKERRQ(ierr);<br>       ierr = PetscObjectComposeFunction((PetscObject)B,"MatPtAP_seqaij_seqmaij_C",MatPtAP_SeqAIJ_SeqMAIJ);CHKERRQ(ierr);<br>+      ierr = PetscObjectComposeFunction((PetscObject)B,"MatPtAP_seqaijperm_seqmaij_C",MatPtAP_SeqAIJ_SeqMAIJ);CHKERRQ(ierr);<br>     } else {<br>       Mat_MPIAIJ  *mpiaij = (Mat_MPIAIJ*)A->data;<br>       Mat_MPIMAIJ *b;</span><br></div>