nzm

p-endogamousDiploid/.gitignore

+user.properties

p-endogamousDiploid/enumerateEndo.py

@@ -2,7 +2,7 @@ import sys
 
 
 TCODE=['A','B','C','D','E','F','G','H']
-TINDX=[0  , 1 , 4,  5,  2, 3, 6, 7]
+TINDX=[0  , 1 , 4 ,  5,  2,  3,  6, 7 ]
 
 NRULE=256
 REFLECTION=[0,4,2,6,1,5,3,7]
@@ -65,8 +65,9 @@ def exportList():
 
 #transECA(213)
 
-def transtionTabTcode(tcodef, tcodeg):
+def transtionTabTcode(tcodepair):
     """ transition table of the diploid in the T-code mode"""
+    (tcodef, tcodeg)=tcodepair
     tabtcode="["
     for i in range(8):
         ch= TCODE[i]
@@ -91,11 +92,53 @@ def run():
     exportList()
 
 def xcmd():
-    scanSpace()
+    
     for couple in CLIST:
         tabf,tabg=wcodeToTransitionTable(couple[0]),wcodeToTransitionTable(couple[1])
         tcodef,tcodeg=Tcode(tabf),Tcode(tabg)
         ttab=transtionTabTcode(tcodef, tcodeg)
         print("\\tabEndo{%d}{%d}{%s}{%s}{%s}\\\\"%(couple[0], couple[1], tcodef, tcodeg, ttab))
 
-xcmd()
+def pres(tcodepair, ch):
\ No newline at end of file
+    return ch in tcodepair[0] or ch in tcodepair[1]
+
+def filterZeroOnlyFP():
+    print("looking for ZERO.FP")
+    for couple in CLIST:
+        tabf,tabg=wcodeToTransitionTable(couple[0]),wcodeToTransitionTable(couple[1])
+        tcodepair=(Tcode(tabf),Tcode(tabg))
+        strttab=transtionTabTcode(tcodepair)
+        ttab=strttab.replace("[","").replace("]","").split(",")
+        #if ttab[0]=='0':
+        condA=pres(tcodepair, 'A')
+        condH=pres(tcodepair, 'H')
+        condDE=pres(tcodepair, 'D') or pres(tcodepair, 'E')
+        condBEC=pres(tcodepair, 'B') or pres(tcodepair, 'E') or pres(tcodepair, 'C')
+        condGDF=pres(tcodepair, 'G') or pres(tcodepair, 'D') or pres(tcodepair, 'F')
+        condBFGC=pres(tcodepair, 'B') or pres(tcodepair, 'F') or pres(tcodepair, 'G') or pres(tcodepair, 'C')
+        wordfp=""
+        if not condA:
+            wordfp+="a"
+        if not condH:
+            wordfp+="b"
+        if not condDE:
+            wordfp+="c"
+        if not condBEC:
+            wordfp+="d"
+        if not condGDF:
+            wordfp+="e"
+        if not condBFGC:
+            wordfp+="f"
+        if wordfp=="":
+            wordfp="n"                
+        print("wFP:%s  %3d-%3d %s,%s   %s"%(wordfp, couple[0], couple[1], tcodepair[0], tcodepair[1], strttab))
+        #if not condA and condH and condDE and condBEC and condGDF and condBFGC:
+        #    print("ZFP:%s  %d-%d %s,%s   %s"%(ttab[1], couple[0], couple[1], tcodepair[0], tcodepair[1], strttab))
+        
+
+
+
+
+scanSpace()
+#xcmd()
+filterZeroOnlyFP()