A small fraction of the PEAK code is copyrighted by Numerical Recipes (Press et al. 1992(, to the extent that I have not been allowed to make it available in digital form, as also mentioned in section 11.2 (item (a) on p. 344). Kind permission has been given, however, to make the copyrighted code available in printed form in the book (in appendix C, pp. 441), with this permission gratefully being acknowledged.
Copyrighted routines are contained in the following three auxiliary files:
nr_routines_v1.f90
* contains the four subroutines four1, gauleg, gaussj, plgndr2
* these four routines are modified versions of four1, gauleg, gaussj, plgndr
as contained in Numerical Recipes (Press et al. 1992)
* the file is available in printed form in appendix C.3 of the book (pp. 464-466)
kasahara_nms_1_aux.f90
* contains the two subroutines jacobi and eigsrt
* these two routines are taken (in essentially unmodified form)
from Numerical Recipes (Press et al. 1992)
* this file is not reprinted in the book since these routines have remained essentially unchanged
add_f_fft.f90
* contains the subroutines fft and four1
* the code for subroutine fft is available within pk_code_s8.f90
* the code for subroutine four1 is part of nr_routines_v1.f90
I do apologize to all of you for the inconveniences that are incurred in setting up PEAK as they arise due to these copyright restrictions.
Nevertheless, such an approach, given by referring to a few copyrighted routines, did appear to me to be a justifyable option for carrying out standard numerical auxiliary tasks (such as computing Fourier transforms, or eigenstructures and inverses of matrices) within PEAK. The inconveniences thus incurred appear justifyable, as sacrificing immediate digital availability of the code is in turn traded for permission to use widely-known, well-tested, and self-contained Numerical Recipes code for such standard numerical auxiliary tasks. Thank you for your understanding and patience.
Sunday, 03 March 2024