MFNSolve#

Solves the matrix function problem. Given a vector \(b\), the vector \(x = f(A)b\) is returned.

Synopsis#

#include "slepcmfn.h" 
PetscErrorCode MFNSolve(MFN mfn,Vec b,Vec x)

Collective

Input Parameters#

mfn - the matrix function solver context
b - the right hand side vector

Output Parameter#

x - the solution (this may be the same vector as b, then b will be overwritten with the answer)

Options Database Keys#

-mfn_view - print information about the solver used
-mfn_view_pre - print information about the solver before the solve starts
-mfn_view_mat - view the matrix
-mfn_view_rhs - view right hand side vector
-mfn_view_solution - view computed solution vector
-mfn_converged_reason - print reason for convergence/divergence, and number of iterations

Notes#

The matrix \(A\) is specified with MFNSetOperator(). The function \(f\) is specified via the FN object obtained with MFNGetFN() or set with MFNSetFN().

All the command-line options listed above admit an optional argument specifying the viewer type and options. For instance, use -mfn_view_mat binary:amatrix.bin to save the matrix to a binary file, or -mfn_view_solution :sol.m:ascii_matlab to save the solution in a file that can be executed in Matlab.