Jed, It never worked in general. It was done originally and worked for BAIJ matrices with block size 4 to run PETSc-fun3d fast that is all it worked for.
I think it is much harder to get working completely then you think. It is more than just changing a few PetscScalar in the code to MatScalar, there are I think lots of situations where copies of stuff (with proper casting) have to be made to be passed in. Barry On Jun 5, 2009, at 8:40 AM, Jed Brown wrote: > > I don't see anything that would be particularly hard to fix, but it > would be a rather repetitive task. It looks like the only public > functions that need a change are MatCreateSeqAIJWithArrays and > MatCreateMPIAIJWithSplitArrays, but there is a lot of internal stuff. > Was there a particular reason MatScalar was deemed unimportant (lack > of > interoperability with any other packages comes to mind)? > > Jed >