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! LCOV_EXCL_START 20 qr(i) = 0 21! LCOV_EXCL_STOP 22 endif 23 write(*,'(A,F12.8)') '',qr(i) 24 enddo 25 do i=1,3 26 if (abs(tau(i))<1.0D-14) then 27! LCOV_EXCL_START 28 tau(i) = 0 29! LCOV_EXCL_STOP 30 endif 31 write(*,'(A,F12.8)') '',tau(i) 32 enddo 33 34 call ceeddestroy(ceed,err) 35 36 end 37!----------------------------------------------------------------------- 38