figs/hierarchical_graphene.fig/gendat


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
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