Initial commitv1.0

4 files changed, 235 insertions, 0 deletions

diff --git a/figs/hierarchical_graphene.fig/Makefile b/figs/hierarchical_graphene.fig/Makefile
new file mode 100644
index 0000000..e5e1f5d
--- /dev/null
+++ b/figs/hierarchical_graphene.fig/Makefile
@@ -0,0 +1,32 @@
+PROJECTNAME=graphene_vector_field
+
+DATS=graphene_vector_field.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+graphene_vector_field.dat:
+	meankondo -C graphene.mk > graphene-numkondo.mk
+	./gendat > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+	rm -f graphene-numkondo.mk
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/hierarchical_graphene.fig/gendat b/figs/hierarchical_graphene.fig/gendat
new file mode 100755
index 0000000..d905203
--- /dev/null
+++ b/figs/hierarchical_graphene.fig/gendat
@@ -0,0 +1,30 @@
+#!/bin/bash
+
+# location of numkondo config
+config=$(dirname $0)/graphene-numkondo.mk
+
+# range and number of points in each direction
+minx=-2
+maxx=2
+nx=20
+
+miny=-1
+maxy=1
+ny=20
+
+# compute each point
+for i in $(seq $nx); do
+  # alpha_0
+  alpha0=$(python -c "print($minx+($maxx-($minx))/($nx-1)*($i-1))")
+  for j in $(seq $ny); do
+    # alpha_1
+    alpha1=$(python -c "print($miny+($maxy-($miny))/($ny-1)*($j-1))")
+
+    # print alphas
+    echo -n "$alpha0 $alpha1 "
+
+    # compute step and select the relevant lines
+    dat=$(numkondo -N 1 -F -I "0:$alpha0,1:$alpha1" "$config" | sed -r '/^[^01]/d;s/^[0-9]*://;s/,//' | tr '\n' ' ')
+    echo "$dat"
+  done
+done
diff --git a/figs/hierarchical_graphene.fig/graphene.mk b/figs/hierarchical_graphene.fig/graphene.mk
new file mode 100644
index 0000000..7bd1cca
--- /dev/null
+++ b/figs/hierarchical_graphene.fig/graphene.mk
@@ -0,0 +1,148 @@
+## This computes the logarithm of the flow equation
+
+#!fields
+# extrenal fields: xy0 where x=alpha, y=sigma
+x:110,120,210,220
+# internal fields: xy1 where x=alpha, y=sigma
+i:111,121,211,221
+# all are Fermions
+f:110,120,210,220,111,121,211,221
+
+&
+#!groups
+# different spins are independent
+(111,211) (121,221)
+
+&
+#!preprocessor_variables
+
+# psi's with internal fields
+psiAU=(1/2)[f110]+[f111],
+psiAD=(1/2)[f120]+[f121],
+psiBU=(1/2)[f210]+[f211],
+psiBD=(1/2)[f220]+[f221],
+
+psiAU-=(1/2)[f-110]+[f-111],
+psiAD-=(1/2)[f-120]+[f-121],
+psiBU-=(1/2)[f-210]+[f-211],
+psiBD-=(1/2)[f-220]+[f-221],
+
+# operators with internal fields
+O0=
+<<$psiAU>*<$psiBU->>+
+<<$psiBU>*<$psiAU->>+
+<<$psiAD>*<$psiBD->>+
+<<$psiBD>*<$psiAD->>,
+
+O1=
+<<$psiAU>*<$psiAU->*<$psiAD>*<$psiAD->>+
+<<$psiBU>*<$psiBU->*<$psiBD>*<$psiBD->>,
+
+O2=
+<<$psiAU>*<$psiAD->*<$psiBD>*<$psiBU->>+
+<<$psiBU>*<$psiBD->*<$psiAD>*<$psiAU->>+
+<<$psiAD>*<$psiAU->*<$psiBU>*<$psiBD->>+
+<<$psiBD>*<$psiBU->*<$psiAU>*<$psiAD->>,
+
+O3=
+<<$psiAU>*<$psiAU->*<$psiBU>*<$psiBU->>+
+<<$psiAD>*<$psiAD->*<$psiBD>*<$psiBD->>,
+
+O4=
+<<$psiAU>*<$psiBU->*<$psiAD>*<$psiBD->>+
+<<$psiBU>*<$psiAU->*<$psiBD>*<$psiAD->>,
+
+O5=
+<<$psiAU>*<$psiAU->*<$psiAD>*<$psiBU->*<$psiBU>*<$psiBD->>+
+<<$psiAD>*<$psiAD->*<$psiAU>*<$psiBD->*<$psiBD>*<$psiBU->>+
+<<$psiBU>*<$psiBU->*<$psiBD>*<$psiAU->*<$psiAU>*<$psiAD->>+
+<<$psiBD>*<$psiBD->*<$psiBU>*<$psiAD->*<$psiAD>*<$psiAU->>,
+
+O6=<<$psiAU>*<$psiAU->*<$psiAD>*<$psiAD->*<$psiBU>*<$psiBU->*<$psiBD>*<$psiBD->>,
+
+
+# psi's without internal fields
+phiAU=[f110],
+phiAD=[f120],
+phiBU=[f210],
+phiBD=[f220],
+
+phiAU-=[f-110],
+phiAD-=[f-120],
+phiBU-=[f-210],
+phiBD-=[f-220],
+
+# operators without internal fields
+E0=
+<<$phiAU>*<$phiBU->>+
+<<$phiBU>*<$phiAU->>+
+<<$phiAD>*<$phiBD->>+
+<<$phiBD>*<$phiAD->>,
+
+E1=
+<<$phiAU>*<$phiAU->*<$phiAD>*<$phiAD->>+
+<<$phiBU>*<$phiBU->*<$phiBD>*<$phiBD->>,
+
+E2=
+<<$phiAU>*<$phiAD->*<$phiBD>*<$phiBU->>+
+<<$phiBU>*<$phiBD->*<$phiAD>*<$phiAU->>+
+<<$phiAD>*<$phiAU->*<$phiBU>*<$phiBD->>+
+<<$phiBD>*<$phiBU->*<$phiAU>*<$phiAD->>,
+
+E3=
+<<$phiAU>*<$phiAU->*<$phiBU>*<$phiBU->>+
+<<$phiAD>*<$phiAD->*<$phiBD>*<$phiBD->>,
+
+E4=
+<<$phiAU>*<$phiBU->*<$phiAD>*<$phiBD->>+
+<<$phiBU>*<$phiAU->*<$phiBD>*<$phiAD->>,
+
+E5=
+<<$phiAU>*<$phiAU->*<$phiAD>*<$phiBU->*<$phiBU>*<$phiBD->>+
+<<$phiAD>*<$phiAD->*<$phiAU>*<$phiBD->*<$phiBD>*<$phiBU->>+
+<<$phiBU>*<$phiBU->*<$phiBD>*<$phiAU->*<$phiAU>*<$phiAD->>+
+<<$phiBD>*<$phiBD->*<$phiBU>*<$phiAD->*<$phiAD>*<$phiAU->>,
+
+E6=<<$phiAU>*<$phiAU->*<$phiAD>*<$phiAD->*<$phiBU>*<$phiBU->*<$phiBD>*<$phiBD->>
+
+&
+#!propagator
+111;211: 1 ,121;221: 1 ,
+211;111: 1 ,221;121: 1
+
+&
+#!input_polynomial
+<%exp<
+  <<[l0]>*<$O0>>+
+  <<[l1]>*<$O1>>+
+  <<[l2]>*<$O2>>+
+  <<[l3]>*<$O3>>+
+  <<[l4]>*<$O4>>+
+  <<[l5]>*<$O5>>+
+  <<[l6]>*<$O6>>
+>>
+
+&
+#!postprocess_flow_equation
+<<8>*<%log_1<$FLOW>>>
+
+&
+#!id_table
+0: <$E0>,
+1: <$E1>,
+2: <$E2>,
+3: <$E3>,
+4: <$E4>,
+5: <$E5>,
+6: <$E6>
+
+&
+#!labels
+
+0:"O0" ,
+1:"O1" ,
+2:"O2" ,
+3:"O3" ,
+4:"O4" ,
+5:"O5" ,
+6:"O6"
diff --git a/figs/hierarchical_graphene.fig/graphene_vector_field.gnuplot b/figs/hierarchical_graphene.fig/graphene_vector_field.gnuplot
new file mode 100644
index 0000000..9eef976
--- /dev/null
+++ b/figs/hierarchical_graphene.fig/graphene_vector_field.gnuplot
@@ -0,0 +1,25 @@
+set ylabel "$\\alpha_1$" norotate
+set xlabel "$\\alpha_0$"
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+unset key
+
+set pointsize 1
+
+# color functions
+
+rgb(r,g,b) = 65536 * int(r) + 256 * int(g) + int(b)
+color(x) = rgb(x*255, 0 , (1-x)*255)
+
+# get min/max
+set yrange [-1000:1000]
+stats "graphene_vector_field.dat" u (log10(sqrt(($3-$1)**2+($4-$2)**2)))
+
+set xrange [-2:2]
+set yrange [-1:1]
+
+rescale=15
+
+plot "graphene_vector_field.dat" u 1:2:(($3-$1)/sqrt(($3-$1)**2+($4-$2)**2)/rescale):(($4-$2)/sqrt(($3-$1)**2+($4-$2)**2)/rescale):(color((log10(sqrt(($3-$1)**2+($4-$2)**2))-STATS_min)/(STATS_max-STATS_min))) with vec filled head linecolor rgb variable