Question

我有QR分解方法的问题。我使用dgeqrf子例程进行分解，但编译器中没有错误，但之后会出现问题。我没有发现错误在哪里。另一个问题是，A = Q * R =＆gt;如果矩阵为零，则分解为零或丢失秩。

program decomposition

!CONTAINS
!subroutine Qrdecomposition(A_mat, R)
real,dimension(2,2)   :: A_mat    !real,dimension(2,2),intent(inout)   
:: A_mat
real,dimension(2,2)   :: R        !real,dimension(2,2),intent(out)     
:: R
real,dimension(2,2)                  :: A
integer                              :: M,N,LDA,LWORK,INFO
real,allocatable, dimension(:,:)     :: TAU
real,allocatable, dimension(:,:)     :: WORK
external   dgeqrf
M=2
N=2
LDA=2
LWORK=2
INFO=0
A_mat(1,1)=4
A_mat(1,2)=1
A_mat(2,1)=3
A_mat(2,2)=1
A=A_mat

call dgeqrf(M,N,A,TAU,WORK,LWORK,INFO)
R=A
print *,R,WORK,LWORK

!end subroutine Qrdecomposition
end program decomposition

Answer 1

我看到你的代码中有三个错误：

1）您忘记了LDA的{{1}}参数，

2）必须明确分配dgeqrf和TAU，

3）所有数组都应该以双精度声明，以与WORK接口一致：

dgeqrf

在某些情况下，Fortran会执行数组的自动分配，但通常不应该依赖它，并且不是这种情况。

编辑 roygvib指出了第3点，见下文。

QR分解Fortran错误

1 个答案: