This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
qawo-problem and a book project
- From: Harald Wiedemann <wiedemann dot harald at litef dot de>
- To: gsl-discuss at sources dot redhat dot com
- Date: Wed, 9 Jul 2003 10:33:50 +0200
- Subject: qawo-problem and a book project
Hallo,
I am writing a book on "Numerical Physics", which will be published next year
by Springer as german text book. Since I want the students
to be able to translate the programs and to play with them I want to include
free c/c++/Fortran-Compilers, graphic programs and free numerical
libraries. GSL will play a dominant role, since the programs are discussed
in the text in their GSL-versions. As far as I understand all this is
no license problem as long as I include the source code of GSL
in the attached CD, right?
Now to my question concerning qawo-integration. Attached is a program
kreisbl2.cpp which shall calculate the interference pattern behind a
circular aperture and it seems to work fine - producing also results
which are not obviously wrong... - until it crashs with a
"Speicherzugriffsfehler", i.e. "segmentation fault" (compilation with
different compilers / compiler options doesn't matter). I modified it to
make it more easier to investigate (notation for linux):
> kreisbl2 kreisbl2.par dummy.tst 259
gives 0.000519 0.456277
(kreisbl2.par containts some parameters, dummy.tst is the result file, and 259
is a parameter which originally was an loop index determing the position
where the intensity is calculated)
> kreisbl2 kreisbl2.par dummy.tst 260
gives segmentation fault.
First I checked whether the function to integrate gives nan (not a number)
at a specific value which could enforce such an error, but this seems not the
case. Though I can live with this problem (just some points less in the
result plot...) it is not nice and I would be appreciated if someone could
help.
Thanks in advance,
Harald Wiedemann