[Scalapack] Issues regarding pdormrq vs. pdormr2

[Björn Adlerborn]

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