#include "slepcvec.h" PetscErrorCode VecCheckOrthogonality(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) |
If matrix B is provided then the check uses the B-inner product, W'*B*V.
If lev is not NULL, it will contain the maximum entry of matrix W'*V - I (in absolute value) omitting the diagonal. Otherwise, the matrix W'*V is printed.