CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Any subrout to compute the orthogonal complement?

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   July 27, 2006, 02:35
Default Any subrout to compute the orthogonal complement?
Posts: n/a
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!
  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
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

All times are GMT -4. The time now is 14:13.