Re: nullspaces
- To: mathgroup at smc.vnet.net
- Subject: [mg68988] Re: [mg68948] nullspaces
- From: Sseziwa Mukasa <mukasa at jeol.com>
- Date: Sat, 26 Aug 2006 02:04:24 -0400 (EDT)
- References: <200608250934.FAA09201@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
On Aug 25, 2006, at 5:34 AM, hespeler wrote:
> Hello
>
> I have a general problem with the computation of nullspaces for
> large dense matrices, e.g. dim(2,108), explicit example is added at
> the end of the mail. Using the internal command B=NullSpace[A] (A
> is an inexact, potentially complex matrix) and trying to reconfirm
> the result by
>
> Chop[A.ConjugateTranspose[NullSpace[B]]]==N[ZeroMatrix[Dimensions[A]
> [[1]],Dimensions[B][[1]]]]
The way NullSpace is defined the correct check is Chop[A.Transpose[B]]
==ZeroMatrix[First[Dimensions[A]],First[Dimensions[B]]] there is no
need to take conjugates or to apply NullSpace again.
for example:
In[122]:=
Block[{a=Table[Random[Complex,{-1-I,1+I}],{2},{108}],b},b=NullSpace[a];
Chop[a.Transpose[b]]==ZeroMatrix[First[Dimensions[a]],First
[Dimensions[b]]]]
Out[122]=
True
Regards,
Ssezi
- References:
- nullspaces
- From: hespeler <hespeler@gmx.de>
- nullspaces