def B(H,R):
# res = cs.SX.sym('res',3)
# res = np.array([0,0,0])
arr = modE.calcB(R,H)
return arr
#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')