import numpy as np
import scipy as sp
from matplotlib import pyplot as plt
import os, optparse, glob, time
from LE_Utils import filename_pars
from LE_SPressure import force_dlin, pdf_WN, plot_wall
from LE_Utils import plot_fontsizes
fsa, fsl, fst = plot_fontsizes()
def main():
"""
Plot components of the bulk constant.
Adapted from LE_BulkConst.py
For dlin only so far.
"""
me = "test_dP: "
t0 = time.time()
parser = optparse.OptionParser(conflict_handler="resolve")
parser.add_option('-s',
'--show',
dest="showfig",
default=False,
action="store_true")
parser.add_option('--nosave',
dest="nosave",
default=False,
action="store_true")
import warnings
from time import time as sysT
from itertools import chain
from sys import argv
from LE_LightBoundarySim import lookup_xmax,calculate_xmin,calculate_xini
from LE_Utils import force_1D_const, force_1D_lin
from LE_Utils import save_data, filename_pars
warnings.filterwarnings("ignore",
"No labelled objects found. Use label='...' kwarg on individual plots.",
UserWarning)
## Global variables
from LE_Utils import plot_fontsizes
fsa,fsl,fst = plot_fontsizes()
def main():
"""
NAME
LE_Pressure.py
PURPOSE
Calculate pressure in vicinity of linear potential for particles driven
by exponentially correlated noise.
EXECUTION
python LE_Pressure.py histfile/directory flags
ARGUMENTS
import numpy as np
from matplotlib import pyplot as plt
from time import time
from LE_SBS import force_dlin, force_nu, force_dnu,\
fxy_infpot
from LE_Utils import plot_fontsizes
plt.rcParams.update({"axes.labelsize": plot_fontsizes()[0]})
R = 5.0
S = 5.0
lam = 0.5
wob = R + lam
wib = S - lam
nu = 10.0
dt = 0.01
force = lambda xy, r: force_dnu(xy, r, R, S, lam, nu)
## FORCE MAGNITUDE -- 2D
fxy = lambda xy, r: force_dlin(xy, r, R, S)
N = 100
x = np.linspace(-2 * R, 2 * R, N)
y = np.linspace(-2 * R, 2 * R, N)
farr = np.array([
fxy(np.array([xi, yi]), np.sqrt(xi * xi + yi * yi)) for xi in x for yi in y
])
farrabs = np.sqrt(farr[:, 0] * farr[:, 0] + farr[:, 1] * farr[:, 1]).reshape(
[N, N])
import numpy as np
import scipy as sp
from matplotlib import pyplot as plt
from platform import system
from datetime import datetime
from LE_Utils import plot_fontsizes
## ============================================================================
"""
Statistics of simulation.
"""
outdir = "Pressure/"+str(datetime.now().strftime("%y%m%d"))+"_CIR_DN_dt"
nosave = 1
annoloc = (0.02,0.87)
plt.rcParams.update({"axes.labelsize": plot_fontsizes()[0]})
## ============================================================================
## ============================================================================
def plot_step_wall(xy,rcoord,R,S,a,dt,vb):
"""
Distributions of spatial steps in wall regions
"""
me = "LE_RunPlot.plot_step_wall: "
fig, axs = plt.subplots(2,2); axs = axs.reshape([axs.size])
##-------------------------------------------------------------------------
## distribution of x step
ax = axs[0]
xstep = np.hstack([np.diff(xy[:,0]),0])