VP1D4f_dict.py

###################################################################################################
##               Dictionary of parameters: https://github.com/cchandre/Vlasov1D                  ##
###################################################################################################

import numpy as xp
from scipy.special import erf

kappa = 1.30834
Tf = 80

integrator_kinetic = 'position-Verlet'
nsteps = 40
integrator_fluid = 'DOP853'
precision = 1e-13

n_casimirs = 3

Lx = 2 * xp.pi
Lv = 10
Nx = 2**9
Nv = 2**12

A = 1e-6
k = xp.pi / Lx
# alpha = 0.1
f_init = lambda x, v: (1 - A * xp.cos(k * x)) * v**2 * xp.exp(-v**2 / 2) / xp.sqrt(2 * xp.pi)
# f_init = lambda x, v: (1 - A * xp.cos(k * x)) * xp.exp(-v**2 / 2) / xp.sqrt(2 * xp.pi) * (1 + erf(alpha * v / xp.sqrt(2)))

output_var = 'E'
output_modes = 8

## 'Compute', 'Plot' and/or 'Save'
Kinetic = ['Compute', 'Plot', 'Save']
Fluid = ['Compute', 'Plot', 'Save']

darkmode = True

###################################################################################################
##                              DO NOT EDIT BELOW                                                ##
###################################################################################################
dict = {'kappa': kappa,
        'Tf': Tf,
        'integrator_kinetic': integrator_kinetic,
        'nsteps': nsteps,
        'integrator_fluid': integrator_fluid,
        'precision': precision,
        'n_casimirs': n_casimirs,
        'Lx': Lx,
        'Lv': Lv,
        'Nx': Nx,
        'Nv': Nv,
        'f_init': f_init,
        'output_var': output_var,
        'output_modes': output_modes,
        'Kinetic': Kinetic,
        'Fluid': Fluid,
        'darkmode': darkmode}
###################################################################################################