使用lapack进行qr分解遇到的问题

Transpose

R需要自己提取

[Fortran] 纯文本查看 复制代码

   complex(8),allocatable::r(:,:)
   allocate(r(nn,nn))
   r=0.d0
   do i=1,nn
      r(i,i:nn)=A(i,i:nn)
   end do

Q需要调用另一个子程序

[Fortran] 纯文本查看 复制代码

   call zungqr(mm,nn,nn,a,lda,tau,work,lwork,info)

这样就可以得到 A = Q * R