|July 27, 2006, 02:35||
Any subrout to compute the orthogonal complement?
Any subrout in Lapack or IMSL or free code to compute the orthogonal complement of m vectors?
Now, there are m (m<n) vectors in a n-dimensional cartesian space at hand, which subroutine in Lapack can give the orthogonal complement of these m vectors. Orthogonal complement mean that the dot product of any vector U in the complement set and any vector v in the m-vector set at hand is zero.
I guess there are two ways to deal with this problem:
1)solve the equation Ax=0, while A is a mxn matrix whose row vectors are the m vectors at hand. But I don't know which subrout in Lapack can solve all the solutions to these underdetermined linear equations.
2)solve the equation AX=lamda*X,A is a nxn matrix whose first m rows are the m vectors the other row vectors are 0 vectors. And then find the eigenvector corresponding to the eigenvalue lamda=0. But I think this method takes a lot of time because it computes many unnecessary eigenvalues lamda/=0.
Any more efficient way? I prefer the first method!
|Thread||Thread Starter||Forum||Replies||Last Post|
|How to compute the gradient of a scalar as a post-processing||ayoros||OpenFOAM Post-Processing||15||August 7, 2015 04:48|
|[GAMBIT] how to generate orthogonal mesh using GAMBIT?||ivanbuz||ANSYS Meshing & Geometry||8||October 21, 2009 03:41|
|How do I compute the normal vector to a surface with FEM||Vasilis||Main CFD Forum||3||July 31, 2009 07:59|
|To compute Average flowfield with oodles||flying||OpenFOAM Running, Solving & CFD||3||May 5, 2009 08:46|
|Orthogonal curvilinear coordinates||Rudolf Pfeiffer||Main CFD Forum||6||December 10, 1998 11:20|