Exemplo n.º 1
0
def B(H,R):
#    res = cs.SX.sym('res',3)
#    res = np.array([0,0,0])
    arr = modE.calcB(R,H)
    
    return arr
Exemplo n.º 2
0
#a = modC.calcB_cyl(np.array([0.04, 0.06]), -(np.pi/2))
#print [round(a[0],5), round(a[1],5)]

''' verify that everything works as it should... '''
r = np.arange(0.05,0.1,0.0005)
a = np.arange(0, np.pi/2, 0.01)
b_cyl = np.zeros((len(a),2))
#b_bar = np.zeros((len(r),))
b_dip = np.zeros((len(r),3))
cnt = 0
startT = time.time()
for i in r:
#    b_bar[cnt] = modE.barMagnet(i)
#    b_cyl[cnt] = modC.calcB_cyl(np.array([i, 0]), np.pi/2)
    b_dip[cnt] = modE.calcB(np.array([i, 0, 0]), np.array([1, 0, 0]))    
    cnt += 1
print "time needed: ", time.time()-startT    
    
cnt = 0
for i in a:
    print "here! ", cnt
    b_cyl[cnt] = modC.calcB_cyl(np.array([r[0], 0]), i)        
    cnt += 1
    
#a = modC.calcB_cyl(np.array([r[0], 0]), np.pi/2)
    
plt.close('all')    
#plt.figure()
#plt.plot(r,b_bar,'r')    
#plt.title('bar Magnet')