CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Any numerical triple integration program is available in Fortran?

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   January 14, 1999, 02:11
Default Any numerical triple integration program is available in Fortran?
Posts: n/a
I have a complicated equation in terms of Sine and Cosine which contains a vector integration. However, the Mathematica software failed to perform the analytical integration. Is there any software other than Mathematica which is capable to perform a complicated analytical integeration? The program quad3d in 'Numerical recipes in Fortran' is not practical in this case. Because, this program takes at least 2 hours to perform one integration from -inf to +inf with our equation and we have to do at least 1 Lakh iterations during our computations. Is there any numerical program available to do a triple integration within a reasonable time in Fortran language? Let me know your suggestions and insights in this regard.

With regards,

Radhakrishnan. K
  Reply With Quote

Old   January 14, 1999, 18:19
Default Re: Any numerical triple integration program is available in Fortran?
Adrin Gharakhani
Posts: n/a
It's not clear from your description what you are integrating, more specifically, how many variable parameters are involved in your integral. In general, if the integrals are expensive (e.g. elliptic integrals) it's a good idea to create tables for them, which will cost you time just once; in the computations all you have to do is to look the tables up (plus simple interpolation, which are fast). If you have many parameters, then you can at least reduce their number by doing a hybrid table/analytic, table/numerical, etc.

Hope I haven't stated the obvious

Adrin Gharakhani
  Reply With Quote

Old   February 15, 1999, 13:37
Default Re: Any numerical triple integration program is available in Fortran?
Posts: n/a

Based on the very limited information you have provided about your problem, I can offer the following suggestions.

1. Since your function is composed of sines and cosines,

I would strongly urge you to attempt a solution through

the use of Fourier techniques. ( an FFT may be applicable

to your problem, and could give you an answer that is

thousands of times faster than a "brute force" approach).

2. SInce you seem to be integrating between (+/-) inf,

you need to introduce a change of variables so as to

integrate between finite limits. This is easily done.

3. If your boundary is of a complex shape, then your best

bet is to employ Gaussian integration techniques along

each direction. It is possible to decompose your volume

integral into surface integrals on the boundary. The

surface integrals in turn may be decomposed into line

integrals along each elementary surface defining your

volume. These line integrals are finally evaluated by

Gaussian integration. I would be happy to provide more

details of this method.
  Reply With Quote

Old   March 4, 1999, 02:03
Default Re: Any numerical triple integration program is available in Fortran?
Radhakrishnan. K
Posts: n/a
Hello Sir,

Many many thanks for your helpful response on my question. Meanwhile, I have solved my problem using some numerical methods.
  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Call Fortran program from OpenFoam tsjb00 OpenFOAM Running, Solving & CFD 6 April 22, 2013 15:26
CFX11 + Fortran compiler ? Mohan CFX 20 March 30, 2011 18:56
Fortran Program with FLUENT UDFs gautamcfd FLUENT 0 August 17, 2010 05:39
Numerical integration Damir Galeev Main CFD Forum 2 June 24, 2001 05:59
control plotting in fortran program zhe zhang Main CFD Forum 2 November 16, 1999 14:10

All times are GMT -4. The time now is 15:58.