[petsc-users] Using slepc MFN

[Guglielmo, Tyler Hardy]

Not sure if there are many slepc users here, but I have a question regarding MFN in slepc.

Essentially I am running a loop over time steps, and at each step I would like to compute a matrix exponential, however the matrix that I am computing the exponential of changes in a non-trivial way at each time step.

It seems like the proper way to implement this is to build and destroy an MFN object at each time step.  However, this doesn’t seem terribly efficient as most Krylov methods I have looked at (Expokit) require a large amount of workspace memory allocation at runtime.

Is there a safe way to change the operator appearing in the MFN object without having to destroy and recreate?

Best,
Tyler

+++++++++++++++++++++++++++++
Tyler Guglielmo
Postdoctoral Researcher
Lawrence Livermore National Lab
Office: 925-423-6186
Cell: 210-480-8000
+++++++++++++++++++++++++++++
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221215/370e0e60/attachment-0001.html>