コード例 #1
0
def profile_driver():
    if 1:
        if 0:
            #NEED TO CHANGE THESE! SPECIFIC TO BILL"S MACHINE
            spinfile=r'C:/Documents and Settings/wflynn/Desktop/spins1.txt.'#'C:/eig_test_Spins.txt'
            interactionfile=r'C:/Documents and Settings/wflynn/Desktop/montecarlo1.txt'#'C:/eig_test_montecarlo.txt'
        if 0:
            spinfile=r'C:/Documents and Settings/wflynn/Desktop/spins.txt'
            interactionfile=r'C:/Documents and Settings/wflynn/Desktop/yang_montecarlo.txt'
        if 1:
            #NEED TO CHANGE THESE! SPECIFIC TO BILL"S MACHINE
            spinfile=r'C:/Documents and Settings/wflynn/Desktop/spinwave_test_spins1.txt'#'C:/eig_test_Spins.txt'
            interactionfile=r'C:/Documents and Settings/wflynn/Desktop/spinwave_test_montecarlo1.txt'
        
        steps=100
        data={}
        data['kx']=1.
        data['ky']=0.
        data['kz']=0.
        direction=data
        Hsave = driver1(spinfile,interactionfile)
        sys.exit()
        eigs = Hsave.eigenvals().keys()
        print 'eigs'
        print eigs
        printing.generate_output(eigs)
        driver2(Hsave,direction,steps,0,2*N.pi)
コード例 #2
0
            #NEED TO CHANGE THESE! SPECIFIC TO BILL"S MACHINE
            spinfile=r'C:/Documents and Settings/wflynn/Desktop/fm_chain_spins_x.txt'#'C:/eig_test_Spins.txt'
            interactionfile=r'C:/Documents and Settings/wflynn/Desktop/fm_chain_montecarlo_x.txt'
        if 0:
            spinfile=r'C:/Documents and Settings/wflynn/Desktop/spins.txt'
            interactionfile=r'C:/Documents and Settings/wflynn/Desktop/yang_montecarlo.txt'
        
        steps=100
        data={}
        data['kx']=1.
        data['ky']=0.
        data['kz']=0.
        direction=data
        Hsave = driver1(spinfile,interactionfile)
        sys.exit()
        eigs = Hsave.eigenvals().keys()
        print 'eigs'
        print eigs
        printing.generate_output(eigs)
        driver2(Hsave,direction,steps,0,2*N.pi)
        #atom_list, jnums, jmats=readfiles.read_interactions(myfilestr,spins)
        #N_atoms=len(atom_list)
        #N_atoms_uc=1
        #print N_atoms
        #atom_list=generate_atoms()
        #calculate_dispersion(atom_list,N_atoms_uc,N_atoms,jmats)
        #print jmats



    
コード例 #3
0
        #interactionfile=r'c:\montecarlo_square.txt'

        #simple cubic fm
        #spinfile=r'C:/Documents and Settings/wflynn/My Documents/workspace/spinwaves/spinwaves/spinwavecalc/tests/spins_sc.txt'
        #interactionfile=r'C:/Documents and Settings/wflynn/My Documents/workspace/spinwaves/spinwaves/spinwavecalc/tests/montecarlo_sc.txt'

        #NEED TO CHANGE THESE! SPECIFIC TO BILL"S MACHINE
        pathfilename = os.path.join('..')
        spinfile = os.path.join(pathfilename, 'cross_section/test_Spins.txt')
        interactionfile = os.path.join(pathfilename,
                                       'cross_section/test_montecarlo.txt')

        steps = 100
        data = {}
        data['kx'] = 1.
        data['ky'] = 0.
        data['kz'] = 0.
        direction = data
        Hsave = driver1(spinfile, interactionfile)
        eigs = Hsave.eigenvals().keys()
        print eigs
        printing.generate_output(eigs)
        driver2(Hsave, direction, steps, 0, 2 * N.pi)
        #atom_list, jnums, jmats=readfiles.read_interactions(myfilestr,spins)
        #N_atoms=len(atom_list)
        #N_atoms_uc=1
        #print N_atoms
        #atom_list=generate_atoms()
        #calculate_dispersion(atom_list,N_atoms_uc,N_atoms,jmats)
        #print jmats