In addition, NAG recommends that before calling any library routine you should read the following reference material (see Section 5):
(a) Essential Introduction
(b) Chapter Introduction
(c) Routine Document
If you intend to use the NAG library within a multithreaded application please refer to the document on Thread Safety.
The libraries supplied with this implementation have been compiled in a manner that facilitates the use of multiple threads. Information on how to use the libraries in a multithreaded application is given in Section 3.1.
Assuming that libnag_noblas.a and the NEC tuned BLAS have been installed in a directory in the search path of the linker, such as /usr/lib, then you can link to the NAG Fortran Library, utilising the NEC tuned BLAS, in the following manner:
f90 -float0 -P stack -Wf,-w,double16 driver.f -lnag_noblas -lblas_64where driver.f is your application program.
If the NEC tuned BLAS are not available, libnag.a should have been installed instead of libnag_noblas.a; if this is in a directory in the search path of the linker, such as /usr/lib, then you can link to the NAG Fortran Library in the following manner:
f90 -float0 -P stack -Wf,-w,double16 driver.f -lnag
If you are developing a multi-threaded application, you should include
the
f90 -float0 -P stack -Wf,-pvctl,vwork=stack driver.f -lnag_noblas -lblas_64or (if the NEC tuned BLAS are not available):
f90 -float0 -P stack -Wf,-pvctl,vwork=stack driver.f -lnagNEC tuned BLAS are part of a product called Mathkeisan/SX and can be obtained by making a request to technical@atcc.necsys.com.
nagexample c06eaewill copy the example program and its data into the files c06eaee.f and c06eaee.d in the current directory and process them to produce the example program results.
In the NAG Fortran Library Manual, routine documents that have been typeset since Mark 12 present the example programs in a generalised form, using bold italicised terms as described in Section 3.3.
The example programs supplied to a site in machine-readable form have been modified as necessary so that they are suitable for immediate execution. Note that the distributed example programs are those used in this implementation and may not correspond exactly with the programs published in the manual. The distributed example programs should be used in preference wherever possible.
The distributed example results are those obtained using the library libnag_noblas.a and the NEC tuned BLAS; using libnag.a and the NAG BLAS may give slightly different results.
real - REAL (REAL*8)
basic precision - single precision
complex - COMPLEX (COMPLEX*16)
additional precision - double precision (REAL*16)
machine precision - the machine precision, see the value
returned by X02AJE in Section 4
Thus a parameter described as real should be declared as
REAL in your program. If a routine accumulates an inner product in
additional
precision, it is using double precision.
In routine documents that have been typeset since Mark 12 additional bold italicised terms are used in the published example programs and they must be interpreted as follows:
real as an intrinsic function name - REAL imag - AIMAG cmplx - CMPLX conjg - CONJG e in constants, e.g. 1.0e-4 - E, e.g. 1.0E-4 e in formats, e.g. e12.4 - E, e.g. E12.4
All references to routines in Chapter F07 - Linear Equations (LAPACK)
and Chapter F08 - Least-squares and Eigenvalue Problems (LAPACK) use the
LAPACK name, not the NAG F07/F08 name. The LAPACK name is precision dependent,
and hence the name appears in a bold italicised typeface.
The typeset examples use the single precision form of the LAPACK name.
To convert this name to its double precision form, change the first character
either from S to D or C to Z as appropriate.
For example:
sgetrf refers to the LAPACK routine name - SGETRF cpotrs - CPOTRS
A02AAF denotes the double precision version A02AAE denotes the single precision version
The names of auxiliary routines have also been modified by interchanging the first three and the last three characters, e.g. C02AFZ has been changed to AFZC02.
In the NAG Fortran Library Manual almost all library routine names end in F. Therefore, when using the manual in conjunction with this single precision implementation, all such names must be understood to refer to the single precision versions with names ending in E. Some routines in the Library require users to specify particular auxiliary routines. Again, when using this implementation it is necessary to specify the amended names.
Chapters C05, D03 and E04 have some routines whose last letter is A rather than F (B in single precision implementations). An A version is always paired with an F routine, the A version being safe to use in a multithreaded environment, but otherwise having identical functionality to the F version. When using the manual in conjunction with this single precision implementation, all names ending in A must be understood to refer to the single precision versions with names ending in B.
The names of COMMON blocks have also been modified, e.g. AD02BJ is renamed BJD02A. This is unlikely to affect the user.
(a) Subroutines are called as such
(b) Functions are declared with the right type
(c) The correct number of arguments are passed
(d) All arguments match in type and structure
These interface blocks have been generated automatically by analysing the source code for the NAG Fortran Library. As a consequence, and because these files have been thoroughly tested, their use is recommended in preference to writing your own declarations.
The NAG Fortran Library Interface Block files are organised by Library chapter. The module names are:
nag_f77_a_chapter nag_f77_c_chapter nag_f77_d_chapter nag_f77_e_chapter nag_f77_f_chapter nag_f77_g_chapter nag_f77_h_chapter nag_f77_m_chapter nag_f77_p_chapter nag_f77_s_chapter nag_f77_x_chapterThese are supplied in pre-compiled form (.mod files) and they can be accessed by specifying the
In order to make use of these modules from existing Fortran 77 code the following changes need to be made:
* S18DEE Example Program Text
* Mark 14 Revised. NAG Copyright 1989.
*******************************************************************
* Add USE statement for relevant chapters *
USE NAG_F77_S_CHAPTER
* *
*******************************************************************
* .. Parameters ..
INTEGER NIN, NOUT
PARAMETER (NIN=5,NOUT=6)
INTEGER N
PARAMETER (N=2)
* .. Local Scalars ..
COMPLEX Z
REAL FNU
INTEGER IFAIL, NZ
CHARACTER*1 SCALE
* .. Local Arrays ..
COMPLEX CY(N)
* .. External Subroutines ..
*******************************************************************
* EXTERNAL declarations need to be removed (and type declarations *
* for functions). *
C EXTERNAL S18DEE
* *
*******************************************************************
* .. Executable Statements ..
WRITE (NOUT,*) 'S18DEE Example Program Results'
* Skip heading in data file
READ (NIN,*)
WRITE (NOUT,*)
WRITE (NOUT,99999) 'Calling with N =', N
WRITE (NOUT,*)
WRITE (NOUT,*)
+' FNU Z SCALE CY(1) CY(2)
+ NZ IFAIL'
WRITE (NOUT,*)
20 READ (NIN,*,END=40) FNU, Z, SCALE
IFAIL = 0
*
CALL S18DEE(FNU,Z,N,SCALE,CY,NZ,IFAIL)
*
WRITE (NOUT,99998) FNU, Z, SCALE, CY(1), CY(2), NZ, IFAIL
GO TO 20
40 STOP
*
99999 FORMAT (1X,A,I2)
99998 FORMAT (1X,F7.4,' (',F7.3,',',F7.3,') ',A,
+ 2(' (',F7.3,',',F7.3,')'),I4,I4)
END
S07AAE F(1) = 1.0E+13
F(2) = 1.0E-14
S10AAE E(1) = 18.50
S10ABE E(1) = 708.0
S10ACE E(1) = 708.0
S13AAE x(hi) = 708.3
S13ACE x(hi) = 2.2E+15
S13ADE x(hi) = 2.2E+15
S14AAE IFAIL = 1 if X > 170.
IFAIL = 2 if X < -170.
IFAIL = 3 if abs(X) < 2.23E-308
S14ABE IFAIL = 2 if X > 2.55E+305
S15ADE x(hi) = 26.6
x(low) = -6.25
S15AEE x(hi) = 6.25
S17ACE IFAIL = 1 if X > 2.2E+15
S17ADE IFAIL = 1 if X > 2.2E+15
IFAIL = 3 if 0.0 < X <= 2.23E-308
S17AEE IFAIL = 1 if abs(X) > 2.2E+15
S17AFE IFAIL = 1 if abs(X) > 2.2E+15
S17AGE IFAIL = 1 if X > 103.8
IFAIL = 2 if X < -2.2E+10
S17AHE IFAIL = 1 if X > 104.1
IFAIL = 2 if X < -2.2E+10
S17AJE IFAIL = 1 if X > 104.1
IFAIL = 2 if X < -1.8E+9
S17AKE IFAIL = 1 if X > 104.1
IFAIL = 2 if X < -1.8E+9
S17DCE IFAIL = 2 if abs (Z) < 5.97E-154
IFAIL = 4 if abs (Z) or FNU+N-1 > 6.71E+7
IFAIL = 5 if abs (Z) or FNU+N-1 > 4.50E+15
S17DEE IFAIL = 2 if imag (Z) > 700.
IFAIL = 3 if abs (Z) or FNU+N-1 > 6.71E+7
IFAIL = 4 if abs (Z) or FNU+N-1 > 4.50E+15
S17DGE IFAIL = 3 if abs (Z) > 1.65E+5
IFAIL = 4 if abs (Z) > 2.72E+10
S17DHE IFAIL = 3 if abs (Z) > 1.65E+5
IFAIL = 4 if abs (Z) > 2.72E+10
S17DLE IFAIL = 2 if abs (Z) < 5.97E-154
IFAIL = 4 if abs (Z) or FNU+N-1 > 6.71E+7
IFAIL = 5 if abs (Z) or FNU+N-1 > 4.50E+15
S18ADE IFAIL = 2 if 0.0 < X <= 2.23E-308
S18AEE IFAIL = 1 if abs(X) > 711.6
S18AFE IFAIL = 1 if abs(X) > 711.6
S18CDE IFAIL = 2 if 0.0 < X <= 2.23E-308
S18DCE IFAIL = 2 if abs (Z) < 5.97E-154
IFAIL = 4 if abs (Z) or FNU+N-1 > 6.71E+7
IFAIL = 5 if abs (Z) or FNU+N-1 > 4.50E+15
S18DEE IFAIL = 2 if real (Z) > 700.
IFAIL = 3 if abs (Z) or FNU+N-1 > 6.71E+7
IFAIL = 4 if abs (Z) or FNU+N-1 > 4.50E+15
S19AAE IFAIL = 1 if abs(x) >= 49.50
S19ABE IFAIL = 1 if abs(x) >= 49.50
S19ACE IFAIL = 1 if X > 997.26
S19ADE IFAIL = 1 if X > 997.26
S21BCE IFAIL = 3 if an argument < 1.579E-205
IFAIL = 4 if an argument >= 3.774E+202
S21BDE IFAIL = 3 if an argument < 2.820E-103
IFAIL = 4 if an argument >= 1.404E+102
X01AAE (PI) = 3.1415926535897932 X01ABE (GAMMA) = 0.5772156649015329
The basic parameters of the model
X02BHE = 2 X02BJE = 53 X02BKE = -1021 X02BLE = 1024 X02DJE = .TRUE.Derived parameters of the floating-point arithmetic
X02AJE = Z'3CA0000000000001' ( 1.11022302462516E-16) X02AKE = Z' 10000000000000' ( 2.22507385850720E-308) X02ALE = Z'7FEFFFFFFFFFFFFF' ( 1.79769313486232E+308) X02AME = Z' 10000000000000' ( 2.22507385850720E-308) X02ANE = Z'2010000000000000' ( 2.98333629248008E-154)Parameters of other aspects of the computing environment
X02AHE = Z'4320000000000000' ( 2.25179981368525E+15 ) X02BBE = 9223372036854775807 X02BEE = 15 X02DAE = .FALSE.
A full online version of the NAG Fortran Library Manual is supplied in the form of Portable Document Format (PDF) files, with an HTML index, in the NAGdoc directory. The introductory material is also provided as HTML files in the NAGdoc directory.
A main index file has been provided (NAGdoc/fl/html/mark20.html) which contains a fully linked contents document pointing to all the available PDF (and where available HTML) files. Use your HTML browser to navigate from here.
In addition the following are provided in the doc directory:
The NAG Response Centres are available for general enquiries from all users and also for technical queries from sites with an annually licensed product or support service.
The Response Centres are open during office hours, but contact is possible by fax, email and phone (answering machine) at all times.
When contacting a Response Centre it helps us deal with your enquiry quickly if you can quote your NAG site reference and NAG product code (in this case FLNE520SA).
The NAG websites provide information about implementation availability, descriptions of products, downloadable software, product documentation and technical reports. The NAG websites can be accessed at
http://www.nag.co.uk/, http://www.nag.com/ (in North America) or http://www.nag-j.co.jp/ (in Japan)
NAG Ltd Wilkinson House Jordan Hill Road OXFORD OX2 8DR NAG Ltd Response Centre United Kingdom email: support@nag.co.uk Tel: +44 (0)1865 511245 Tel: +44 (0)1865 311744 Fax: +44 (0)1865 310139 Fax: +44 (0)1865 310139 NAG Inc 1431 Opus Place, Suite 220 Downers Grove IL 60515-1362 NAG Inc Response Center USA email: infodesk@nag.com Tel: +1 630 971 2337 Tel: +1 630 971 2345 Fax: +1 630 971 2706 Fax: +1 630 971 2706 Nihon NAG KK Hatchobori Frontier Building 2F 4-9-9 Hatchobori Chuo-ku Tokyo 104-0032 Japan email: help@nag-j.co.jp Tel: +81 (0)3 5542 6311 Fax: +81 (0)3 5542 6312