1!----------------------------------------------------------------------- 2 program test 3 4 include 'ceedf.h' 5 6 integer ceed,err 7 real*8 qr(12), tau(3) 8 9 character arg*32 10 11 qr = (/ 1, -1, 4, 1, 4, -2, 1, 4, 2, 1, -1, 0 /) 12 13 call getarg(1,arg) 14 15 call ceedinit(trim(arg)//char(0),ceed,err) 16 call ceedqrfactorization(ceed,qr,tau,4,3,err); 17 do i=1,12 18 if (abs(qr(i))<1.0D-14) then 19 qr(i) = 0 20 endif 21 write(*,'(A,F12.8)') '',qr(i) 22 enddo 23 do i=1,3 24 if (abs(tau(i))<1.0D-14) then 25 tau(i) = 0 26 endif 27 write(*,'(A,F12.8)') '',tau(i) 28 enddo 29 30 call ceeddestroy(ceed,err) 31 32 end 33!----------------------------------------------------------------------- 34