#!/usr/bin/env python # -*- coding: utf8 -*- import sys if len(sys.argv) > 1: mode = sys.argv[1] else: mode = 'both' import elasticite as el e = el.EdgeGrid(mode=mode) el.main(e)
z = montage(z, z_s['line_fresnelastique'])
###########################################################################
z = montage(z, revert(z_s['line_onde_dense']))
# z = montage(z, revert(z_s['line_onde_solo']))
###########################################################################
# check that there is not overflow @ 30 fps
el.check(e, z)
###########################################################################
# save the file
np.save(filename, z)
return z_s
if __name__ == "__main__":
import sys
if len(sys.argv) > 1: mode = sys.argv[1]
else: mode = 'both'
filename = 'mat/master.npy'
e = el.EdgeGrid(N_lame=25,
grid_type='line',
mode=mode,
verb=False,
filename=filename)
if mode == 'writer':
z_s = master(e, filename)
else:
# running the code
el.main(e)
#!/usr/bin/env python
# -*- coding: utf8 -*-
"""
Sur une ligne de lames, on fait tourner les lames avec un mouvement émergent.
"""
import elasticite as el
EdgeGrid = el.EdgeGrid
if __name__ == "__main__":
import sys
if len(sys.argv) > 1: mode = sys.argv[1]
else: mode = 'both'
e = el.EdgeGrid(N_lame=20, grid_type='line', mode=mode, verb=False)
el.main(e)
"""
Dans une grille en tore (pacman) privilégie les co-linéarités à angle droit.
"""
import elasticite as el
import numpy as np
class EdgeGrid(el.EdgeGrid):
def champ(self):
force = np.zeros_like(self.lames[2, :])
damp = lambda t: 0.001
noise = lambda t: 2* np.exp((np.cos(2*np.pi*(t-0.) / 6.)-1.)/ 1.**2)
colin_t = lambda t: -2.*np.exp((np.cos(2*np.pi*(t-4.) / 6.)-1.)/ 1.**2)
colin_d = lambda d: np.exp(-d/.05)
#colin_d = lambda d: -1. + 2*np.exp(-np.log((d+1.e-12)/.1)**2/2/1.5**2)
delta_angle = self.angle_relatif()-np.pi/2.
#delta_angle *= np.sign(delta_angle)
force += colin_t(self.t) * np.sum(np.sin(4*delta_angle)*colin_d(self.distance(do_torus=True)), axis=1)
force += noise(self.t)*np.pi*np.random.randn(self.N_lame)
force -= damp(self.t) * self.lames[3, :]/self.dt
return 1. * force
if __name__ == "__main__":
import sys
if len(sys.argv)>1: mode = sys.argv[1]
else: mode = 'both'
e = el.EdgeGrid(N_lame=72, mode=mode, verb=False)
el.main(e)