#include "slepcbv.h" PetscErrorCode BVOrthogonalizeColumn(BV bv,PetscInt j,PetscScalar *H,PetscReal *norm,PetscBool *lindep)Collective
bv | - the basis vectors context | |
j | - index of column to be orthogonalized |
H | - (optional) coefficients computed during orthogonalization | |
norm | - (optional) norm of the vector after being orthogonalized | |
lindep | - (optional) flag indicating that refinement did not improve the quality of orthogonalization |
Leading columns V[0..l-1] also participate in the orthogonalization, as well as the constraints. If H is given, it must have enough space to store j-l+1 coefficients (the last coefficient will contain the value norm, unless the norm argument is NULL).
If a non-standard inner product has been specified with BVSetMatrix(), then the vector is B-orthogonalized, using the non-standard inner product defined by matrix B. The output vector satisfies V[j]'*B*V[0..j-1] = 0.
This routine does not normalize the resulting vector, see BVOrthonormalizeColumn().
In the case of an indefinite inner product, the lindep parameter is not computed (set to false).