Example plots for Frank Elastic Constants¶

[1]:

from FrankElastic.plot_frank import *

[2]:

%matplotlib inline
import os.path as path

Plot correlations¶

[3]:

if not path.exists("FrankElastic/Sdata.p"):
    get_data_for_Splot("FrankElastic/Sdata.p", nthreads=25, N=0.05, npts=25)

[4]:

plot_S("FrankElastic/Sdata.p")

Plots of Frank elastic constants¶

[5]:

if not path.exists("FrankElastic/highN.p"):
    get_Multi_frank_data(saveas="highN.p", gammaN=False, nthreads=25, npts=150, maxgamma=100)

[6]:

if not path.exists("FrankElastic/lowN.p"):
    get_Multi_frank_data(saveas="lowN.p", gammaN=True, nthreads=25, npts=48, maxgamma=35)

[7]:

Multi_frank_plot(load="FrankElastic/highN.p")

[8]:

Multi_frank_plot(load="FrankElastic/lowN.p")

/home/users/qmac/combinedCodeBase/wlcsim/wlcsim/FrankElastic/rigidrod.py:119: RuntimeWarning: invalid value encountered in double_scalars
  a_prime[ii] = 15*gammas[ii]*eM[0][0,0]/(4*np.sqrt(np.pi)*eM[0][0,2])

[9]:

seperate_bend_twist_splay("FrankElastic/highN.p")

[10]:

if not path.exists("FrankElastic/Nsweep.p"):
    get_Nsweep_data(saveas="FrankElastic/Nsweep.p", nthreads=25)

[11]:

plotNsweep("FrankElastic/Nsweep.p")