In addition, NAG recommends that before calling any Library routine you should read the following reference material (see Section 6):
(a) Essential Introduction
(b) Chapter Introduction
(c) Routine Document
The libraries supplied with this implementation have not been compiled in a manner that facilitates the use of multiple threads.
Fortran 90/95 users are advised that the compiled *.mod files (the interface blocks) have been compiled with the Intel Fortran Compiler 8.1 and are intended for use with that compiler. Users may have to compile the interface blocks themselves if they wish to use them with a different compiler.
When the DLLs are used with a non-Intel compiler, please note that two input/output systems are in use: those of Intel for library routines and of course the compiler's own input/output routines for the calling program. This means that programs like the E04UDF example program cannot read the data from just one file. This is because the program reads some of the data using its input/output system. When the option setting routine tries to read the data file, the Intel input/output routines are used. The two input/output systems are completely disjoint and so in particular Intel has no knowledge of the position in the data file that the program input/output system has reached. The problem is circumvented by having two separate data files. Routines affected by this are mainly the option setting routines in chapters H02 and E04.
If using the MKL-based NAG Fortran Library DLL FLDLL214Z_mkl.dll, it is important to ensure that the compiler's floating-point exception handling is set to allow full IEEE 754 floating-point exception behaviour (non-stop arithmetic), otherwise, in certain circumstances, an arithmetic exception may be thrown and the program may terminate abnormally. For the Intel compiler the option is /fpe:3, which is the default behaviour. For the NAGWare f95 compiler the option is -ieee=full, which is not the default behaviour for f95. Please check your compiler documentation for the correct option settings. Note that the operation of the self-contained NAG DLL FLDLL214Z_nag.dll is not affected by the floating-point exception settings.
c:\Program Files\NAG\FL21\fldll214zlIf this folder does not exist, please consult the system manager (or the person who did the installation). In some of the following subsections, this folder is referred to as install dir.
We also assume that the default shortcut for the Library command prompt is placed in the Start Menu under:
Start|All Programs|NAG|FL21|
If this shortcut does not exist, please consult the system manager (or the person who did the installation). (Other shortcuts created as part of the Library installation procedure are also assumed to be in this location.)
To ensure that the NAG DLLs are accessible at runtime, the install dir\bin folder must be on the path. If the MKL version of the DLL is to be used, the install dir\mkl_ia32_8.0\bin folder must also be on the path, but should appear later in the path than the install dir\bin folder, since the NAG versions of a few BLAS/LAPACK routines are included in FLDLL214Z_mkl.dll to avoid problems with the vendor versions. (See Section 5 for details.)
To check the accessibility of the NAG DLLs, run the program NAG_Fortran_DLL_info.exe which is available from the Start Menu shortcut
NAG Fortran Library DLLs. Check AccessibilitySee Section 4.2 for details of this utility.
See Section 4.1.1.1 below for information on setting environment variables from a command prompt. The PATH, LIB and INCLUDE environment variables may already have been set globally as part of the installation or this may be done via the Control Panel. (On Windows XP, from Control Panel select System | Advanced | Environment Variables.) Either the user or the system variables may be edited, although Administrator privileges will be required to edit the system ones. Edit the PATH environment variable to include
c:\Program Files\NAG\FL21\fldll214zl\batch; c:\Program Files\NAG\FL21\fldll214zl\bin; c:\Program Files\NAG\FL21\fldll214zl\mkl_ia32_8.0\bin; existing pathadd or edit the LIB environment variable to include
c:\Program Files\NAG\FL21\fldll214zl\lib; c:\Program Files\NAG\FL21\fldll214zl\mkl_ia32_8.0\lib; any existing library pathadd or edit the INCLUDE environment variable to include
c:\Program Files\NAG\FL21\fldll214zl\nag_interface_blocks; any existing include pathsubstituting the correct folder where the NAG Fortran DLLs are installed if necessary.
In this DLL implementation, for convenience, the MKL symbols are exported directly from the NAG import library FLDLL214Z_mkl.lib, so it is not necessary to specify the MKL import libraries libguide40.lib and mkl_s_dll.lib as well. However, if the MKL import libraries are specified, it is important that the NAG import library preceeds them, i.e. the order should be
FLDLL214Z_mkl.lib libguide40.lib mkl_s_dll.libbecause certain parts of the MKL should not be used (see Section 5).
Information on calling the NAG Fortran DLLs from various different environments is given below. More information on calling NAG Fortran or C DLLs may also be available on the NAG web site at
http://www.nag.co.uk/numeric/Num_DLLhelp.asp
The NAG Fortran DLLs have been built using Intel Fortran 8.1. To call the DLLs from a program compiled with Intel Fortran 9, you may need to move or rename the file libmmd.dll in the install dir\bin folder, so that the correct Intel Fortran runtime DLL is picked up.
4.1.1.1. From a command window
To access this implementation from a command window some environment variables need to be set (if this has not been done globally - see above).
The shortcut:
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs. Command Prompt
may be used to start a command prompt window with the correct settings for the INCLUDE, LIB and PATH environment variables for the Library and the supplied MKL.
If the shortcut is not used, you can set the environment variables by running the batch file envvars.bat for this implementation. The default location of this file is:
c:\Program Files\NAG\FL21\fldll214zl\batch\envvars.batIf the file is not in the default location, you can locate it by searching for the file envvars.bat containing fldll214zl.
You may then compile and link to the NAG Fortran Library on the command line using one of the following commands:
ifort /iface:cvf /libs:dll driver.f FLDLL214Z_mkl.lib ifort /iface:cvf /libs:dll driver.f FLDLL214Z_nag.libwhere driver.f is your application program.
The first command will use the library without the NAG version of the BLAS/LAPACK procedures (FLDLL214Z_mkl.dll) with the MKL libraries mkl_*.dll and libguide40.dll.
The second command will use the library with the NAG version of the BLAS/LAPACK procedures (FLDLL214Z_nag.lib).
Please note that the Intel Visual Fortran compiler environment variables must be set in the command window. For more details refer to the User's Guide for the compiler.
4.1.1.2. From MS Visual Studio .NET
The instructions given here are for Visual Studio .NET 2003 with Intel Fortran Compiler 8.1. Other versions may vary.
To ensure that the NAG DLLs are accessible at runtime, the PATH environment variable must be set such that the location of the NAG Fortran DLLs, specifically the folder install dir\bin, is on the path. If the MKL version of the DLL is required, the location of the MKL DLLs, install dir\mkl_ia32_8.0\bin must also be on the path, but should appear after the install dir\bin folder.
Once Visual Studio has been opened, it is possible to set up the folders for use with Intel Fortran in this and all subsequent projects which use this compiler. One way to do so is:
The interface block is simply accessed by inserting a USE statement as described in greater detail in Section 4.6.
Whilst the above changes will apply to every Intel Fortran project, the following tasks need to be performed for each individual Intel Fortran project.
The Library is designed to be run in fully optimised mode, so to avoid any warning messages, you might decide to set the active configuration to Release. You can do this from the Toolbar or alternatively via the Build|Configuration Manager menus. If you work in Debug mode, you may receive a warning message about conflicting C runtime libraries. This is usually harmless.
There are a number of ways to add the NAG Library to the project. We describe just two; choose the one that most suits you.
If the Solution Explorer window is open then make sure that group project (the first line) is NOT selected. From the Project menu, choose the project Properties item. (Alternatively right click on a specific single project in the Solution Explorer and choose Properties.)
From the form, click Linker in the leftmost panel and then choose Input. The right hand panel will now have an Additional Dependencies entry, and you need to type FLDLL214Z_mkl.lib in this location to use the FLDLL214Z_mkl.lib library and MKL. If you wish to use the self-contained NAG Library then you need to add FLDLL214Z_nag.lib instead.
You will also need to set the calling convention to "CVF". From the Properties form, click Fortran and then choose External Procedures. Click on the Calling Convention entry in the right hand panel and select CVF from the drop-down list. Selecting this option will automatically change the entry for String Length Argument Passing to After Individual String Argument when you click on Apply or OK; this is the convention required by this implementation of the NAG Library.
We also recommend that you use the single-threaded DLL version of the runtime library. To select this option, from the Properties form, click Fortran and then choose Libraries. Click on the Runtime Library entry in the right hand panel and select Single-threaded DLL from the drop-down list.
Press the OK button to accept the changes and close the form.
The Properties information may also be accessed via the Toolbar. With the project selected in Solution Explorer, choose the Properties Window button on the Toolbar. In the ensuing window choose then the rightmost Property Pages icon. As in the paragraphs above, from the form click Linker in the leftmost panel and then choose Input. The right hand panel will now have an Additional Dependencies entry, and you need to type the names of the relevant libraries in this location. Under the Fortran section, click on External Procedures and select CVF for the calling convention, and click on Libraries and select the Single-threaded DLL runtime library, as described above. Press the OK button.
The project should now compile and link using the appropriate choice from the Build menu.
To run a program that does not require input or output redirected from standard input or standard output, from within the Microsoft Development Environment, the program may be executed via the Debug menu (by selecting Start Without Debugging, for example).
If a data file needs to be attached to the standard input or the output of a program needs to be redirected to the standard output, we recommend that you run the executable from a command prompt window to avoid the limitations of the current version of Visual Studio .NET.
Assuming that the folder containing the libraries has been added to the LIB environment variable, you may compile and link to the NAG Fortran Library on the command line in the following manner:
f90 driver.f FLDLL214Z_mkl.libor
f90 driver.f FLDLL214Z_nag.libwhere driver.f is your application program.
f95 -ieee=full -f77 -w=obs -o driver.exe driver.f "install dir\lib\FLDLL214Z_mkl.lib"or
f95 -f77 -w=obs -o driver.exe driver.f "install dir\lib\FLDLL214Z_nag.lib"where driver.f is your application program and driver.exe is the executable produced. Note the use of the -ieee=full option with FLDLL214Z_mkl.lib - see Section 3 above.
The full pathname of the FLDLL214Z_mkl.lib or FLDLL214Z_nag.lib file must be specified and must be enclosed within quotes if it contains spaces.
install dir\samples\absoft_fortran_examples
There are four issues to consider when using the NAG library
with Absoft Pro Fortran:
VAL(LOC(RELABS)), VAL(7)(see lines marked "CHANGE2" in absoft.for), where RELABS is declared and assigned the value 'Default' in the lines marked "NEW" in that file)
Assuming that the folder containing the libraries has been added to the LIB environment variable, you may compile and link to the NAG Fortran Library on the command line using one of the following commands:
f77 driver.for FLDLL214Z_mkl.lib f77 driver.for FLDLL214Z_nag.lib f95 driver.f FLDLL214Z_mkl.lib f95 driver.f FLDLL214Z_nag.libwhere driver.for or driver.f is your application program.
No source file changes are necessary to call the DLLs from FTN95. However, since FTN95 uses a variant of the cdecl calling convention, the compiler has to be told that the routines in the DLLs are to be called using the CVF calling convention. This can be accomplished using the /IMPORT_LIB command line switch as follows:
ftn95 driver.f /import_lib "install dir\bin\FLDLL214Z_mkl.dll"or
ftn95 driver.f /import_lib "install dir\bin\FLDLL214Z_nag.dll"The full pathname of install dir should be specified to the DLLs and should be enclosed within quotes if it contains spaces. The effect of this is to assume that all exported names in the DLL are CVF STDCALL and that any use of them should use the CVF STDCALL calling convention. External names passed via the argument list to a routine in a NAG DLL are automatically adjusted for whether or not they occur in the same source.
The NAG libraries may be added to the SLINK command line as usual e.g.
slink driver.obj "install dir\bin\FLDLL214Z_mkl.dll"or
slink driver.obj "install dir\bin\FLDLL214Z_nag.dll"As with compilation, the full path to the DLLs should be specified here, within quotes if the pathname contains spaces. It is worth emphasising that the linker should link directly against the DLLs, not the *.lib files.
It is also possible to compile and link using commands such as
ftn95 /f_stdcall driver.f slink driver.obj FLDLL214Z_nag.dll
Plato3 is the Integrated Development Environment (IDE) that is provided with the more recent versions of FTN95. To use Plato3 for a project involving a NAG routine:
open(6,file='c:\test.res')
in the main program before any write statements to channel 6.
Assuming that the LIB and PATH environment variables have been set up appropriately for your installation of the NAG Library, the command for linking to the Mark 21 DLLs using pgf77 is:
pgf77 -driver.f FLDLL214Z_nag.lib -o driver.exe(for the self-contained variant of the Library), or
pgf77 driver.f FLDLL214Z_mkl.lib -o driver.exe(for the MKL-based variant of the Library). The pgf90 compiler may be used instead of pgf77 in the same way.
This has been tested using version 6.0-8 of the PGI compilers under Windows XP.
install dir\samples\lahey_fortran_exampleThe advice given here has been tested using Lahey Fortran version 7.1.
There are a few simple changes that must be made to a standard Fortran program to allow the NAG DLLs to be used by Lahey Fortran:
DLL_IMPORT D02CJF
DLL_IMPORT D02CJW
DLL_IMPORT OUT
and in the user-supplied subroutine or function insert a DLL_EXPORT
statement i.e.
DLL_EXPORT OUT
should be inserted in subroutine OUT.
Remember to declare all these subroutines and functions as EXTERNAL and also to declare the type of any functions used. The names of imported functions are case sensitive; this means that NAG names must be in upper case.
The second family of changes concern the treatment of character arguments which must be adjusted to suit the convention used by the NAG DLLs. Character arguments must be stripped of the hidden length argument that Lahey places at the end of all the arguments; this is accomplished by passing the value of the address of the argument as follows: VAL(POINTER(char_arg)). Then, to conform to the NAG DLL standard, the length argument needs to be added immediately following the character argument. Both arguments are passed by value. Here is an example:
CALL D02CJF (X, XEND2, N, Y, FCN, TOL, VAL(POINTER('Default')),
+ VAL(LEN('Default')), OUT, G, W, IFAIL)
Use the compiler switch "-ml bc" to compile all routines in the program. The Lahey compiler uses this switch to specify that the stdcall calling convention be used. The import library file for the DLL should be in a location specified in the LIB environment variable. Alternatively, the programmer can specify the location using the -LIBPATH linker option. Specifying the location of the library on the compiler line is a third possibility i.e.
lf95 d02cjfe.f "install dir\lib\FLDLL214Z_mkl.lib" -ml bc
or
lf95 d02cjfe.f "install dir\lib\FLDLL214Z_nag.lib" -ml bc
As ever, make sure that the DLLs are on the PATH.
Examples of the use of the DLLs from C and C++ are given in the install dir\samples\c_examples and install dir\samples\cpp_examples folders. (Note that if the e04ccfcppclass.sln file is loaded into Visual Studio 2005, the Conversion Wizard will update the files in the project/solution as necessary.)
A document, techdoc.html, giving more detailed advice on calling the DLLs from C is available in install dir\headers. There is also a shortcut to this document on the Start Menu under
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs.
C & C++ Header Files Information
by default.
Key information:
cl driver.c FLDLL214Z_mkl.libor
cl driver.c FLDLL214Z_nag.libwhere driver.c is your application program. This assumes that the folder containing the header files has been added to the INCLUDE environment variable. If not, you could use:
cl driver.c FLDLL214Z_mkl.lib /I"install dir\headers"or
cl driver.c FLDLL214Z_nag.lib /I"install dir\headers"
The following instuctions apply to Visual Studio .NET 2003 and Visual Studio 2005. Later versions may vary.
If you are working under the Visual Studio IDE, set the following values to enable linking to work. Under the project's Properties, select Configuration Properties | Linker | Input and add FLDLL214Z_nag.lib or FLDLL214Z_mkl.lib to the Additional Dependencies field. If the LIB environment variable has not been set elsewhere, select Configuration Properties | Linker | General and add install dir\lib to the Additional Library Directories field.
Assuming that the folder containing the libraries has been added to the LIB environment variable, you may compile and link your C application program to the NAG Fortran Library on the command line in the following manner:
icl driver.c FLDLL214Z_mkl.lib /I"install dir\headers"or
icl driver.c FLDLL214Z_nag.lib /I"install dir\headers"where driver.c is your application program.
Assuming that the folder containing the libraries has been added to the LIB environment variable, you may compile and link your C application program to the NAG Fortran Library on the command line in the following manner:
acc driver.c FLDLL214Z_mkl.libor
acc driver.c FLDLL214Z_nag.libwhere driver.c is your application program.
impdef name.def "install dir\bin\name.dll" implib name.lib name.defwhere name denotes the name of the NAG DLL e.g. FLDLL214Z_nag. The first statement constructs a module definition file, name.def, and the second takes this module definition file and constructs an import library, name.lib. Do not be alarmed by warning messages from IMPLIB. These arise from the number of alternative symbols exported from the DLL in order to provide convenient alternatives for different users.
Assuming that the folder containing the imported libraries has been added to the LIB environment variable, you may compile and link your C application program to the NAG Fortran Library on the command line in the following manner:
bcc32 -I"install dir\headers" driver.c name.libwhere driver.c is your application program.
Alternatively you may add the location of the NAG header files to the configuration file bcc32.cfg. For more details please see the compiler documentation. If you have amended the configuration file then you may simply type:
bcc32 driver.c name.lib
In order to avoid the MKL floating-point exception problems mentioned in Section 3, we recommend that you use the self-contained NAG DLL (FLDLL214Z_nag.dll) when calling the NAG Library from the Borland C/C++ compiler.
Examples of use of the DLLs from within Excel are given in the install dir\samples\excel_examples folder. The folder install dir\samples\excel_examples\linear_algebra contains the file xls_demo.txt. This file gives some hints about using NAG DLLs from within Excel spreadsheets.
Key information:
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs.
VB6 Declare Statements
by default.
Examples of use of the DLLs from Visual Basic are given in the install dir\samples\vb_examples folder.
Key information:
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs.
VB6 Declare Statements
by default.
Key information:
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs.
VB.NET Declare Statements
by default.
Examples of use of the DLLs from Delphi are given in the install dir\samples\delphi_examples folder. These have been tested with Delphi 2006 for Microsoft Win32. The file readme.txt in the install dir\samples\delphi\e04ucf\console folder indicates how to use the DLLs from Delphi in a console window.
Key information:
In order to avoid the MKL floating-point exception problems mentioned in Section 3, we recommend that you use the self-contained NAG DLL (FLDLL214Z_nag.dll) when calling the NAG Library from Delphi.
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs. Check Accessibilityshortcut on the Start Menu, the DLLs will be found if the global PATH environment variable has been set to include their location. If it is run from a Command Prompt window, then the environment local to that window will prevail. (See Section 4.1 for more information on environment variables.)
For each of the DLLs that it loads successfully, NAG_Fortran_DLL_info.exe will print the location of the DLL, and the implementation details as provided by a call to the routine A00AAF. It also calls A00ACF to check the availability of a valid license key if required by the implementation and reports if this is not found; "Status OK" indicates that a valid key was either found or not required.
The example programs are most easily accessed by the batch files nag_example.bat or nag_example_mkl.bat.
The batch files need the environment variable NAG_FLDLL214ZL.
As mentioned in Section 4.1.1.1, the installation procedure provides a shortcut which starts a Command Prompt with local environment variables. The environment variables include NAG_FLDLL214ZL. This shortcut is placed in the Start Menu under
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs. Command PromptIf the shortcut is not used, you need to set this environment variable. You can set this environment variable by running the batch file envvars.bat for this implementation. The default location of this file is:
c:\Program Files\NAG\FL21\fldll214zl\batch\envvars.batIf the file is not in the default location, you can locate it by searching for the file envvars.bat containing fldll214zl.
nag_example_mkl.bat will provide you with a copy of an example program (and its data, if any), compile the program and link it with the library FLDLL214Z_mkl.lib and the MKL. Finally, the executable program will be run. The example program concerned is specified by the argument to nag_example_mkl.bat, e.g.
nag_example_mkl e04ucfwill copy the example program and its data into the files e04ucfe.f and e04ucfe.d in the current folder and process them to produce the example program results in the file e04ucfe.r.
Alternatively you could use:
nag_example e04ucf
The difference between nag_example_mkl.bat and nag_example.bat is that while nag_example_mkl.bat uses the library FLDLL214Z_mkl.lib and the MKL libraries, nag_example.bat uses the self-contained library FLDLL214Z_nag.lib. The executables created will use FLDLL214Z_mkl.dll and FLDLL214Z_nag.dll respectively.
The example programs are supplied in machine-readable form. They are suitable for immediate usage. Note that the distributed example programs are those used in this implementation and may not correspond exactly with the programs published in the Library Manual. The distributed example programs should be used in preference wherever possible.
The distributed example results are those obtained with the NAG self-contained library FLDLL214Z_nag.dll, (using the NAG BLAS and LAPACK routines). Running the examples with MKL BLAS or LAPACK may give slightly different results.
real or double precision - DOUBLE PRECISION (REAL*8)
basic precision - double precision
complex or complex*16 - COMPLEX*16
additional precision - quadruple precision (REAL*16,COMPLEX*32)
machine precision - the machine precision, see the value
returned by X02AJF in Section 5
Thus a parameter described as
real or double precision
should be declared as DOUBLE PRECISION in your program. If a
routine accumulates an inner product in additional precision,
it is using software to simulate quadruple precision.
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.
(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) compiled with the Intel Fortran Compiler 8.1, for use with that compiler. The interface blocks are also supplied in source form (.f90 files), which may be compiled for use with other compilers.
If you use the Library command prompt shortcut or set the environment variables by running the batch file envvars.bat for this implementation (see Section 4.1.1.1), you can use any of the commands described in Section 4.1.1.1 to access these modules since the environment variable INCLUDE will be set.
In order to make use of these modules from existing Fortran 77 code, the following changes need to be made:
The above steps need to be done for each unit (main program, function or subroutine) in your code.
These changes are illustrated by showing the conversion of the Fortran 77 version of the example program for NAG Fortran Library routine D01DAF. Please note that this is not exactly the same as the example program that is distributed with this implementation. Each change is surrounded by comments boxed with asterisks.
* D01DAF Example Program Text
*****************************************************
* Add USE statements for relevant chapters *
USE NAG_F77_D_CHAPTER
* *
*****************************************************
* .. Parameters ..
INTEGER NOUT
PARAMETER (NOUT=6)
* .. Local Scalars ..
DOUBLE PRECISION ABSACC, ANS, YA, YB
INTEGER IFAIL, NPTS
* .. External Functions ..
DOUBLE PRECISION FA, FB, P1, P2A, P2B
EXTERNAL FA, FB, P1, P2A, P2B
* .. External Subroutines ..
******************************************************
* EXTERNAL declarations need to be removed. *
C EXTERNAL D01DAF
* *
******************************************************
* .. Executable Statements ..
WRITE (NOUT,*) 'D01DAF Example Program Results'
YA = 0.0D0
YB = 1.0D0
ABSACC = 1.0D-6
WRITE (NOUT,*)
WRITE (NOUT,*) 'First formulation'
IFAIL = 1
*
CALL D01DAF(YA,YB,P1,P2A,FA,ABSACC,ANS,NPTS,IFAIL)
*
WRITE (NOUT,99999) 'Integral =', ANS
WRITE (NOUT,99998) 'Number of function evaluations =', NPTS
IF (IFAIL.GT.0) WRITE (NOUT,99997) 'IFAIL = ', IFAIL
WRITE (NOUT,*)
WRITE (NOUT,*) 'Second formulation'
IFAIL = 1
*
CALL D01DAF(YA,YB,P1,P2B,FB,ABSACC,ANS,NPTS,IFAIL)
*
WRITE (NOUT,99999) 'Integral =', ANS
WRITE (NOUT,99998) 'Number of function evaluations =', NPTS
IF (IFAIL.GT.0) WRITE (NOUT,99997) 'IFAIL = ', IFAIL
STOP
*
99999 FORMAT (1X,A,F9.4)
99998 FORMAT (1X,A,I5)
99997 FORMAT (1X,A,I2)
END
*
DOUBLE PRECISION FUNCTION P1(Y)
* .. Scalar Arguments ..
DOUBLE PRECISION Y
* .. Executable Statements ..
P1 = 0.0D0
RETURN
END
*
DOUBLE PRECISION FUNCTION P2A(Y)
* .. Scalar Arguments ..
DOUBLE PRECISION Y
* .. Intrinsic Functions ..
INTRINSIC SQRT
* .. Executable Statements ..
P2A = SQRT(1.0D0-Y*Y)
RETURN
END
*
DOUBLE PRECISION FUNCTION FA(X,Y)
* .. Scalar Arguments ..
DOUBLE PRECISION X, Y
* .. Executable Statements ..
FA = X + Y
RETURN
END
*
DOUBLE PRECISION FUNCTION P2B(Y)
*****************************************************
* Add USE statements for relevant chapters *
USE NAG_F77_X_CHAPTER
* *
*****************************************************
* .. Scalar Arguments ..
DOUBLE PRECISION Y
* .. External Functions ..
******************************************************
* Function Type declarations need to be removed. *
C DOUBLE PRECISION X01AAF
* *
******************************************************
******************************************************
* EXTERNAL declarations need to be removed. *
C EXTERNAL X01AAF
* *
******************************************************
* .. Executable Statements ..
P2B = 0.5D0*X01AAF(0.0D0)
RETURN
END
*
DOUBLE PRECISION FUNCTION FB(X,Y)
* .. Scalar Arguments ..
DOUBLE PRECISION X, Y
* .. Intrinsic Functions ..
INTRINSIC COS, SIN
* .. Executable Statements ..
FB = Y*Y*(COS(X)+SIN(X))
RETURN
END
In this implementation calls to the NAG version of the following Basic Linear Algebra Subprograms (BLAS) and linear algebra routines (LAPACK) are included in the library FLDLL214Z_mkl.lib to avoid problems with the vendor version:
DBDSQR ZBDSQR ZGEEVX
S07AAF F(1) = 1.0D+13
F(2) = 1.0D-14
S10AAF E(1) = 1.8500D+1
S10ABF E(1) = 7.080D+2
S10ACF E(1) = 7.080D+2
S13AAF x(hi) = 7.083D+2
S13ACF x(hi) = 1.0D+16
S13ADF x(hi) = 1.0D+17
S14AAF IFAIL = 1 if X > 1.70D+2
IFAIL = 2 if X < -1.70D+2
IFAIL = 3 if abs(X) < 2.23D-308
S14ABF IFAIL = 2 if X > 2.55D+305
S15ADF x(hi) = 2.66D+1
x(low) = -6.25D+0
S15AEF x(hi) = 6.25D+0
S17ACF IFAIL = 1 if X > 1.0D+16
S17ADF IFAIL = 1 if X > 1.0D+16
IFAIL = 3 if 0.0D+00 < X <= 2.23D-308
S17AEF IFAIL = 1 if abs(X) > 1.0D+16
S17AFF IFAIL = 1 if abs(X) > 1.0D+16
S17AGF IFAIL = 1 if X > 1.038D+2
IFAIL = 2 if X < -5.6D+10
S17AHF IFAIL = 1 if X > 1.041D+2
IFAIL = 2 if X < -5.6D+10
S17AJF IFAIL = 1 if X > 1.041D+2
IFAIL = 2 if X < -1.8D+9
S17AKF IFAIL = 1 if X > 1.041D+2
IFAIL = 2 if X < -1.8D+9
S17DCF IFAIL = 2 if abs (Z) < 3.93D-305
IFAIL = 4 if abs (Z) or FNU+N-1 > 3.27D+4
IFAIL = 5 if abs (Z) or FNU+N-1 > 1.07D+9
S17DEF IFAIL = 2 if imag (Z) > 7.00D+2
IFAIL = 3 if abs (Z) or FNU+N-1 > 3.27D+4
IFAIL = 4 if abs (Z) or FNU+N-1 > 1.07D+9
S17DGF IFAIL = 3 if abs (Z) > 1.02D+3
IFAIL = 4 if abs (Z) > 1.04D+6
S17DHF IFAIL = 3 if abs (Z) > 1.02D+3
IFAIL = 4 if abs (Z) > 1.04D+6
S17DLF IFAIL = 2 if abs (Z) < 3.93D-305
IFAIL = 4 if abs (Z) or FNU+N-1 > 3.27D+4
IFAIL = 5 if abs (Z) or FNU+N-1 > 1.07D+9
S18ADF IFAIL = 2 if 0.0D+00 < X <= 2.23D-308
S18AEF IFAIL = 1 if abs(X) > 7.116D+2
S18AFF IFAIL = 1 if abs(X) > 7.116D+2
S18CDF IFAIL = 2 if 0.0D+00 < X <= 2.23D-308
S18DCF IFAIL = 2 if abs (Z) < 3.93D-305
IFAIL = 4 if abs (Z) or FNU+N-1 > 3.27D+4
IFAIL = 5 if abs (Z) or FNU+N-1 > 1.07D+9
S18DEF IFAIL = 2 if real (Z) > 7.00D+2
IFAIL = 3 if abs (Z) or FNU+N-1 > 3.27D+4
IFAIL = 4 if abs (Z) or FNU+N-1 > 1.07D+9
S19AAF IFAIL = 1 if abs(x) >= 4.95000D+1
S19ABF IFAIL = 1 if abs(x) >= 4.95000D+1
S19ACF IFAIL = 1 if X > 9.9726D+2
S19ADF IFAIL = 1 if X > 9.9726D+2
S21BCF IFAIL = 3 if an argument < 1.579D-205
IFAIL = 4 if an argument >= 3.774D+202
S21BDF IFAIL = 3 if an argument < 2.820D-103
IFAIL = 4 if an argument >= 1.404D+102
X01AAF (PI) = 3.1415926535897932D+00 X01ABF (GAMMA) = 0.5772156649015329D+00
The basic parameters of the model
X02BHF = 2 X02BJF = 53 X02BKF = -1021 X02BLF = 1024 X02DJF = .TRUE.
Derived parameters of the floating-point arithmetic
X02AJF = 1.11022302462516D-16 X02AKF = 2.22507385850721D-308 X02ALF = 1.79769313486231D+308 X02AMF = 2.22507385850721D-308 X02ANF = 2.22507385850721D-308
Parameters of other aspects of the computing environment
X02AHF = 1.42724769270596D+45 X02BBF = 2147483647 X02BEF = 15 X02DAF = .FALSE.
The routine X03AAW changes the floating-point control word such that it sets the working precision to double precision (53-bit mantissa) and sets the rounding mode to nearest. X03AAW takes a single INTEGER argument, which is used to store the floating-point control word value on input and return it to the calling program, so this routine is also used to retrieve the original (i.e. on entry) value of the control word.
The routine X03AAX sets the floating-point control word to the value specified in its single INTEGER argument. It is typically used to restore the floating-point control word to its original value after a call to X03AAW, but may, of course, be used to set a different value.
The Library is designed to operate in double precision (53-bit) mode, not the extended precision mode also possible on the chip. A normal Intel Fortran program will operate, by default, in this mode, but other environments may re-set the floating-point control word so that the chip operates in extended precision mode. Excel is one such environment. To obtain consistent behaviour the user may wish to use X03AAW directly before entering any other Library routine in order to restore the mode of operation to double precision.
A corollary of the behaviour of these two routines is that they may also be used as a "get and set" pair, with X03AAW used to get the value of the floating-point control word on entry, and X03AAX used to set a new value (or restore the original value), but note the "side-effect" of calling X03AAW, namely that this will set the control word as described above.
The Library Manual is supplied in the form of Portable Document Format (PDF) files, with an HTML index, in the manual folder (either installed locally or on the distribution CD). The introductory material is also provided as HTML files in the manual folder.
A main index file has been provided (manual\html\mark21.html) which contains a fully linked contents document pointing to all the available PDF (and where available HTML) files. This index file is available from the Start Menu under
Start|All Programs|NAG|FL21|NAG Fortran Library Manualby default. Use your HTML browser to navigate from here.
In addition the following are provided:
Start|All Programs|NAG|FL21|NAG Fortran Library DLLs. Users' Noteby default.
The Library Manual is also available as an HTML Help file. In order to display the mathematical expressions correctly, you will need to install Design Science's MathPlayer plug-in. MathPlayer may optionally be installed with this NAG Library product, or may be downloaded from Design Science's web site:
http://www.dessci.com/en/products/mathplayer/download.htm
The HTML Help version of the Library Manual is available from the Start Menu under
Start|All Programs|NAG|FL21|NAG Fortran Library Manual HTML HelpIf the NAG Library materials have been installed on a network drive, you may need to copy the HTML Help file (nagdoc_fl21.chm) to a local drive. If you still have trouble seeing the correct help file materials, e.g. if you see a message such as
Navigation to the webpage was canceledthe file might have been blocked by security updates in Windows or Internet Explorer 7. Right-click on nagdoc_fl21.chm in Windows Explorer and select Properties from the pop-up menu. If there is an Unblock button at the bottom click on it to unblock the file and then click on OK to close the Properties dialog box.
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 FLDLL214ZL).
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 the following URLs:
http://www.nag.co.uk/, http://www.nag.com/ or http://www.nag-j.co.jp/
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: support@nag.com Tel: +1 630 971 2337 Tel: +1 630 971 2337 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