debug显示写入位置冲突

huhelong

debug显示和运行显示见图片，不知道如何解决了

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

查看代码

打印?

01	program dggevexample

02	implicit none

03	integer,parameter ::n=3

04	real*8 :: a(n,n),b(n,n),alpha(n),beta(n)

05	real*8 :: vl(n,n),vr(n,n)

06	integer ::lda,ldb,ldvl,ldvr,lwork,info

07	real*8,allocatable :: work(:)

08	!Initialize input matrices a and b

09	a=reshape((/1.0,2.0,3.0,2.0,4.0,5.0,3.8,5.8,6.0/),[n,n])

10	b=reshape((/1.0,0.0,0.0,0.0,2.6,0.0,0.0,0.0,3.0/),[n,n])

11	!Set leading dimensions for input matrices

12

lda =n

13

ldb =n

14	!Allocate workspace

15	lwork = 4*n

16	allocate(work(lwork))

17	!Call DGGEV to compute eigenvalues and eigenvectors

18	call dggev('v','v',n,a,lda,b,ldb, alpha, beta,vl,ldvl,vr,ldvr,work,lwork,info)

19	!Print results

20

21	if(info.eq.0.d0)then

22

23	write(*,*)'Eigenvalues:'

24	write(*,'(3f10.6)')alpha

25	write(*,*)'Eigenvectors:'

26	write(*,'(3f10.6)') vl(:,1),vl(:,2),vl(:,3)

27

else

28	write(*,*)'Error in DGGEV:',info

29

end if

30	!Deallocate workspace

31	deallocate(work)

32	end program dggevexample

li913

dggev需要17个参数，你给了16个。

huhelong

谢谢您的指导，问题已解决