* integration
2002-12-31 9:55 integration Francisco J Molina
@ 2002-07-03 22:04 ` Francisco J Molina
2002-12-31 9:55 ` integration Brian Gough
1 sibling, 0 replies; 11+ messages in thread
From: Francisco J Molina @ 2002-07-03 22:04 UTC (permalink / raw)
To: gsl-discuss
I am using gsl_integration_qagi to integrate a function similar to a
Gaussian ( it is always positive and tends to zero as x tends to +/-
infinity.
Since this function is never 0, I wonder if it is possible to modify the
source code to use a 60-point Gauss-Kronrod rule of QAGS?
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
2002-12-31 9:55 ` integration Brian Gough
@ 2002-07-04 15:07 ` Brian Gough
0 siblings, 0 replies; 11+ messages in thread
From: Brian Gough @ 2002-07-04 15:07 UTC (permalink / raw)
To: Francisco J Molina; +Cc: gsl-discuss
Francisco J Molina writes:
> I am using gsl_integration_qagi to integrate a function similar to a
> Gaussian ( it is always positive and tends to zero as x tends to +/-
> infinity.
>
> Since this function is never 0, I wonder if it is possible to modify the
> source code to use a 60-point Gauss-Kronrod rule of QAGS?
The integration rule is fixed at compile time, but if you copy the
file integration/qags.c and its #include's into your application you
can modify the reference to the integration rule and override the
version in the library.
Brian
^ permalink raw reply [flat|nested] 11+ messages in thread
* integration
@ 2002-12-31 9:55 Francisco J Molina
2002-07-03 22:04 ` integration Francisco J Molina
2002-12-31 9:55 ` integration Brian Gough
0 siblings, 2 replies; 11+ messages in thread
From: Francisco J Molina @ 2002-12-31 9:55 UTC (permalink / raw)
To: gsl-discuss
I am using gsl_integration_qagi to integrate a function similar to a
Gaussian ( it is always positive and tends to zero as x tends to +/-
infinity.
Since this function is never 0, I wonder if it is possible to modify the
source code to use a 60-point Gauss-Kronrod rule of QAGS?
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
2002-12-31 9:55 integration Francisco J Molina
2002-07-03 22:04 ` integration Francisco J Molina
@ 2002-12-31 9:55 ` Brian Gough
2002-07-04 15:07 ` integration Brian Gough
1 sibling, 1 reply; 11+ messages in thread
From: Brian Gough @ 2002-12-31 9:55 UTC (permalink / raw)
To: Francisco J Molina; +Cc: gsl-discuss
Francisco J Molina writes:
> I am using gsl_integration_qagi to integrate a function similar to a
> Gaussian ( it is always positive and tends to zero as x tends to +/-
> infinity.
>
> Since this function is never 0, I wonder if it is possible to modify the
> source code to use a 60-point Gauss-Kronrod rule of QAGS?
The integration rule is fixed at compile time, but if you copy the
file integration/qags.c and its #include's into your application you
can modify the reference to the integration rule and override the
version in the library.
Brian
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
2003-08-20 11:07 ` integration Brian Gough
@ 2003-08-20 11:52 ` Gert Van den Eynde
0 siblings, 0 replies; 11+ messages in thread
From: Gert Van den Eynde @ 2003-08-20 11:52 UTC (permalink / raw)
To: Brian Gough, Axel Hutt; +Cc: gsl-discuss
On Wednesday 20 August 2003 13:07, Brian Gough wrote:
> The only thing I can contribute here is that all of the quadpack
> routines assume smooth functions. Apart from that, I don't know much
> about integro-differential equations.
A small addendum: QUADPACK (the original FORTRAN code) allows functions with
end-point singularities (routine DQAGS(E)). You would do that anyway when
integrating a function with singularties: locate the singularities and split
up the integration interval accordingly. With the routine DQAGP(E) you can
even provide information on the kind of singularty.I suppose the GSL routine
is the C translation of DQAGS. Of course, the function inside the integration
boundaries should be smooth (numerical integration is for a big part based on
the approximation of the function by orthogonal polynomials).
For an overview of all original QUADPACK routines, have a look here
http://netlib.org/quadpack/doc
Hope this helps,
gert
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
2003-08-15 9:45 integration Axel Hutt
@ 2003-08-20 11:07 ` Brian Gough
2003-08-20 11:52 ` integration Gert Van den Eynde
0 siblings, 1 reply; 11+ messages in thread
From: Brian Gough @ 2003-08-20 11:07 UTC (permalink / raw)
To: Axel Hutt; +Cc: gsl-discuss
Axel Hutt writes:
> I apply the GSL-routine qag and related for the integration part in
> integro-differential equations. The routines are okay, but quite
> slow for rather non-smooth data (like random initial values). Does
> anybody also apply GSL-routines for integration for
> integro-differential equations and if yes, which routines do you
> apply? I also try some hand-made Taylor-expansion-routines, which
> however show less stable behaviour.
The only thing I can contribute here is that all of the quadpack
routines assume smooth functions. Apart from that, I don't know much
about integro-differential equations.
^ permalink raw reply [flat|nested] 11+ messages in thread
* integration
@ 2003-08-15 9:45 Axel Hutt
2003-08-20 11:07 ` integration Brian Gough
0 siblings, 1 reply; 11+ messages in thread
From: Axel Hutt @ 2003-08-15 9:45 UTC (permalink / raw)
To: gsl-discuss
Dear all,
I apply the GSL-routine qag and related for the integration
part in integro-differential equations. The routines are okay,
but quite slow for rather non-smooth data (like random initial
values).
Does anybody also apply GSL-routines for integration for
integro-differential equations and if yes, which routines do you
apply?
I also try some hand-made Taylor-expansion-routines, which however
show less stable behaviour.
Axel
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
1999-03-08 7:56 ` integration Aric Hagberg
@ 1999-03-09 12:43 ` Brian Gough
0 siblings, 0 replies; 11+ messages in thread
From: Brian Gough @ 1999-03-09 12:43 UTC (permalink / raw)
To: gsl-discuss
Good point, Aric, we should start with the straightforward methods. If
you'd like to write something along those lines that would be really
useful.
Brian
Aric Hagberg wrote:
>
> DASSL is quite a complicated piece of software. To properly
> implement an ODE solver like that - actually I believe that does
> differential+algebraic equations too - we need things
> like nonlinear equation solvers, finite differencing,
> linear equation solvers, etc.
>
> It is certainly possible to write an ODE solver for GSL
> right now using Runge-Kutta or some other explicit method without
> needing those things.
>
> Aric
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
1999-03-05 12:15 ` integration Brian Gough
@ 1999-03-08 7:56 ` Aric Hagberg
1999-03-09 12:43 ` integration Brian Gough
0 siblings, 1 reply; 11+ messages in thread
From: Aric Hagberg @ 1999-03-08 7:56 UTC (permalink / raw)
To: gsl-discuss; +Cc: d-dooling
DASSL is quite a complicated piece of software. To properly
implement an ODE solver like that - actually I believe that does
differential+algebraic equations too - we need things
like nonlinear equation solvers, finite differencing,
linear equation solvers, etc.
It is certainly possible to write an ODE solver for GSL
right now using Runge-Kutta or some other explicit method without
needing those things.
Aric
^ permalink raw reply [flat|nested] 11+ messages in thread
* Re: integration
1999-03-04 13:33 integration David J. Dooling
@ 1999-03-05 12:15 ` Brian Gough
1999-03-08 7:56 ` integration Aric Hagberg
0 siblings, 1 reply; 11+ messages in thread
From: Brian Gough @ 1999-03-05 12:15 UTC (permalink / raw)
To: d-dooling; +Cc: gsl-discuss
Hi,
I am working on a C implementation of QUADPACK for the integration (i.e.
quadrature) section. We haven't done any work on an ODE solvers
directory yet but something like DASSL is a candidate for that.
regards
Brian Gough
^ permalink raw reply [flat|nested] 11+ messages in thread
* integration
@ 1999-03-04 13:33 David J. Dooling
1999-03-05 12:15 ` integration Brian Gough
0 siblings, 1 reply; 11+ messages in thread
From: David J. Dooling @ 1999-03-04 13:33 UTC (permalink / raw)
To: gsl-discuss
Hello, I stumbled across the GSL today and started looking through the
manual. In the manual there were sections for integration, but they
contained no text. Therefore, I was wondering if anyone could answer
a couple questions.
What sorts of integration capabilities does the GSL plan to implement?
I would guess you are planning on implementing simple integration
schemes for functions, e. g., Runge-Kutta and slightly more
sophisticated techniques, and a brief look at the code suggests this
is accurate.
Are there any plans to implement more robust differential equation
solvers capable of solving systems of ODE's, stiff ODE's, or
differential/algebraic equations (like DASSL)?
Thanks.
DAVID.
~~~~~~
d-dooling@nwu.edu
David J. Dooling
Dept. of Chemical Engineering
Northwestern University
http://winnie.chem-eng.nwu.edu/students/dooling.html
^ permalink raw reply [flat|nested] 11+ messages in thread
end of thread, other threads:[~2003-08-20 11:52 UTC | newest]
Thread overview: 11+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2002-12-31 9:55 integration Francisco J Molina
2002-07-03 22:04 ` integration Francisco J Molina
2002-12-31 9:55 ` integration Brian Gough
2002-07-04 15:07 ` integration Brian Gough
-- strict thread matches above, loose matches on Subject: below --
2003-08-15 9:45 integration Axel Hutt
2003-08-20 11:07 ` integration Brian Gough
2003-08-20 11:52 ` integration Gert Van den Eynde
1999-03-04 13:33 integration David J. Dooling
1999-03-05 12:15 ` integration Brian Gough
1999-03-08 7:56 ` integration Aric Hagberg
1999-03-09 12:43 ` integration Brian Gough
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).