Changeset 3411


Ignore:
Timestamp:
May 29, 2018 1:36:05 PM (5 years ago)
Author:
vondreele
Message:

another transform constrint fix to atom nos.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIconstrGUI.py

    r3410 r3411  
    12601260                continue
    12611261            name = xnames[ix-1]
    1262             IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%s'%(npId,name,iat))]
     1262            IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%d'%(npId,name,ia))]
    12631263            DepCons = []
    12641264            for iop,opval in enumerate(XOpr[ix-1]):
    12651265                if opval:
    1266                     DepCons.append([opval,G2obj.G2VarObj('%d::%s:%d'%(opId,xnames[iop],ia))])
     1266                    DepCons.append([opval,G2obj.G2VarObj('%d::%s:%s'%(opId,xnames[iop],iat))])
    12671267            if len(DepCons) == 1:
    12681268                constraints['Phase'].append([IndpCon,DepCons[0],None,None,'e'])
     
    12721272                constraints['Phase'].append([IndpCon]+DepCons+[0.0,None,'c'])
    12731273        for name in ['Afrac','AUiso']:
    1274             IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%s'%(npId,name,iat))]
    1275             DepCons = [1.0,G2obj.G2VarObj('%d::%s:%d'%(opId,name,ia))]
     1274            IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%d'%(npId,name,ia))]
     1275            DepCons = [1.0,G2obj.G2VarObj('%d::%s:%s'%(opId,name,iat))]
    12761276            constraints['Phase'].append([IndpCon,DepCons,None,None,'e'])
    12771277#        for iu,Uid in enumerate(Uids):
     
    12821282            IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%s'%(npId,Uid[2],ia))]
    12831283            DepCons = []
    1284             for iat in range(3):
    1285                 for ibt in range(3):
    1286                     if abs(Trans[Uid[0],iat]) > 1.e-4 and abs(Trans[Uid[1],ibt]) > 1.e-4:
    1287                         parm = '%d::%s:%d'%(opId,unames[ibt][iat],ia)
     1284            for iau in range(3):
     1285                for ibu in range(3):
     1286                    if abs(Trans[Uid[0],iau]) > 1.e-4 and abs(Trans[Uid[1],ibu]) > 1.e-4:
     1287                        parm = '%d::%s:%s'%(opId,unames[ibu][iau],iat)
    12881288                        if not parm in varyList:
    12891289                            varyList.append(parm)
    1290                         DepCons.append([Trans[ibt,iat]/detTrans,G2obj.G2VarObj(parm)])
     1290                        DepCons.append([Trans[ibu,iau]/detTrans,G2obj.G2VarObj(parm)])
    12911291            if len(DepCons) == 1:
    12921292                constraints['Phase'].append([IndpCon,DepCons[0],None,None,'e'])
Note: See TracChangeset for help on using the changeset viewer.