Beispiel #1
0
def testline():
    from CoilViewer import seglistviewer
    from RingArray import fieldplot
    Rring=1.0
    seglist=[[2,(-5.0,0.0,0.0),(5.0,0.0,0.0)]]
    seglistviewer(seglist)
    filename="/tmp/tmp.txt"



    fieldfind(Rring,seglist,filename)
    fieldplot(filename,XYZS=[False,True,True,True])
Beispiel #2
0
def commandlineapp():
    import sys
    from RingArray import fieldplot

    coilfile=sys.argv[1]
    fieldfile=sys.argv[2]

    print coilfile
    if coilfile.split(".")[-1] == "xml":
        c=smoothcoil().XML_load(coilfile,globals())
        print "opened using xml"
    elif coilfile.split(".")[-1] == "cPickle":
        c=smoothcoil().load(coilfile)
        print "opened using cPickle"
      
    Rring=c.Rformer
    seglist=c.seglist
    seglistviewer(seglist)
    fieldfind(Rring,seglist,fieldfile)
    fieldplot(fieldfile,XYZS=[True,False,True,True])
Beispiel #3
0
def testring():
    from sensecoiltypes import HelixSeg
    from CoilViewer import seglistviewer
    from RingArray import fieldplot
    numsegs=30
    Rring=1.0

    dphi=2*pi/float(numsegs)
    seglist=[]
    phistart=0
    phiend=dphi
    for i in range(numsegs):
        seglist.append(HelixSeg(2,phistart,-Rring,phiend,-Rring,Rring).xyzSeg())
        seglist.append(HelixSeg(3,phiend,Rring,phistart,Rring,Rring).xyzSeg())
        phistart=phistart+dphi
        phiend=phiend+dphi

    seglistviewer(seglist)
    filename="/tmp/tmp.txt"
    fieldfind(Rring,seglist,filename)
    fieldplot(filename,XYZS=[False,True,True,True])