! F01ELF Example Program Text ! Mark 24 Release. NAG Copyright 2012. Module f01elfe_mod ! F01ELF Example Program Module: ! Parameters and User-defined Routines ! .. Use Statements .. Use nag_library, Only: nag_wp ! .. Implicit None Statement .. Implicit None Contains Subroutine fcos2(iflag,nz,z,fz,iuser,ruser) ! .. Use Statements .. Use nag_library, Only: nag_wp ! .. Scalar Arguments .. Integer, Intent (Inout) :: iflag Integer, Intent (In) :: nz ! .. Array Arguments .. Complex (Kind=nag_wp), Intent (Out) :: fz(nz) Complex (Kind=nag_wp), Intent (In) :: z(nz) Real (Kind=nag_wp), Intent (Inout) :: ruser(*) Integer, Intent (Inout) :: iuser(*) ! .. Intrinsic Procedures .. Intrinsic :: cos ! .. Executable Statements .. Continue fz(1:nz) = cos((2.0E0_nag_wp,0.0E0_nag_wp)*z(1:nz)) Return End Subroutine fcos2 End Module f01elfe_mod Program f01elfe ! F01ELF Example Main Program ! .. Use Statements .. Use nag_library, Only: f01elf, nag_wp, x04caf Use f01elfe_mod, Only: fcos2 ! .. Implicit None Statement .. Implicit None ! .. Parameters .. Integer, Parameter :: nin = 5, nout = 6 ! .. Local Scalars .. Real (Kind=nag_wp) :: imnorm Integer :: i, ifail, iflag, lda, n ! .. Local Arrays .. Real (Kind=nag_wp), Allocatable :: a(:,:) Real (Kind=nag_wp) :: ruser(1) Integer :: iuser(1) ! .. Executable Statements .. Write (nout,*) 'F01ELF Example Program Results' Write (nout,*) ! Skip heading in data file Read (nin,*) Read (nin,*) n lda = n Allocate (a(lda,n)) ! Read A from data file Read (nin,*)(a(i,1:n),i=1,n) ! Find f( A ) ifail = 0 Call f01elf(n,a,lda,fcos2,iuser,ruser,iflag,imnorm,ifail) ! Print solution ifail = 0 Call x04caf('G','N',n,n,a,lda,'F(A) = COS(2A)',ifail) ! Print the norm of the imaginary part to check it is small Write (nout,*) Write (nout,Fmt='(1X,A,F6.2)') 'Imnorm =', imnorm End Program f01elfe