possible typos in descriptions

Post here if you want to report a bug to the LAPACK team

possible typos in descriptions

Postby Yana Gurieva » Wed Jun 08, 2011 2:24 am

Functions c/zlarfg and c/zlarfgp have now the following formula in Purpose description
H**H * ( alpha ) = ( beta ), H**H * H = I.
( x ) ( 0 )

Should not it be
H * ( alpha ) = ( beta ), H**H * H = I.
( x ) ( 0 )

instead?
Yana Gurieva
 
Posts: 7
Joined: Mon Apr 11, 2011 1:26 am

Re: possible typos in descriptions

Postby Julien Langou » Wed Jun 08, 2011 3:04 am

Hi Yana,

There is a difference between the real case and the complex case. In the real case, the elementary Householder reflections ("H") are symmetric. ( H = H**T.) In the complex case, the elementary Householder reflections ("H") are not Hermitian. ( H is not H**H.)

In LAPACK clarfg, H is defined as the Householder transformation such that

Code: Select all
[            ]**H * [ alpha ] = [ beta ]
[            ]      [   x1  ]   [  0   ]
[     H1     ]      [   x2  ]   [  0   ]
[            ]      [   x3  ]   [  0   ]
[            ]      [   x4  ]   [  0   ]


The **H is necessary.

Cheers,
Julien.
Julien Langou
 
Posts: 734
Joined: Thu Dec 09, 2004 12:32 pm
Location: Denver, CO, USA


Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest