- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
include
"link_f90_static.h"USE
QDAWO_INTUSE
UMACH_INTUSE
CONST_INTINTEGER
IWEIGH, NOUTREAL
A, ABS, B, ERRABS, ERREST, ERROR, &EXACT, F, OMEGA, PI,
RESULTINTRINSIC
ABSEXTERNAL
F! Get output unit number
CALL
UMACH (2, NOUT)! Set limits of integration
A = 0.0
B = 1.0
! Weight function = sin(10.*pi*x)
IWEIGH = 2
PI = CONST(PI)
OMEGA = 10.*PI
! Set error tolerances
ERRABS = 0.0
CALL
QDAWO (F, A, B, IWEIGH, OMEGA, RESULT, ERRABS=ERRABS, &ERREST=ERREST)
! Print results
EXACT = -0.1281316
ERROR =
ABS(RESULT-EXACT)WRITE
(NOUT,99999) RESULT, EXACT, ERREST, ERROR99999
FORMAT ( Computed =, F8.3, 13X, Exact =, F8.3, /, /, &Error estimate =, 1PE10.3, 6X, Error =, 1PE10.3)
END
!
REAL
FUNCTION F (X)REAL
XREAL
ALOGINTRINSIC
ALOGIF
(X .EQ. 0.) THENF = 0.0
ELSE
F =
ALOG(X)END IF
RETURN
END
Link Copied
5 Replies
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
this library has many error...but this is a library that i just copy form IMSL library, can anyone tell me what`s wrong, and what is the function of 99999
THX
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
i can fix it when use " instead of '....but i want to know what is the function of 99999
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Its just the label of the FORMAT statement refered to in the WRITE on the line above.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
You don't say what the errors are.
Please read carefully the section of the Building Applications manual on IMSL, in the Using Libraries chapter. Make sure that you have added the proper IMSL INCLUDE and LIB folders to the appropriate directory lists under Tools..Options..Intel Fortran.
Please read carefully the section of the Building Applications manual on IMSL, in the Using Libraries chapter. Make sure that you have added the proper IMSL INCLUDE and LIB folders to the appropriate directory lists under Tools..Options..Intel Fortran.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
....this library has many error?
I've used IMSL for several decades; I've never seen significant problems with it. If you're getting them, my first presumption (and it should be yours, too) is that there is something wrong with how the library is being used.
As said by Steve, you're going to have to be a bit more forthcoming with information about the nature of the errors. For example, are the errors in the compile step, the link step, or run time errors, e.g., log(x) is undefined in real numbers if x <= 0 and undefined, always, when x = 0, so trying to take log(-1) will cause a run time error.
The 99999 is just the label for the format statement.
I've used IMSL for several decades; I've never seen significant problems with it. If you're getting them, my first presumption (and it should be yours, too) is that there is something wrong with how the library is being used.
As said by Steve, you're going to have to be a bit more forthcoming with information about the nature of the errors. For example, are the errors in the compile step, the link step, or run time errors, e.g., log(x) is undefined in real numbers if x <= 0 and undefined, always, when x = 0, so trying to take log(-1) will cause a run time error.
The 99999 is just the label for the format statement.
![](/skins/images/7FC17B7B85029576C25F1E43CE255B51/responsive_peak/images/icon_anonymous_message.png)
Reply
Topic Options
- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page