#!/usr/bin/python2 # -*- coding:utf-8 -*- from __future__ import division import numpy as np # ---- this section is usefull everytime ---- import os, sys sys.path.append('../libs') from evaluation import Plotter filename = os.path.splitext(os.path.basename(__file__))[0] p = Plotter(show=True, save=True, pgf=False, name=filename, directory='') # ------------------------------------------- # ==== DEFINITIONS ==== # ==== FUNCTIONS ==== # ==== CALCULATION ==== # ==== PLOTTING ==== p.new(title=u'Titlename') p.plot(np.linspace(0, 10, 100), np.linspace(0, 10, 100)**2, '+-', label=u'Label % *.2f' % (6, 0.42)) p.make()
#!/usr/bin/python2 # -*- coding:utf-8 -*- from __future__ import division import numpy as np # ---- this section is usefull everytime ---- import os, sys sys.path.append('../libs') from evaluation import Plotter filename = os.path.splitext(os.path.basename(__file__))[0] p = Plotter(show = True, save = True, pgf = False, name=filename, directory = '') # ------------------------------------------- # ==== DEFINITIONS ==== # ==== FUNCTIONS ==== # ==== CALCULATION ==== # ==== PLOTTING ==== p.new(title=u'Titlename') p.plot(np.linspace(0,10,100), np.linspace(0,10,100)**2, '+-', label=u'Label % *.2f' % (6,0.42)) p.make()
#!/usr/bin/python # -*- coding: utf-8 -*- from __future__ import division import numpy as np # ---- this section is usefull everytime ---- import os, sys sys.path.append("../libs") from evaluation import Plotter filename = os.path.splitext(os.path.basename(__file__))[0] p = Plotter(show=True, save=False, pgf=False, name=filename, directory="") # ------------------------------------------- name, x_init, v_init, m, g = np.load("solar_system.npy") # print name # print x_init # print v_init # print m # print g def compute_forces(x, m): global g _l, c = x.shape F = np.zeros((c, c, 2)) for i in range(c): for j in range(c):
#!/usr/bin/python # -*- coding: utf-8 -*- from __future__ import division import numpy as np # ---- this section is usefull everytime ---- import os, sys sys.path.append('../libs') from evaluation import Plotter filename = os.path.splitext(os.path.basename(__file__))[0] p = Plotter(show=True, save=False, pgf=False, name=filename, directory='') # ------------------------------------------- name, x_init, v_init, m, g = np.load('solar_system.npy') #print name #print x_init #print v_init #print m #print g def compute_forces(x, m): global g _l, c = x.shape F = np.zeros((c, c, 2)) for i in range(c): for j in range(c): if j > i: F[i, j, :] = g * m[i] * m[j] * (
#!/usr/bin/python2 # -*- coding:utf-8 -*- from __future__ import division import numpy as np # ---- this section is usefull everytime ---- import os, sys sys.path.append('../libs') from evaluation import Plotter filename = os.path.splitext(os.path.basename(__file__))[0] p = Plotter(show = True, save = True, pgf = False, name=filename, directory = '') # ------------------------------------------- # ==== DEFINITIONS ==== m = 2. g = 9.81 dt = 0.1 # ==== FUNCTIONS ==== def compute_forces(x,v,vw = np.array([0.,0.]), gamma = 0.): global m, g return np.array([0.,-m*g])-gamma*(v-vw) def step_euler(x,v,f): global m, dt v = v + f/m *dt x = x + v*dt return x, v