petsc-3.14.6 2021-03-30

MatCreateLMVMSR1

Creates a limited-memory Symmetric-Rank-1 approximation matrix used for a Jacobian. L-SR1 is symmetric by construction, but is not guaranteed to be positive-definite.

Synopsis

#include "petscksp.h" 
PetscErrorCode MatCreateLMVMSR1(MPI_Comm comm, PetscInt n, PetscInt N, Mat *B)

The provided local and global sizes must match the solution and function vectors used with MatLMVMUpdate() and MatSolve(). The resulting L-SR1 matrix will have storage vectors allocated with VecCreateSeq() in serial and VecCreateMPI() in parallel. To use the L-SR1 matrix with other vector types, the matrix must be created using MatCreate() and MatSetType(), followed by MatLMVMAllocate(). This ensures that the internal storage and work vectors are duplicated from the correct type of vector.

Collective

Input Parameters

	comm	- MPI communicator, set to PETSC_COMM_SELF
	n	- number of local rows for storage vectors
	N	- global size of the storage vectors

Output Parameter

- the matrix

It is recommended that one use the MatCreate(), MatSetType() and/or MatSetFromOptions() paradigm instead of this routine directly.

Options Database Keys

-mat_lmvm_num_vecs

- maximum number of correction vectors (i.e.: updates) stored

Level

intermediate

Location

src/ksp/ksp/utils/lmvm/sr1/sr1.c
Index of all KSP routines
Table of Contents for all manual pages
Index of all manual pages