ScaLAPACK Archives

[Scalapack] Issues regarding pdormrq vs. pdormr2

Hi again.
Forgot to say that descH=descT=
  descH(1) = 1
  descH(2) = 0
  descH(3) = 26
  descH(4) = 26
  descH(5) = 14
  descH(6) = 14
  descH(7) = 0
  descH(8) = 0
  descH(9) = 26

2011/9/12 Bj?rn Adlerborn <bjorn.adlerborn@Domain.Removed>

Hi!
After calling PDGERQF with the following parameters
(14, 14, T, 11, 11, descT, TAU, WORK(IPWORK), LDWORK1, IERR)

(i used LDWORK1 = -1 for getting the minimum amount of workspace, likewise
for the rest of the calls)

and i then call pdormrq with
('R', 'T', 14, 14, 14, T, 11, 11, descT, TAU, H, 11, 11, descH,
WORK(IPWORK), LDWORK1, IERR)

and compared with a call to pdormr2 with
('R', 'T', 14, 14, 14, T, 11, 11, descT, TAU, H, 11, 11, descH,
WORK(IPWORK), LDWORK1, IERR)

i get two diffrent results.
Seems like the first call (pdormrq) only applies the reflectors from 15:26,
and the later from 11:26 (which is correct).
Is it indended to be like this, i.e, that it only preforms blocked
applications of the reflections from first RQ call ?




--
Kindly regards Bj?rn

+46702048160





-- 
Mvh Bj?rn

0702048160
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
http://lists.eecs.utk.edu/mailman/private/scalapack/attachments/20110913/6e3c8a90/attachment.html
 

<Prev in Thread] Current Thread [Next in Thread>


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or