When was the last time MatScalar was supported?

Barry Smith bsmith at mcs.anl.gov
Fri Jun 5 11:54:46 CDT 2009


     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.


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

More information about the petsc-dev mailing list