Pierre Jolivet <[email protected]> writes: >> On 29 May 2017, at 06:49, Jed Brown <[email protected]> wrote: >> >> You can't get a callback from this "external Krylov solver"? > > You mean, a callback to let the external solver use the memory allocated by > PETSc for the Mats? Unfortunately, I need the underlying memory to be > contiguous between all Mats, and I don't think it is possible with the > current PETSc operators for matrix creation.
You can allocate the full block, then pass portions of the (column-aligned) array in MatCreateDense().
signature.asc
Description: PGP signature
