diff --git a/p-endogamousDiploid/ECAlist-format.py b/p-endogamousDiploid/ECAlist-format.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/p-endogamousDiploid/enumerateEndo.py b/p-endogamousDiploid/enumerateEndo.py
index 06d082f63f4dbc627261f8c4368f9ce477a0851a..cc1e8b1639be506d23699369f34328e7bcbdb9c7 100644
--- a/p-endogamousDiploid/enumerateEndo.py
+++ b/p-endogamousDiploid/enumerateEndo.py
@@ -2,9 +2,8 @@ import sys
NRULE=256
-
REFLECTION=[0,4,2,6,1,5,3,7]
-
+VERBOSE=False
CLIST=[]
seen= [[False]* NRULE]*NRULE
@@ -18,7 +17,8 @@ def scanSpace():
for f in range(NRULE):
fft=transECA(f)
fr,fc,frc=fft[1], fft[2], fft[3]
- print ("%d,%d %d,%d %d,%d"%(f, fr, f, fc, f, frc))
+ if VERBOSE:
+ print ("%d,%d %d,%d %d,%d"%(f, fr, f, fc, f, frc))
mark(f,fr) ; mark(f, fc) ; mark(f, frc)
def reflect(tab):
@@ -43,13 +43,21 @@ def exportList():
filenameout= "list-endogamous.txt"
f=open(filenameout,"w")
for couple in CLIST:
- f.write("%3d %3d\n"%(couple[0],couple[1]))
+ f.write("%d %d\n"%(couple[0],couple[1])) # ; or \n
#transECA(213)
-print('-'*20)
-scanSpace()
-print(CLIST)
-print('-'*20)
-print(len(CLIST))
-exportList()
\ No newline at end of file
+def run():
+ print('-'*20)
+ scanSpace()
+ print(CLIST)
+ print('-'*20)
+ print(len(CLIST))
+ exportList()
+
+def xcmd():
+ scanSpace()
+ for couple in CLIST:
+ print("\\tabEndo{%d}{%d}\\\\"%(couple))
+
+xcmd()
\ No newline at end of file
diff --git a/p-endogamousDiploid/list-endogamous.txt b/p-endogamousDiploid/list-endogamous.txt
index 8abebe931c39e4097166a1cbb52bbe97e209c3ec..2216a4d92ae5138cbed600a5d33297067a45c193 100644
--- a/p-endogamousDiploid/list-endogamous.txt
+++ b/p-endogamousDiploid/list-endogamous.txt
@@ -1,121 +1,121 @@
- 0 255
- 1 127
- 2 16
- 2 191
- 2 247
- 3 17
- 3 63
- 3 119
- 4 223
- 5 95
- 6 20
- 6 159
- 6 215
- 7 21
- 7 31
- 7 87
- 8 64
- 8 239
- 8 253
- 9 65
- 9 111
- 9 125
- 10 80
- 10 175
- 10 245
- 11 81
- 11 47
- 11 117
- 12 68
- 12 207
- 12 221
- 13 69
- 13 79
- 13 93
- 14 84
- 14 143
- 14 213
- 15 85
- 18 183
- 19 55
- 22 151
- 24 66
- 24 231
- 24 189
- 25 67
- 25 103
- 25 61
- 26 82
- 26 167
- 26 181
- 27 83
- 27 39
- 27 53
- 28 70
- 28 199
- 28 157
- 29 71
- 30 86
- 30 135
- 30 149
- 32 251
- 33 123
- 34 48
- 34 187
- 34 243
- 35 49
- 35 59
- 35 115
- 36 219
- 37 91
- 38 52
- 38 155
- 38 211
- 40 96
- 40 235
- 40 249
- 41 97
- 41 107
- 41 121
- 42 112
- 42 171
- 42 241
- 43 113
- 44 100
- 44 203
- 44 217
- 45 101
- 45 75
- 45 89
- 46 116
- 46 139
- 46 209
- 50 179
- 54 147
- 56 98
- 56 227
- 56 185
- 57 99
- 58 114
- 58 163
- 58 177
- 60 102
- 60 195
- 60 153
- 62 118
- 62 131
- 62 145
- 72 237
- 73 109
- 74 88
- 74 173
- 74 229
- 76 205
- 78 92
- 78 141
- 78 197
- 90 165
- 94 133
+0 255
+1 127
+2 16
+2 191
+2 247
+3 17
+3 63
+3 119
+4 223
+5 95
+6 20
+6 159
+6 215
+7 21
+7 31
+7 87
+8 64
+8 239
+8 253
+9 65
+9 111
+9 125
+10 80
+10 175
+10 245
+11 81
+11 47
+11 117
+12 68
+12 207
+12 221
+13 69
+13 79
+13 93
+14 84
+14 143
+14 213
+15 85
+18 183
+19 55
+22 151
+24 66
+24 231
+24 189
+25 67
+25 103
+25 61
+26 82
+26 167
+26 181
+27 83
+27 39
+27 53
+28 70
+28 199
+28 157
+29 71
+30 86
+30 135
+30 149
+32 251
+33 123
+34 48
+34 187
+34 243
+35 49
+35 59
+35 115
+36 219
+37 91
+38 52
+38 155
+38 211
+40 96
+40 235
+40 249
+41 97
+41 107
+41 121
+42 112
+42 171
+42 241
+43 113
+44 100
+44 203
+44 217
+45 101
+45 75
+45 89
+46 116
+46 139
+46 209
+50 179
+54 147
+56 98
+56 227
+56 185
+57 99
+58 114
+58 163
+58 177
+60 102
+60 195
+60 153
+62 118
+62 131
+62 145
+72 237
+73 109
+74 88
+74 173
+74 229
+76 205
+78 92
+78 141
+78 197
+90 165
+94 133
104 233
106 120
106 169
diff --git a/p-endogamousDiploid/oldydat2gfx.py b/p-endogamousDiploid/oldydat2gfx.py
new file mode 100644
index 0000000000000000000000000000000000000000..ccbba74ede403db4b56c4752ea3e5a1db8f7578d
--- /dev/null
+++ b/p-endogamousDiploid/oldydat2gfx.py
@@ -0,0 +1,68 @@
+import re
+import matplotlib.pyplot as plt
+import csv
+
+
+DIPLOID="0-255"
+FILETOREAD="datfiles/asymptotic-%s.dat"%DIPLOID
+print("starting pdf gen for file:%s ..."%FILETOREAD)
+
+def fileoutput():
+ return 'aphaireo-ampVSsig-%s.pdf'%FILETOREAD
+
+def read(file_path):
+
+ xVal = []
+ yVal = []
+ print ("reading dat file")
+ # open file in read mode
+ with open(file_path, 'r') as read_obj:
+ count=0
+ # iterate over each line as a ordered dictionary
+ for line in read_obj:
+ # row variable is a dictionary that represents a row in csv
+ xpar = count
+ count+=1
+ ypar = int("99")
+ xVal.append(xpar)
+ yVal.append(ypar)
+ return (xVal,yVal)
+
+def parsePar(expParse,line):
+ strparsed= re.search(expParse, line)
+ return int( strparsed.group(1) )
+
+
+
+#main : dataRep has already been read
+def processDiploid():
+
+ xval,yval=read(FILETOREAD)
+ #xvalR= tuple(100.0*(v / totalgraphweight) for v in xval) #100.0 is for percent
+
+ fig, ax = plt.subplots(1, figsize=(8, 6))
+
+ # limit lines (draw first)
+ #Yline= tuple([VDOWN,VUP] for v in xval)
+ #ax.plot(xvalR, Yline, color ='red', linewidth=1.5, linestyle="--")
+
+
+ #ax.plot(xvalR2, yval2, color="blue", linewidth=2.5, linestyle="-")
+ ax.plot(xval, yval, color="blue", linewidth=2.5, linestyle="-")
+
+ plt.xlabel("lambda (in%)")
+ plt.ylabel("density ")
+ plt.xlim(0, 100)
+ plt.ylim(0, 1.)
+
+ fig.show()
+
+ plt.savefig(fileoutput()) # produce PDF
+
+
+
+#### REPORT
+
+
+################# RUN ###############
+processDiploid()
diff --git a/p-endogamousDiploid/table-endo.tex b/p-endogamousDiploid/table-endo.tex
new file mode 100644
index 0000000000000000000000000000000000000000..9faa7958ca61d9dff2a85be7b7fe1eb2a5bf14d7
--- /dev/null
+++ b/p-endogamousDiploid/table-endo.tex
@@ -0,0 +1,168 @@
+\tabEndo{0}{255}\\
+\tabEndo{1}{127}\\
+\tabEndo{2}{16}\\
+\tabEndo{2}{191}\\
+\tabEndo{2}{247}\\
+\tabEndo{3}{17}\\
+\tabEndo{3}{63}\\
+\tabEndo{3}{119}\\
+\tabEndo{4}{223}\\
+\tabEndo{5}{95}\\
+\tabEndo{6}{20}\\
+\tabEndo{6}{159}\\
+\tabEndo{6}{215}\\
+\tabEndo{7}{21}\\
+\tabEndo{7}{31}\\
+\tabEndo{7}{87}\\
+\tabEndo{8}{64}\\
+\tabEndo{8}{239}\\
+\tabEndo{8}{253}\\
+\tabEndo{9}{65}\\
+\tabEndo{9}{111}\\
+\tabEndo{9}{125}\\
+\tabEndo{10}{80}\\
+\tabEndo{10}{175}\\
+\tabEndo{10}{245}\\
+\tabEndo{11}{81}\\
+\tabEndo{11}{47}\\
+\tabEndo{11}{117}\\
+\tabEndo{12}{68}\\
+\tabEndo{12}{207}\\
+\tabEndo{12}{221}\\
+\tabEndo{13}{69}\\
+\tabEndo{13}{79}\\
+\tabEndo{13}{93}\\
+\tabEndo{14}{84}\\
+\tabEndo{14}{143}\\
+\tabEndo{14}{213}\\
+\tabEndo{15}{85}\\
+\tabEndo{18}{183}\\
+\tabEndo{19}{55}\\
+\tabEndo{22}{151}\\
+\tabEndo{24}{66}\\
+\tabEndo{24}{231}\\
+\tabEndo{24}{189}\\
+\tabEndo{25}{67}\\
+\tabEndo{25}{103}\\
+\tabEndo{25}{61}\\
+\tabEndo{26}{82}\\
+\tabEndo{26}{167}\\
+\tabEndo{26}{181}\\
+\tabEndo{27}{83}\\
+\tabEndo{27}{39}\\
+\tabEndo{27}{53}\\
+\tabEndo{28}{70}\\
+\tabEndo{28}{199}\\
+\tabEndo{28}{157}\\
+\tabEndo{29}{71}\\
+\tabEndo{30}{86}\\
+\tabEndo{30}{135}\\
+\tabEndo{30}{149}\\
+\tabEndo{32}{251}\\
+\tabEndo{33}{123}\\
+\tabEndo{34}{48}\\
+\tabEndo{34}{187}\\
+\tabEndo{34}{243}\\
+\tabEndo{35}{49}\\
+\tabEndo{35}{59}\\
+\tabEndo{35}{115}\\
+\tabEndo{36}{219}\\
+\tabEndo{37}{91}\\
+\tabEndo{38}{52}\\
+\tabEndo{38}{155}\\
+\tabEndo{38}{211}\\
+\tabEndo{40}{96}\\
+\tabEndo{40}{235}\\
+\tabEndo{40}{249}\\
+\tabEndo{41}{97}\\
+\tabEndo{41}{107}\\
+\tabEndo{41}{121}\\
+\tabEndo{42}{112}\\
+\tabEndo{42}{171}\\
+\tabEndo{42}{241}\\
+\tabEndo{43}{113}\\
+\tabEndo{44}{100}\\
+\tabEndo{44}{203}\\
+\tabEndo{44}{217}\\
+\tabEndo{45}{101}\\
+\tabEndo{45}{75}\\
+\tabEndo{45}{89}\\
+\tabEndo{46}{116}\\
+\tabEndo{46}{139}\\
+\tabEndo{46}{209}\\
+\tabEndo{50}{179}\\
+\tabEndo{54}{147}\\
+\tabEndo{56}{98}\\
+\tabEndo{56}{227}\\
+\tabEndo{56}{185}\\
+\tabEndo{57}{99}\\
+\tabEndo{58}{114}\\
+\tabEndo{58}{163}\\
+\tabEndo{58}{177}\\
+\tabEndo{60}{102}\\
+\tabEndo{60}{195}\\
+\tabEndo{60}{153}\\
+\tabEndo{62}{118}\\
+\tabEndo{62}{131}\\
+\tabEndo{62}{145}\\
+\tabEndo{72}{237}\\
+\tabEndo{73}{109}\\
+\tabEndo{74}{88}\\
+\tabEndo{74}{173}\\
+\tabEndo{74}{229}\\
+\tabEndo{76}{205}\\
+\tabEndo{78}{92}\\
+\tabEndo{78}{141}\\
+\tabEndo{78}{197}\\
+\tabEndo{90}{165}\\
+\tabEndo{94}{133}\\
+\tabEndo{104}{233}\\
+\tabEndo{106}{120}\\
+\tabEndo{106}{169}\\
+\tabEndo{106}{225}\\
+\tabEndo{108}{201}\\
+\tabEndo{110}{124}\\
+\tabEndo{110}{137}\\
+\tabEndo{110}{193}\\
+\tabEndo{122}{161}\\
+\tabEndo{126}{129}\\
+\tabEndo{128}{254}\\
+\tabEndo{130}{144}\\
+\tabEndo{130}{190}\\
+\tabEndo{130}{246}\\
+\tabEndo{132}{222}\\
+\tabEndo{134}{148}\\
+\tabEndo{134}{158}\\
+\tabEndo{134}{214}\\
+\tabEndo{136}{192}\\
+\tabEndo{136}{238}\\
+\tabEndo{136}{252}\\
+\tabEndo{138}{208}\\
+\tabEndo{138}{174}\\
+\tabEndo{138}{244}\\
+\tabEndo{140}{196}\\
+\tabEndo{140}{206}\\
+\tabEndo{140}{220}\\
+\tabEndo{142}{212}\\
+\tabEndo{146}{182}\\
+\tabEndo{152}{194}\\
+\tabEndo{152}{230}\\
+\tabEndo{152}{188}\\
+\tabEndo{154}{210}\\
+\tabEndo{154}{166}\\
+\tabEndo{154}{180}\\
+\tabEndo{156}{198}\\
+\tabEndo{160}{250}\\
+\tabEndo{162}{176}\\
+\tabEndo{162}{186}\\
+\tabEndo{162}{242}\\
+\tabEndo{164}{218}\\
+\tabEndo{168}{224}\\
+\tabEndo{168}{234}\\
+\tabEndo{168}{248}\\
+\tabEndo{170}{240}\\
+\tabEndo{172}{228}\\
+\tabEndo{172}{202}\\
+\tabEndo{172}{216}\\
+\tabEndo{184}{226}\\
+\tabEndo{200}{236}\\