def relaxedScan3D(xTips, yTips, zTips):
    ntips = len(zTips)
    print(" zTips : ", zTips)
    rTips = np.zeros((ntips, 3))
    rs = np.zeros((ntips, 3))
    fs = np.zeros((ntips, 3))
    rTips[:, 0] = 1.0
    rTips[:, 1] = 1.0
    rTips[:, 2] = zTips
    fzs = np.zeros((len(zTips), len(yTips), len(xTips)))
    for ix, x in enumerate(xTips):
        print("relax ix:", ix)
        rTips[:, 0] = x
        for iy, y in enumerate(yTips):
            rTips[:, 1] = y
            itrav = PP.relaxTipStroke(rTips, rs, fs) / float(len(zTips))
            fzs[:, iy, ix] = fs[:, 2].copy()
    return fzs
def relaxedScan3D( xTips, yTips, zTips ):
	ntips = len(zTips); 
	print " zTips : ",zTips
	rTips = np.zeros((ntips,3))
	rs    = np.zeros((ntips,3))
	fs    = np.zeros((ntips,3))
	rTips[:,0] = 1.0
	rTips[:,1] = 1.0
	rTips[:,2] = zTips 
	fzs    = np.zeros(( len(zTips), len(yTips ), len(xTips ) ));
	for ix,x in enumerate( xTips  ):
		print "relax ix:", ix
		rTips[:,0] = x
		for iy,y in enumerate( yTips  ):
			rTips[:,1] = y
			itrav = PP.relaxTipStroke( rTips, rs, fs ) / float( len(zTips) )
			fzs[:,iy,ix] = fs[:,2].copy()
	return fzs
	del FFel
'''
plt.figure(figsize=( 5*nslice,5 )); plt.title( ' FF total ' )
for i in range(nslice):
	plt.subplot( 1, nslice, i+1 )
	plt.imshow( FF[i,:,:,2], origin='image', interpolation='nearest' )
'''

print " # ============  Relaxed Scan 3D "

for ix,x in enumerate( xTips  ):
	print "relax ix:", ix
	rTips[:,0] = x
	for iy,y in enumerate( yTips  ):
		rTips[:,1] = y
		itrav = PP.relaxTipStroke( rTips, rs, fs ) / float( len(zTips) )
		fzs[:,iy,ix] = fs[:,2].copy()
		#print itrav
		#if itrav > 100:
		#	print " bad convergence > %i iterations per pixel " % itrav
		#	print " exiting "
		#	break
		

print " # ============  convert Fz -> df "

dfs = PP.Fz2df( fzs, dz = dz, k0 = PP.params['kCantilever'], f0=PP.params['f0Cantilever'], n=int(PP.params['Amplitude']/dz) )

print " # ============  Plot Relaxed Scan 3D "

#slices = range( PP.params['plotSliceFrom'], PP.params['plotSliceTo'], PP.params['plotSliceBy'] )
Exemplo n.º 4
0
    del FFel
'''
plt.figure(figsize=( 5*nslice,5 )); plt.title( ' FF total ' )
for i in range(nslice):
	plt.subplot( 1, nslice, i+1 )
	plt.imshow( FF[i,:,:,2], origin='image', interpolation='nearest' )
'''

print " # ============  Relaxed Scan 3D "

for ix, x in enumerate(xTips):
    print "relax ix:", ix
    rTips[:, 0] = x
    for iy, y in enumerate(yTips):
        rTips[:, 1] = y
        itrav = PP.relaxTipStroke(rTips, rs, fs) / float(len(zTips))
        fzs[:, iy, ix] = fs[:, 2].copy()
        #print itrav
        #if itrav > 100:
        #	print " bad convergence > %i iterations per pixel " % itrav
        #	print " exiting "
        #	break

print " # ============  convert Fz -> df "

dfs = PP.Fz2df(fzs,
               dz=dz,
               k0=PP.params['kCantilever'],
               f0=PP.params['f0Cantilever'],
               n=int(PP.params['Amplitude'] / dz))