[petsc-dev] I think TS should have a TSSetMassMatrix()

Barry Smith bsmith at mcs.anl.gov
Sun Jul 23 13:42:13 CDT 2017

   Then explicit methods could automatically use a solve to apply the mass matrix inverse for each "right hand side" application and implicit methods could automatically "drop in" the mass matrix with the shift when a right hand side Jacobian is given also. Instead of requiring the user to handle this themselves. Are there any gotcha's that I am missing or is this relatively straightforward?

   Why wouldn't we have this interface? 

   Note that adjoints are slightly tricky with a mass matrix without this API.

    Almost for sure MFEM must have this type of API.

More information about the petsc-dev mailing list