def Helmholtz(): current = 1.0 # amps Rring = 1.0 RingList = [[current, Rring / 2.0, Rring], [current, -Rring / 2.0, Rring]] fieldfile = "Helmfield.txt" fieldfind(Rring, RingList, fieldfile) fieldplot(fieldfile)
def SimpleSeg(): from CollectionViewer import collviewer_natural current=1.0 L=1.0 LineSegList=[LineSeg(Vector((0,0,0),type="CARTESIAN"),Vector((0,0,L),type="CARTESIAN"))] #collviewer_natural(linelist=[[l,(0,0,0)] for l in LineSegList]) fieldfile="SingSegfield.txt" fieldfind(L,LineSegList,fieldfile) fieldplot(fieldfile)
def Maxwell(): current = 1.0 Rring = 1.0 # [turns,z,radius] cylind = Cylinder(Vector((0, 0, -10), type="CARTESIAN"), Vector((0, 0, 10), type="CARTESIAN"), Rring) HelixSegList = RingMaker(sqrt(3) / 2.0, 1, 16, cylind) HelixSegList.extend(RingMaker(-sqrt(3) / 2.0, -1, 16, cylind)) fieldfile = "Maxwellfield.txt" fieldfind(Rring, HelixSegList, fieldfile) fieldplot(fieldfile)
def MultiSimpleSeg(): from CollectionViewer import collviewer_natural from scipy import arange current=1.0 numsegs=10 L=1.0 dL=L/float(numsegs) starts=arange(0,L,dL) LineSegList=[LineSeg(Vector((0,0,a),type="CARTESIAN"),Vector((0,0,a+dL),type="CARTESIAN")) for a in starts] #collviewer_natural(linelist=[[l,(0,0,0)] for l in LineSegList]) fieldfile="MultiSingSegfield.txt" fieldfind(L,LineSegList,fieldfile) fieldplot(fieldfile)