IPOrthogonalize

Orthogonalize a vector with respect to a set of vectors.

Synopsis

#include "slepcip.h" 
PetscErrorCode IPOrthogonalize(IP ip,PetscInt n,PetscTruth *which,Vec *V,Vec v,PetscScalar *H,PetscReal *norm,PetscTruth *lindep,Vec work,PetscScalar* swork)
Collective on IP

Input Parameters

ip - the inner product (IP) context
n - number of columns of V
which - logical array indicating columns of V to be used
V - set of vectors
work - workspace vector
swork - workspace array

Input/Output Parameter

v - (input) vector to be orthogonalized and (output) result of orthogonalization

Output Parameter

H - coefficients computed during orthogonalization
norm - norm of the vector after being orthogonalized
lindep - flag indicating that refinement did not improve the quality of orthogonalization

Notes

This function applies an orthogonal projector to project vector v onto the orthogonal complement of the span of the columns of V.

On exit, v0 = [V v]*H, where v0 is the original vector v.

This routine does not normalize the resulting vector.

See Also

IPSetOrthogonalization(), IPBiOrthogonalize()

Location: src/ip/iporthog.c
Index of all IP routines
Table of Contents for all manual pages
Index of all manual pages