Пример #1
0
    SYS.run()
    stop = time.clock()
    NO_OCL_DURATION = (stop - start) / 1000.0
    NO_OCL_OUT = SYS.fields["fibre"]

    sys = System(DOMAIN)
    sys.add(Gaussian("gaussian", peak_power=1.0, width=1.0))
    sys.add(OpenclFibre(TS, dorf="float", length=LENGTH, total_steps=STEPS))

    start = time.clock()
    sys.run()
    stop = time.clock()
    OCL_DURATION = (stop - start) / 1000.0
    OCL_OUT = sys.fields["ocl_fibre"]

    NO_OCL_POWER = temporal_power(NO_OCL_OUT)
    OCL_POWER = temporal_power(OCL_OUT)
    DELTA_POWER = NO_OCL_POWER - OCL_POWER

    MEAN_RELATIVE_ERROR = np.mean(np.abs(DELTA_POWER))
    MEAN_RELATIVE_ERROR /= np.amax(temporal_power(NO_OCL_OUT))

    print("Run time without OpenCL: %e" % NO_OCL_DURATION)
    print("Run time with OpenCL: %e" % OCL_DURATION)
    print("Mean relative error: %e" % MEAN_RELATIVE_ERROR)

    # Expect both plots to appear identical:
    double_plot(SYS.domain.t, NO_OCL_POWER, SYS.domain.t, OCL_POWER,
                x_label=labels["t"], y_label=labels["P_t"],
                X_label=labels["t"], Y_label=labels["P_t"])
Пример #2
0
system = System(Domain(bit_width=200.0, samples_per_bit=4096, channels=2))
system.add(Gaussian(width=1.0, peak_power=1.0, channel=0))
system.add(Gaussian(width=1.0, peak_power=0.5, channel=1))
system.add(
    Fibre('fibre',
          length=40.0,
          gamma=[1.0, 1.2],
          beta=[[0.0, 0.0, 0.0, 0.0], [0.0, 0.125, 0.0, 0.0]],
          sim_type='wdm',
          total_steps=400,
          method='RK4IP'))
system.run()

A_fs = system.fields['fibre']

P_nu0 = spectral_power(A_fs[0], True)
P_nu1 = spectral_power(A_fs[1], True)

double_plot(system.domain.nu,
            P_nu0,
            system.domain.nu,
            P_nu1,
            labels["nu"],
            labels["P_nu"],
            labels["nu"],
            labels["P_nu"],
            x_range=(181.1, 204.1),
            X_range=(181.1, 204.1),
            filename="7-2")
