Python ProbeParticle.relaxedScan3Dの例

コード例 #1

ファイル: relaxed_scan.py プロジェクト: ProkopHapala/ProbeParticleModel

def main():
	for iq,Q in enumerate( Qs ):
		FF = FFLJ + FFel * Q
		PP.setFF_Pointer( FF )
		for ik,K in enumerate( Ks ):
			dirname = "Q%1.2fK%1.2f" %(Q,K)
			os.makedirs( dirname )
			PP.setTip( kSpring = np.array((K,K,0.0))/-PP.eVA_Nm )
			fzs = PP.relaxedScan3D( xTips, yTips, zTips )
			PP.saveXSF( dirname+'/OutFz.xsf', headScan, lvecScan, fzs )
			for iA,Amp in enumerate( Amps ):
				AmpStr = "/Amp%2.2f" %Amp
				print "Amp= ",AmpStr
				os.makedirs( dirname+AmpStr )
				dfs = PP.Fz2df( fzs, dz = dz, k0 = PP.params['kCantilever'], f0=PP.params['f0Cantilever'], n=Amp/dz )
				PP.plotImages( dirname+AmpStr+"/df", dfs, slices = range( 0, len(dfs) ) )

コード例 #2

ファイル: relaxed_scan.py プロジェクト: jhv5a/ProbeParticleModel

def main():
    for iq, Q in enumerate(Qs):
        FF = FFLJ + FFel * Q
        PP.setFF_Pointer(FF)
        for ik, K in enumerate(Ks):
            dirname = "Q%1.2fK%1.2f" % (Q, K)
            os.makedirs(dirname)
            PP.setTip(kSpring=np.array((K, K, 0.0)) / -PP.eVA_Nm)
            fzs = PP.relaxedScan3D(xTips, yTips, zTips)
            PP.saveXSF(dirname + '/OutFz.xsf', headScan, lvecScan, fzs)
            for iA, Amp in enumerate(Amps):
                AmpStr = "/Amp%2.2f" % Amp
                print("Amp= ", AmpStr)
                os.makedirs(dirname + AmpStr)
                dfs = PP.Fz2df(fzs,
                               dz=dz,
                               k0=PP.params['kCantilever'],
                               f0=PP.params['f0Cantilever'],
                               n=Amp / dz)
                PP.plotImages(dirname + AmpStr + "/df",
                              dfs,
                              slices=list(range(0, len(dfs))))

コード例 #3

ファイル: ppr.new.py プロジェクト: jhv5a/ProbeParticleModel

xTips, yTips, zTips, lvecScan = PP.prepareScanGrids()

#Ks   = [ 0.25, 0.5, 1.0 ]
#Qs   = [ -0.2, 0.0, +0.2 ]
#Amps = [ 2.0 ]

for iq, Q in enumerate(Qs):
    FF = FFLJ + FFel * Q
    PP.setFF_Pointer(FF)
    for ik, K in enumerate(Ks):
        dirname = "Q%1.2fK%1.2f" % (Q, K)
        os.makedirs(dirname)
        PP.setTip(kSpring=np.array((K, K, 0.0)) / -PP.eVA_Nm)
        #GU.saveVecFieldXsf( 'FFtot', FF, lvec, head )
        fzs = PP.relaxedScan3D(xTips, yTips, zTips)
        GU.saveXSF(dirname + '/OutFz.xsf', fzs, lvecScan, GU.XSF_HEAD_DEFAULT)
        for iA, Amp in enumerate(Amps):
            AmpStr = "/Amp%2.2f" % Amp
            print("Amp= ", AmpStr)
            os.makedirs(dirname + AmpStr)
            dz = PP.params['scanStep'][2]
            dfs = PP.Fz2df(fzs,
                           dz=dz,
                           k0=PP.params['kCantilever'],
                           f0=PP.params['f0Cantilever'],
                           n=Amp / dz)
            extent = (xTips[0], xTips[-1], yTips[0], yTips[-1])
            PPPlot.plotImages(dirname + AmpStr + "/df",
                              dfs,
                              slices=list(range(0, len(dfs))),

コード例 #4

ファイル: relaxed_scan.py プロジェクト: ProkopHapala/ProbeParticleModel

#Ks   = [ 0.25, 0.5, 1.0 ]
#Qs   = [ -0.2, 0.0, +0.2 ]
#Amps = [ 2.0 ]



for iq,Q in enumerate( Qs ):
	FF = FFLJ + FFel * Q
	PP.setFF_Pointer( FF )
	for ik,K in enumerate( Ks ):
		dirname = "Q%1.2fK%1.2f" %(Q,K)
		os.makedirs( dirname )
		PP.setTip( kSpring = np.array((K,K,0.0))/-PP.eVA_Nm )
		#GU.saveVecFieldXsf( 'FFtot', FF, lvec, head )
		fzs = PP.relaxedScan3D( xTips, yTips, zTips )
		GU.saveXSF( dirname+'/OutFz.xsf', fzs, lvecScan, GU.XSF_HEAD_DEFAULT )
		for iA,Amp in enumerate( Amps ):
			AmpStr = "/Amp%2.2f" %Amp
			print "Amp= ",AmpStr
			os.makedirs( dirname+AmpStr )
			dz  = PP.params['scanStep'][2]
			dfs = PP.Fz2df( fzs, dz = dz, k0 = PP.params['kCantilever'], f0=PP.params['f0Cantilever'], n=Amp/dz )
			extent=( xTips[0], xTips[-1], yTips[0], yTips[-1] )
			PPPlot.plotImages( dirname+AmpStr+"/df", dfs, slices = range( 0, len(dfs) ), extent=extent )

print " ***** ALL DONE ***** "

#plt.show()