slepc-main 2024-11-09

VecCheckOrthonormality

Checks (or prints) the level of (bi-)orthonormality of a set of vectors.

#include "slepcvec.h" 
PetscErrorCode VecCheckOrthonormality(Vec V[],PetscInt nv,Vec W[],PetscInt nw,Mat B,PetscViewer viewer,PetscReal *lev)

Collective

	V	- a set of vectors
	nv	- number of V vectors
	W	- an alternative set of vectors (optional)
	nw	- number of W vectors
	B	- matrix defining the inner product (optional)
	viewer	- optional visualization context

lev

- level of orthogonality (optional)

This function is equivalent to VecCheckOrthonormality(), but in addition it checks that the diagonal of W'*V (or W'*B*V) is equal to all ones.

developer