|
[Sponsors] |
Any subrout to compute the orthogonal complement? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 27, 2006, 02:35 |
Any subrout to compute the orthogonal complement?
|
#1 |
Guest
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! |
|
|
|
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 | 16 | March 21, 2018 06:02 |
[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 10:20 |