Changeset 3961 for trunk/GSASIIplot.py


Ignore:
Timestamp:
May 8, 2019 6:16:01 PM (3 years ago)
Author:
vondreele
Message:

set self.tutorialPath = '.\G2tutorial' if others fail. Don't know if this is best?
Add S301 to trigonal to 3m1 list of generalized mustrain terms - seems to improve things but plot is weird. Might need more work.
use get to fix missing Z problem
set up getting better map density from density slice instead of 3/4D full map

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r3959 r3961  
    76757675
    76767676    # PlotStructure initialization here
    7677     global mcsaXYZ,mcsaTypes,mcsaBonds,txID,contourSet
     7677    global mcsaXYZ,mcsaTypes,mcsaBonds,txID,contourSet,Zslice
    76787678    global cell, Vol, Amat, Bmat, A4mat, B4mat
    76797679    txID = 0
     
    89048904                    RenderPlane(plane,color)
    89058905        if drawingData.get('showSlice',False):
    8906             global contourSet
     8906            global contourSet,Zslice
    89078907            if len(D4mapData.get('rho',[])):        #preferentially select 4D map if there
    89088908                modQ = np.array(generalData['SuperVec'][0])
     
    89268926            if FourD:
    89278927                SXYZT = np.vstack((SXYZ.T,np.inner(SXYZ,modQ)+G2frame.tau)).T
    8928                 Z = np.reshape(map_coordinates(rho,(SXYZT%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))
     8928                Zslice = np.reshape(map_coordinates(rho,(SXYZT%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))
    89298929            else:
    8930                 Z = np.reshape(map_coordinates(rho,(SXYZ%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))
    8931             Z = np.where(Z<=Rmax,Z,Rmax)
     8930                Zslice = np.reshape(map_coordinates(rho,(SXYZ%1.*rho.shape).T,order=1,mode='wrap'),(npts,npts))
     8931            Z = np.where(Zslice<=Rmax,Zslice,Rmax)
    89328932            plt.rcParams['figure.facecolor'] = (1.,1.,1.,.5)
    89338933            oldSize = plt.rcParams['figure.figsize']
Note: See TracChangeset for help on using the changeset viewer.