Makefile (4989B)
1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes 2 #JULIA=julia --procs 1 --color=yes 3 PREFIX=psd_width 4 SCRIPT=../simulation.jl 5 PLOTSCRIPT=../plot.jl 6 NRUNS=10 7 ids = \ 8 $(PREFIX)_mu0.3_sigma_c0kPa_wider.pdf \ 9 $(PREFIX)_mu0.3_sigma_c0kPa_wide.pdf \ 10 $(PREFIX)_mu0.3_sigma_c0kPa_normal.pdf \ 11 $(PREFIX)_mu0.3_sigma_c0kPa_narrow.pdf \ 12 $(PREFIX)_mu0.3_sigma_c0kPa_narrower.pdf \ 13 $(PREFIX)_mu0.0_sigma_c200kPa_wider.pdf \ 14 $(PREFIX)_mu0.0_sigma_c200kPa_wide.pdf \ 15 $(PREFIX)_mu0.0_sigma_c200kPa_normal.pdf \ 16 $(PREFIX)_mu0.0_sigma_c200kPa_narrow.pdf \ 17 $(PREFIX)_mu0.0_sigma_c200kPa_narrower.pdf 18 19 .PHONY: all 20 all: $(ids) $(PREFIX)-survived_fraction.png $(PREFIX).png 21 22 all-plots: $(PREFIX)_fits.txt 23 24 $(PREFIX)_fits.txt: $(PLOTSCRIPT) 25 $(RM) $@ 26 for id in $(ids); do \ 27 echo $$id; \ 28 $(JULIA) $< \ 29 --nruns=$(NRUNS) \ 30 "$$id" >> $@;\ 31 done; 32 33 $(PREFIX)-survived_fraction.png: $(ids) 34 gm montage \ 35 $(ids:.pdf=.pdf-survived_fraction.pdf.png) \ 36 -geometry +4+4 \ 37 -trim \ 38 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \ 39 $@ 40 41 $(PREFIX).png: $(ids) 42 gm montage \ 43 $(ids:.pdf=.pdf.png) \ 44 -geometry +4+4 \ 45 -trim \ 46 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \ 47 $@ 48 49 $(PREFIX)_mu0.3_sigma_c0kPa_wider.pdf: $(SCRIPT) 50 for seed in `seq 1 $(NRUNS)`; do \ 51 $(JULIA) $< \ 52 --mu_s=0.3 --mu_d=0.3 \ 53 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 54 --tensile_strength=0000e3 \ 55 --gamma_n=0. --gamma_t=0. \ 56 --r_min=2.0e2 --r_max=1.75e3 \ 57 --seed=$$seed \ 58 "$@"; \ 59 done 60 $(JULIA) $(PLOTSCRIPT) \ 61 --nruns=$(NRUNS) \ 62 "$@" 63 64 $(PREFIX)_mu0.3_sigma_c0kPa_wide.pdf: $(SCRIPT) 65 for seed in `seq 1 $(NRUNS)`; do \ 66 $(JULIA) $< \ 67 --mu_s=0.3 --mu_d=0.3 \ 68 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 69 --tensile_strength=0000e3 \ 70 --gamma_n=0. --gamma_t=0. \ 71 --r_min=4.0e2 --r_max=1.55e3 \ 72 --seed=$$seed \ 73 "$@"; \ 74 done 75 $(JULIA) $(PLOTSCRIPT) \ 76 --nruns=$(NRUNS) \ 77 "$@" 78 79 $(PREFIX)_mu0.3_sigma_c0kPa_normal.pdf: $(SCRIPT) 80 for seed in `seq 1 $(NRUNS)`; do \ 81 $(JULIA) $< \ 82 --mu_s=0.3 --mu_d=0.3 \ 83 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 84 --tensile_strength=0000e3 \ 85 --gamma_n=0. --gamma_t=0. \ 86 --r_min=6.0e2 --r_max=1.35e3 \ 87 --seed=$$seed \ 88 "$@"; \ 89 done 90 $(JULIA) $(PLOTSCRIPT) \ 91 --nruns=$(NRUNS) \ 92 "$@" 93 94 $(PREFIX)_mu0.3_sigma_c0kPa_narrow.pdf: $(SCRIPT) 95 for seed in `seq 1 $(NRUNS)`; do \ 96 $(JULIA) $< \ 97 --mu_s=0.3 --mu_d=0.3 \ 98 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 99 --tensile_strength=0000e3 \ 100 --gamma_n=0. --gamma_t=0. \ 101 --r_min=8.0e2 --r_max=1.15e3 \ 102 --seed=$$seed \ 103 "$@"; \ 104 done 105 $(JULIA) $(PLOTSCRIPT) \ 106 --nruns=$(NRUNS) \ 107 "$@" 108 109 $(PREFIX)_mu0.3_sigma_c0kPa_narrower.pdf: $(SCRIPT) 110 for seed in `seq 1 $(NRUNS)`; do \ 111 $(JULIA) $< \ 112 --mu_s=0.3 --mu_d=0.3 \ 113 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 114 --tensile_strength=0000e3 \ 115 --gamma_n=0. --gamma_t=0. \ 116 --r_min=1.0e3 --r_max=1.0e3 \ 117 --seed=$$seed \ 118 "$@"; \ 119 done 120 $(JULIA) $(PLOTSCRIPT) \ 121 --nruns=$(NRUNS) \ 122 "$@" 123 124 125 $(PREFIX)_mu0.0_sigma_c200kPa_wider.pdf: $(SCRIPT) 126 for seed in `seq 1 $(NRUNS)`; do \ 127 $(JULIA) $< \ 128 --mu_s=0.0 --mu_d=0.0 \ 129 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 130 --tensile_strength=0200e3 \ 131 --gamma_n=0. --gamma_t=0. \ 132 --r_min=2.0e2 --r_max=1.75e3 \ 133 --seed=$$seed \ 134 "$@"; \ 135 done 136 $(JULIA) $(PLOTSCRIPT) \ 137 --nruns=$(NRUNS) \ 138 "$@" 139 140 $(PREFIX)_mu0.0_sigma_c200kPa_wide.pdf: $(SCRIPT) 141 for seed in `seq 1 $(NRUNS)`; do \ 142 $(JULIA) $< \ 143 --mu_s=0.0 --mu_d=0.0 \ 144 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 145 --tensile_strength=0200e3 \ 146 --gamma_n=0. --gamma_t=0. \ 147 --r_min=4.0e2 --r_max=1.55e3 \ 148 --seed=$$seed \ 149 "$@"; \ 150 done 151 $(JULIA) $(PLOTSCRIPT) \ 152 --nruns=$(NRUNS) \ 153 "$@" 154 155 $(PREFIX)_mu0.0_sigma_c200kPa_normal.pdf: $(SCRIPT) 156 for seed in `seq 1 $(NRUNS)`; do \ 157 $(JULIA) $< \ 158 --mu_s=0.0 --mu_d=0.0 \ 159 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 160 --tensile_strength=0200e3 \ 161 --gamma_n=0. --gamma_t=0. \ 162 --r_min=6.0e2 --r_max=1.35e3 \ 163 --seed=$$seed \ 164 "$@"; \ 165 done 166 $(JULIA) $(PLOTSCRIPT) \ 167 --nruns=$(NRUNS) \ 168 "$@" 169 170 $(PREFIX)_mu0.0_sigma_c200kPa_narrow.pdf: $(SCRIPT) 171 for seed in `seq 1 $(NRUNS)`; do \ 172 $(JULIA) $< \ 173 --mu_s=0.0 --mu_d=0.0 \ 174 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 175 --tensile_strength=0200e3 \ 176 --gamma_n=0. --gamma_t=0. \ 177 --r_min=8.0e2 --r_max=1.15e3 \ 178 --seed=$$seed \ 179 "$@"; \ 180 done 181 $(JULIA) $(PLOTSCRIPT) \ 182 --nruns=$(NRUNS) \ 183 "$@" 184 185 $(PREFIX)_mu0.0_sigma_c200kPa_narrower.pdf: $(SCRIPT) 186 for seed in `seq 1 $(NRUNS)`; do \ 187 $(JULIA) $< \ 188 --mu_s=0.0 --mu_d=0.0 \ 189 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 190 --tensile_strength=0200e3 \ 191 --gamma_n=0. --gamma_t=0. \ 192 --r_min=1.0e3 --r_max=1.0e3 \ 193 --seed=$$seed \ 194 "$@"; \ 195 done 196 $(JULIA) $(PLOTSCRIPT) \ 197 --nruns=$(NRUNS) \ 198 "$@" 199 200 .PHONY: clean 201 clean: 202 $(RM) $(PREFIX)_*.run 203 $(RM) -rf $(PREFIX)_*pdf*/ 204 $(RM) -rf test*pdf*/ 205 for f in $(PREFIX)_*.vts; do $(RM) $$f; done 206 for f in $(PREFIX)_*.vtp; do $(RM) $$f; done 207 for f in $(PREFIX)_*.vtu; do $(RM) $$f; done 208 $(RM) $(PREFIX)_*.pdf 209 $(RM) $(PREFIX)_*.png 210 $(RM) $(PREFIX)_*.txt 211 $(RM) test*.pdf 212 $(RM) test*.txt