LMESetSolution#

Sets the placeholder for the solution of the matrix equation, as a low-rank matrix.

Synopsis#

#include "slepclme.h" 
PetscErrorCode LMESetSolution(LME lme,Mat X)

Collective

Input Parameters#

  • lme - the matrix function context

  • X - the solution matrix

Notes#

The matrix equation takes the general form AXE+DXB=C, where the solution matrix is of low rank and is written in factored form X = UDV’. This function provides a Mat object of type MATLRC that stores U, V and (optionally) D. These factors will be computed during LMESolve().

In equation types whose solution X is symmetric, such as Lyapunov, X must be created with V=U (or V=NULL).

If the user provides X with this function, then the solver will return a solution with rank at most the number of columns of U. Alternatively, it is possible to let the solver choose the rank of the solution, by setting X to NULL and then calling LMEGetSolution() after LMESolve().

See Also#

LMEGetSolution(), LMESetRHS(), LMESetProblemType(), LMESolve()

Level#

intermediate

Location#

src/lme/interface/lmesetup.c

Examples#

src/lme/tutorials/ex32.c


Index of all LME routines Table of Contents for all manual pages Index of all manual pages