This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

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



Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]