slepc-3.22.1 2024-10-28

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 A*X*E+D*X*B=C, where the solution matrix is of low rank and is written in factored form X = U*D*V'. 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().

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