Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeyq_v1_0.tar.gz(6149 Kbytes)
Manuscript Title: SPIP: A computer program implementing the Interaction Picture method for simulation of light-wave propagation in optical fibre
Authors: Stéphane Balac, Arnaud Fernandez
Program title: SPIP
Catalogue identifier: AEYQ_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 199(2016)139
Programming language: C.
Computer: Desktop computer.
Operating system: Linux, MS Windows.
RAM: 8 Giga bytes
Keywords: Generalized non-linear Schrödinger equation, Interaction Picture method, Embedded Runge-Kutta method, Adaptive step-size control.
Classification: 4.12, 18.

External routines: FFTW, a C subroutine library for computing the discrete Fourier transform, see [1] and http://www.fftw.org. Gnuplot, a portable command-line driven graphing utility, see [2] and http://www.gnuplot.info.

Nature of problem:
The program solves the Generalized Non-Linear Schrödinger Equation (GNLSE) which occurs in the field of non-linear optics as a model of wave propagation in fibre optics.

Solution method:
The GNLSE is solved by the Interaction Picture method coupled with an embedded Runge-Kutta scheme of order 4. The program includes a costless adaptive step-size control strategy taking advantage of the features of an embedded Runge-Kutta scheme designed for delivering a local error estimate at no extra-cost compared to the standard 4th order Runge-Kutta scheme.

Running time:
Highly dependent on the fibre length and accuracy required for the results. Typically between half a minute and several dozens of minutes.

References:
[1] M. Frigo and S.G. Johnson. The design and implementation of FFTW3. P. IEEE, 2(93):216-231, 2005.
[2] P.K. Janert. Gnuplot in Action, Understanding Data with Graphs. Manning Publications (2009).