In [8]:

    

%matplotlib notebook
import interface
import ipyvolume

smod = interface.interface(0)
smod.init_interface(0)
smod.get_modified_interface(1, 1.9, c=8e-2)
#ipyvolume.quickvolshow(smod.sn.f_field, lighting=True, level=[0.4, 0.6, 0.9])


    

    




0.06097674369812012 RHOEN INIT
0.11967134475708008 RHOIN INIT
0.17737889289855957 RHOEA INIT
0.23481011390686035 RHOIA INIT
0.3080177307128906 NEUTRAL INTERFACE INITIALIZED
0.36444997787475586 ANION INTERFACE INITIALIZED
40.20291519165039 SPHERE
40.31625247001648 RHOTOT






    




/home/mat/.conda/envs/environ/lib/python3.7/site-packages/mkl_fft/_numpy_fft.py:1044: FutureWarning: Using a non-tuple sequence for multidimensional indexing is deprecated; use `arr[tuple(seq)]` instead of `arr[seq]`. In the future this will be interpreted as an array index, `arr[np.array(seq)]`, which will result either in an error or a different result.
  output = mkl_fft.rfftn_numpy(a, s, axes)






    




(192, 192, 192)
smod written to file
65.93408346176147 SMOD

In [21]:

    

smod.get_modified_interface(1, 0, c=0.18)


    

    




283.9699902534485 SPHERE
284.08764004707336 RHOTOT






    




/home/mat/.conda/envs/environ/lib/python3.7/site-packages/mkl_fft/_numpy_fft.py:1044: FutureWarning: Using a non-tuple sequence for multidimensional indexing is deprecated; use `arr[tuple(seq)]` instead of `arr[seq]`. In the future this will be interpreted as an array index, `arr[np.array(seq)]`, which will result either in an error or a different result.
  output = mkl_fft.rfftn_numpy(a, s, axes)






    




(192, 192, 192)
smod written to file
310.33499479293823 SMOD

In [2]:

    

ipyvolume.quickvolshow(smod.sconv.f_field, lighting=True, level=[0.4, 0.6, 0.9])


    

In [6]:

    

ipyvolume.quickvolshow(smod.smodified.f_field, lighting=True, level=[0.4, 0.6, 0.9])


    

    




/home/mat/.conda/envs/environ/lib/python3.7/site-packages/ipyvolume/serialize.py:66: RuntimeWarning: invalid value encountered in true_divide
  gradient = gradient / np.sqrt(gradient[0]**2 + gradient[1]**2 + gradient[2]**2)






    






 
 

In [5]:

    

#smod.rhotot.plotnormalizedplane(2, 100)
smod.rhoconv.plotnormalizedplane(2, 100)
#smod.sconv.plotnormalizedplane(2, 100)
#smod.sa.plotnormalizedplane(2, 100)
#smod.sn.plotsimpleplane(2, 100)
#smod.sphere.plotsimpleplane(2, 100)


    

In [17]:

    

#smod.rhotot.plotnormalizedplane(2, 100)
#smod.rhoconv.plotnormalizedplane(2, 100)
smod.sconv.plotnormalizedplane(2, 100)
#smod.sa.plotnormalizedplane(2, 100)
#smod.sn.plotsimpleplane(2, 100)
#smod.sphere.plotsimpleplane(2, 100)


    

In [22]:

    

smod.sconv.plotnormalizedplane(2, 100)


    

In [11]:

    

import line
import numpy as np
smod.sn.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))

smod.smodified.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))
smod.sa.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))


    

In [12]:

    

import line
import numpy as np
smod.sn.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))

smod.smodified.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))
smod.sa.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))


    

In [2]:

    

smod.rhotot.plotsimpleplane(2, 100)


    

In [3]:

    

ipyvolume.quickvolshow(smod.rhotot.f_field, lighting=True, level=[0.4, 0.6, 0.9])


    

In [12]:

    

import utils_functions as uf
x = uf.sigmoid(uf.SigType.erf, 0, 1, 0, 0.1)
print(x)


    

    




None

In [17]:

    

a = smod.sn.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])))
print(a)


    

    




None

In [23]:

    

import matplotlib.pyplot as plt
smod.sn.plotline(line.line(pos=np.array([13, 13, 0]), dir=np.array([0, 0, 1])), 0, True)
plt.show()


    

In [ ]: