Intel® Fortran Compiler
Build applications that can scale for the future with optimized code designed for Intel® Xeon® and compatible processors.

Integration Library Routine

pas1359
Beginner
1,433 Views
Under VAX FORTRAN77 there was a library routine for integration named "HSPRBL". Is there a similar rountine that I can call out using Visual FORTRAN and what is it? Thanks.
0 Kudos
9 Replies
Steven_L_Intel1
Employee
1,433 Views
There was no such routine in the VAX Fortran product. I can't find any references to such a routine in the various libraries that I can think of. You may want to see if the IMSL libraries, in Visual Fortran Professional Edition, have something that will work for you.

Steve
0 Kudos
pas1359
Beginner
1,433 Views
Are there any Numerical Integration Routines as part of the Visual Fortran Library that I can call out in a Program or do I have to generate my own? The VAX which I still use on occasions, does interpret the "HSPRBL" call out in the Program as an Integration routine. I am asking the question, because I would like to transfer the VAX program to Visual Fortran.
0 Kudos
Steven_L_Intel1
Employee
1,433 Views
If you can call a HSPRBL routine on the VAX, it's because you're linking against some non-Digital/Compaq library that provides it. I suggest checking a link map to see where it comes from.

I suggest you look at the IMSL F90MP and FNL documentation at http://www.vni.com/products/imsl/functionality.html to see if there's a routine that meets your needs. I'm not a mathematician, so I don't quite know what to look for. (I saw some reference to "elliptical integration"..) If you find a routine in these libraries, it would be in Visual Fortran Professional Edition's IMSL libraries.

Steve
0 Kudos
pas1359
Beginner
1,433 Views
Thank you
0 Kudos
Jugoslav_Dujic
Valued Contributor II
1,433 Views
FWIW, searches on IMSL, Netlib and Google yielded zero results :-(.
0 Kudos
Intel_C_Intel
Employee
1,433 Views
I assume you mean numerical integration, or quadrature. If so, IMSL MATH library, chapter 4 will probably do.
0 Kudos
llynisa
Beginner
1,433 Views
If you want to integrate a function that is specified numerically by a set of unequally-spaced data points, the NAG library D01GAF routine uses a Gill & Miller third-order finite difference method. I have developed an equivalent version.

Alan
0 Kudos
gfthomas8
Novice
1,433 Views
Alan:

Any chance of me acquiring your implementation of the Gill-Miller approach?

TIA,
Gerry T.
gfthomas@sympatico.ca
0 Kudos
llynisa
Beginner
1,432 Views
Gerry,

I have sent you a copy of my Gill & Miller function Integdata(x,y,n,errest,ierr). If anyone else would like a copy, ask me at alan.cruttenden@baesystems.com

Alan
0 Kudos
Reply