Пример #3
0
        Output information on Gaussian.
        """
        output_string = [
            'position = {0:f}', 'width = {1:f} ps', 'fwhm = {2:f} ps',
            'peak_power = {3:f} W', 'offset_nu = {4:f} THz', 'm = {5:d}',
            'C = {6:f}', 'initial_phase = {7:f} rad', 'channel = {8:d}'
        ]

        return "\n".join(output_string).format(self.position, self.width,
                                               self.calculate_fwhm(),
                                               self.peak_power, self.offset_nu,
                                               self.m, self.C,
                                               self.initial_phase,
                                               self.channel)


if __name__ == "__main__":
    """ Plot a default Gaussian in temporal and spectral domain """
    from pyofss import Domain, System, Gaussian
    from pyofss import temporal_power, spectral_power
    from pyofss import double_plot, labels

    sys = System(Domain(bit_width=500.0))
    sys.add(Gaussian())
    sys.run()

    double_plot(sys.domain.t, temporal_power(sys.field), sys.domain.nu,
                spectral_power(sys.field, True), labels["t"], labels["P_t"],
                labels["nu"], labels["P_nu"])
Пример #4
0
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian
from pyofss import phase, chirp, double_plot, labels

system = System(Domain(bit_width=10.0))
t = system.domain.t
nu = system.domain.nu
window_nu = system.domain.window_nu

system.add(Gaussian(initial_phase=3.0, width=1.0))
system.run()

double_plot(t, phase(system.field), t, chirp(system.field, window_nu),
            labels["t"], labels["phi"], labels["t"], labels["chirp"],
            filename="1 - phase_offset")

system.clear(True)
system.add(Gaussian(offset_nu=0.5, width=1.0))
system.run()

double_plot(t, phase(system.field), t, chirp(system.field, window_nu),
            labels["t"], labels["phi"], labels["t"], labels["chirp"],
            filename="2 - frequency_offset")

system.clear(True)
system.add(Gaussian(width=1.0, C=0.5))
system.run()

double_plot(t, phase(system.field), t, chirp(system.field, window_nu),
Пример #5
0
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=400.0, samples_per_bit=2048))
system.add(Gaussian(peak_power=1.0, width=30.0, C=-20.0))
system.add(
    Fibre(length=72.0, beta=[0.0, 0.0, 1.0, 0.0], gamma=1.0, total_steps=500))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t,
            P_t,
            system.domain.nu,
            P_nu_normalised,
            labels['t'],
            labels['P_t'],
            labels['nu'],
            labels['P_nu'],
            x_range=(50.0, 350.0),
            X_range=(192.6, 193.6),
            filename="4-13")
Пример #6
0
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=200.0, samples_per_bit=2048))
system.add(Gaussian(peak_power=1.0, width=1.0))
system.add(Fibre(length=5.0, beta=[0.0, 0.0, 0.0, 1.0],
                 gamma=1.0, total_steps=100))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t, P_t, system.domain.nu, P_nu_normalised,
            labels['t'], labels['P_t'], labels['nu'], labels['P_nu'],
            x_range=(95.0, 120.0), X_range=(192.6, 193.7), filename="4-15")
from pyofss import Domain, System, Gaussian
from pyofss import phase, chirp, double_plot, labels

system = System(Domain(bit_width=10.0))
t = system.domain.t
nu = system.domain.nu
window_nu = system.domain.window_nu

system.add(Gaussian(initial_phase=3.0, width=1.0))
system.run()

double_plot(t,
            phase(system.field),
            t,
            chirp(system.field, window_nu),
            labels["t"],
            labels["phi"],
            labels["t"],
            labels["chirp"],
            filename="1 - phase_offset")

system.clear(True)
system.add(Gaussian(offset_nu=0.5, width=1.0))
system.run()

double_plot(t,
            phase(system.field),
            t,
            chirp(system.field, window_nu),
            labels["t"],
            labels["phi"],
Пример #8
0
    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import spectral_power, double_plot, labels

nu_0 = 193.1
nu_1 = 1.2 * nu_0

offset_nu = 0.2 * 193.1

system = System(Domain(bit_width=30.0, samples_per_bit=8192, channels=2))
system.add(Gaussian(width=1.0, peak_power=100.0, channel=0))
system.add(Gaussian(width=1.0, peak_power=1.0, channel=1, offset_nu=offset_nu))
system.add(Fibre('fibre', length=0.4, gamma=[1.0, 1.2],
           beta=[[0.0, 0.0, 0.0, 0.0], [0.0, 10.0, 0.0, 0.0]],
           sim_type='wdm', method='ARK4IP'))
system.run()

A_fs = system.fields['fibre']

P_nu0 = spectral_power(A_fs[0], True)
P_nu1 = spectral_power(A_fs[1], True)

double_plot(system.domain.nu, P_nu0, system.domain.nu, P_nu1,
            labels["nu"], labels["P_nu"], labels["nu"], labels["P_nu"],
            x_range=(nu_0 - 8.0, nu_0 + 8.0), X_range=(nu_1 - 8.0, nu_1 + 8.0),
            filename="7-7")
Пример #9
0
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=200.0, samples_per_bit=2048))
system.add(Gaussian(peak_power=1.0, width=1.0))
system.add(Fibre(length=5.0, beta=[0.0, 0.0, 0.0, 1.0],
                 gamma=1.0, total_steps=100))
system.run()

print system.domain

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)
P_lambda_normalised = P_nu_normalised

t = system.domain.t
nu = system.domain.nu
Lambda = system.domain.Lambda

double_plot(nu, P_nu_normalised, Lambda, P_lambda_normalised, labels['nu'],
            labels['P_nu'], labels['Lambda'], labels['P_lambda'],
            x_range=(192.6, 193.7), X_range=(1547.72, 1556.55),
            filename="frequency_and_wavelength")
Пример #10
0
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=30.0, samples_per_bit=4096))
system.add(Gaussian(peak_power=1.0, width=1.0))
system.add(Fibre(length=0.1, beta=[0.0, 0.0, 0.0, 1.0],
                 gamma=100.0, total_steps=200))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t, P_t, system.domain.nu, P_nu_normalised,
            labels['t'], labels['P_t'], labels['nu'], labels['P_nu'],
            x_range=(10.0, 20.0), X_range=(190.0, 196.2), filename="4-16")
Пример #11
0
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=30.0, samples_per_bit=4096))
system.add(Gaussian(peak_power=1.0, width=1.0))
system.add(Fibre(length=0.1, beta=[0.0, 0.0, 0.0, 1.0],
                 gamma=100.0, total_steps=200))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t, P_t, system.domain.nu, P_nu_normalised,
            labels['t'], labels['P_t'], labels['nu'], labels['P_nu'],
            x_range=(-5.0, 5.0), X_range=(190.0, 196.2), filename="4-16")
Пример #12
0
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Sech, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=800.0, samples_per_bit=4096))
system.add(Sech(peak_power=1.0, width=30.0))
system.add(Fibre(length=72.0, beta=[0.0, 0.0, 1.0, 0.0],
                 gamma=1.0, total_steps=200))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t, P_t, system.domain.nu, P_nu_normalised,
            labels['t'], labels['P_t'], labels['nu'], labels['P_nu'],
            x_range=(-150.0, 150.0), X_range=(192.8, 193.4), filename="4-12")
Пример #13
0
    def __str__(self):
        """
        :return: Information string
        :rtype: string

        Output information on Gaussian.
        """
        output_string = [
            'position = {0:f}', 'width = {1:f} ps', 'fwhm = {2:f} ps',
            'peak_power = {3:f} W', 'offset_nu = {4:f} THz', 'm = {5:d}',
            'C = {6:f}', 'initial_phase = {7:f} rad', 'channel = {8:d}']

        return "\n".join(output_string).format(
            self.position, self.width, self.calculate_fwhm(), self.peak_power,
            self.offset_nu, self.m, self.C, self.initial_phase, self.channel)

if __name__ == "__main__":
    """ Plot a default Gaussian in temporal and spectral domain """
    from pyofss import Domain, System, Gaussian
    from pyofss import temporal_power, spectral_power
    from pyofss import double_plot, labels

    sys = System(Domain(bit_width=500.0))
    sys.add(Gaussian())
    sys.run()

    double_plot(sys.domain.t, temporal_power(sys.field),
                sys.domain.nu, spectral_power(sys.field, True),
                labels["t"], labels["P_t"], labels["nu"], labels["P_nu"])
Пример #14
0
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Gaussian, Fibre
from pyofss import spectral_power, double_plot, labels

system = System(Domain(bit_width=200.0, samples_per_bit=4096, channels=2))
system.add(Gaussian(width=1.0, peak_power=1.0, channel=0))
system.add(Gaussian(width=1.0, peak_power=0.5, channel=1))
system.add(Fibre('fibre', length=40.0, gamma=[1.0, 1.2],
           beta=[[0.0, 0.0, 0.0, 0.0], [0.0, 0.125, 0.0, 0.0]],
           sim_type='wdm', total_steps=400, method='RK4IP'))
system.run()

A_fs = system.fields['fibre']

P_nu0 = spectral_power(A_fs[0], True)
P_nu1 = spectral_power(A_fs[1], True)

double_plot(system.domain.nu, P_nu0, system.domain.nu, P_nu1,
            labels["nu"], labels["P_nu"], labels["nu"], labels["P_nu"],
            x_range=(181.1, 204.1), X_range=(181.1, 204.1), filename="7-2")
Пример #15
0
        phase = self.initial_phase
        phase -= 2.0 * pi * self.offset_nu * domain.t

        sechh = 1./np.cosh(t_normalised)
        sechh = np.where(sechh != 0, np.power(sechh, 1+1j*self.C), 0.)
        magnitude = sqrt(self.peak_power)*sechh

        if domain.channels > 1:
            self.field[self.channel] += magnitude * exp(1j * phase)
        else:
            self.field += magnitude * exp(1j * phase)

        return self.field

if __name__ == "__main__":
    """ Plot a default Diss_soliton in temporal and spectral domain """
    from pyofss import Domain, System, Diss_soliton
    from pyofss import temporal_power, spectral_power, inst_freq
    from pyofss import double_plot, labels

    sys = System(Domain(bit_width=500.0))
    sys.add(Diss_soliton())
    sys.run()

    double_plot(sys.domain.t, temporal_power(sys.field),
                sys.domain.nu, spectral_power(sys.field, True),
                labels["t"], labels["P_t"], labels["nu"], labels["P_nu"],
                inst_freq = inst_freq(sys.field, sys.domain.dt), y2_label=labels["inst_nu"])

Пример #16
0
from pyofss.domain import nu_to_omega
from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, double_plot, labels

domain = Domain(bit_width=20.0, samples_per_bit=8192, channels=2)

nu_0 = 193.1
nu_1 = 1.2 * nu_0

offset_nu = 0.2 * 193.1
offset = 2.5 / domain.bit_width

system = System(domain)
system.add(Gaussian(width=1.0, peak_power=1000.0, channel=0))
system.add(Gaussian(width=1.0, peak_power=0.1, channel=1,
                    position=0.5 - offset, offset_nu=offset_nu))
system.add(Fibre('fibre', length=0.2, gamma=[0.1, 0.12],
           beta=[[0.0, 0.0, 1.0, 0.0], [0.0, 10.0, 1.0, 0.0]],
           centre_omega=(nu_to_omega(nu_0), nu_to_omega(nu_1)),
           sim_type='wdm', method='ARK4IP'))
system.run()

A_fs = system.fields['fibre']

P_t0 = temporal_power(A_fs[0])
P_t1 = temporal_power(A_fs[1])

double_plot(system.domain.t, P_t0, system.domain.t, P_t1,
            labels["t"], labels["P_t"], labels["t"], labels["P_t"],
            x_range=(6.0, 14.0), X_range=(5.0, 10.0), filename="7-8")
Пример #17
0
from pyofss import Domain, System, Gaussian, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=200.0, samples_per_bit=2048))
system.add(Gaussian(peak_power=1.0, width=1.0))
system.add(
    Fibre(length=5.0, beta=[0.0, 0.0, 0.0, 1.0], gamma=1.0, total_steps=100))
system.run()

print system.domain

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)
P_lambda_normalised = P_nu_normalised

t = system.domain.t
nu = system.domain.nu
Lambda = system.domain.Lambda

double_plot(nu,
            P_nu_normalised,
            Lambda,
            P_lambda_normalised,
            labels['nu'],
            labels['P_nu'],
            labels['Lambda'],
            labels['P_lambda'],
            x_range=(192.6, 193.7),
            X_range=(1547.72, 1556.55),
            filename="frequency_and_wavelength")
Пример #18
0
          beta=[[0.0, 0.0, 1.0, 0.0], [0.0, 10.0, 1.0, 0.0]],
          centre_omega=(nu_to_omega(nu_0), nu_to_omega(nu_1)),
          sim_type='wdm',
          method='ARK4IP'))
system.run()

A_fs = system.fields['fibre']

P_t0 = temporal_power(A_fs[0])
P_t1 = temporal_power(A_fs[1])

double_plot(system.domain.t,
            P_t0,
            system.domain.t,
            P_t1,
            labels["t"],
            labels["P_t"],
            labels["t"],
            labels["P_t"],
            filename="7-6a")

P_nu0 = spectral_power(A_fs[0], True)
P_nu1 = spectral_power(A_fs[1], True)

double_plot(system.domain.nu,
            P_nu0,
            system.domain.nu,
            P_nu1,
            labels["nu"],
            labels["P_nu"],
            labels["nu"],
Пример #19
0
from pyofss import temporal_power, spectral_power, double_plot, labels

domain = Domain(bit_width=4.0, samples_per_bit=4096)

s = 0.01
width = 1.0 / (s * domain.centre_omega)
gamma = 100.0 / (width ** 2)

P_ts = []
P_nus = []
length = [20.0 / gamma, 40.0 / gamma]

for l in length:
    system = System(domain)
    system.add(Gaussian(peak_power=1.0, width=width))
    system.add(Fibre(length=l, gamma=gamma, total_steps=200,
                     self_steepening=True, beta=[0.0, 0.0, 1.0]))
    system.run()

    field = system.fields['fibre']
    P_ts.append(temporal_power(field))
    P_nus.append(spectral_power(field, True))

double_plot(system.domain.t, P_ts[0], system.domain.nu, P_nus[0],
            labels["t"], labels["P_t"], labels["nu"], labels["P_nu"],
            x_range=(-0.5, 0.5), X_range=(146.1, 240.1), filename="4-21a")

double_plot(system.domain.t, P_ts[1], system.domain.nu, P_nus[1],
            labels["t"], labels["P_t"], labels["nu"], labels["P_nu"],
            x_range=(-1.0, 1.0), X_range=(146.1, 240.1), filename="4-21b")
Пример #20
0
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
"""

from pyofss import Domain, System, Sech, Fibre
from pyofss import temporal_power, spectral_power, double_plot, labels

system = System(Domain(bit_width=800.0, samples_per_bit=4096))
system.add(Sech(peak_power=1.0, width=30.0))
system.add(Fibre(length=72.0, beta=[0.0, 0.0, 1.0, 0.0],
                 gamma=1.0, total_steps=200))
system.run()

P_t = temporal_power(system.fields['fibre'])
P_nu_normalised = spectral_power(system.fields['fibre'], True)

double_plot(system.domain.t, P_t, system.domain.nu, P_nu_normalised,
            labels['t'], labels['P_t'], labels['nu'], labels['P_nu'],
            x_range=(250.0, 550.0), X_range=(192.8, 193.4), filename="4-12")