#include "slepcbv.h" PetscErrorCode BVOrthogonalize(BV V,Mat R)Collective on BV
V | - basis vectors |
V | - the modified basis vectors | |
R | - a sequential dense matrix (or NULL) |
If V has leading columns, then they are not modified (are assumed to be already orthonormal) and the corresponding part of R is not referenced.
Can pass NULL if R is not required.
The method to be used for block orthogonalization can be set with BVSetOrthogonalization(). If set to GS, the computation is done column by column with successive calls to BVOrthogonalizeColumn().
If V is rank-deficient or very ill-conditioned, that is, one or more columns are (almost) linearly dependent with respect to the rest, then the algorithm may break down or result in larger numerical error. Linearly dependent columns are essentially replaced by random directions, and the corresponding diagonal entry in R is set to (nearly) zero.
Location: src/sys/classes/bv/interface/bvorthog.c
Index of all BV routines
Table of Contents for all manual pages
Index of all manual pages