Makefile (6991B)
1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes 2 #JULIA=julia --procs 1 --color=yes 3 PREFIX=mu_sigmac 4 SCRIPT=../simulation.jl 5 PLOTSCRIPT=../plot.jl 6 NRUNS=10 7 ids = \ 8 $(PREFIX)_mu0.0_sigma_c0kPa.pdf \ 9 $(PREFIX)_mu0.1_sigma_c0kPa.pdf \ 10 $(PREFIX)_mu0.15_sigma_c0kPa.pdf \ 11 $(PREFIX)_mu0.2_sigma_c0kPa.pdf \ 12 $(PREFIX)_mu0.25_sigma_c0kPa.pdf \ 13 $(PREFIX)_mu0.3_sigma_c0kPa.pdf \ 14 $(PREFIX)_mu0.35_sigma_c0kPa.pdf \ 15 $(PREFIX)_mu0.4_sigma_c0kPa.pdf \ 16 $(PREFIX)_mu0.45_sigma_c0kPa.pdf \ 17 $(PREFIX)_mu0.0_sigma_c010kPa.pdf \ 18 $(PREFIX)_mu0.0_sigma_c100kPa.pdf \ 19 $(PREFIX)_mu0.0_sigma_c200kPa.pdf \ 20 $(PREFIX)_mu0.0_sigma_c300kPa.pdf \ 21 $(PREFIX)_mu0.0_sigma_c400kPa.pdf \ 22 $(PREFIX)_mu0.0_sigma_c500kPa.pdf \ 23 $(PREFIX)_mu0.0_sigma_c600kPa.pdf \ 24 $(PREFIX)_mu0.0_sigma_c700kPa.pdf \ 25 $(PREFIX)_mu0.0_sigma_c800kPa.pdf 26 27 .PHONY: all 28 all: $(ids) $(PREFIX)-survived_fraction.png $(PREFIX).png 29 30 all-plots: $(PREFIX)_fits.txt 31 32 $(PREFIX)_fits.txt: $(PLOTSCRIPT) 33 $(RM) $@ 34 for id in $(ids); do \ 35 echo $$id; \ 36 $(JULIA) $< \ 37 --nruns=$(NRUNS) \ 38 "$$id" >> $@;\ 39 done; 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)-survived_fraction.png: $(ids) 50 gm montage \ 51 $(ids:.pdf=.pdf-survived_fraction.pdf.png) \ 52 -geometry +4+4 \ 53 -trim \ 54 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \ 55 $@ 56 57 test.pdf: $(SCRIPT) 58 for seed in `seq 1 2`; do \ 59 $(JULIA) --track-allocation=all $< \ 60 --mu_s=0.0 --mu_d=0.0 \ 61 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 62 --tensile_strength=0000e3 \ 63 --seed=$$seed \ 64 "$@"; \ 65 done 66 $(JULIA) $(PLOTSCRIPT) \ 67 --nruns=2 \ 68 "$@" 69 70 $(PREFIX)_mu0.0_sigma_c0kPa.pdf: $(SCRIPT) 71 for seed in `seq 1 $(NRUNS)`; do \ 72 $(JULIA) $< \ 73 --mu_s=0.0 --mu_d=0.0 \ 74 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 75 --tensile_strength=0000e3 \ 76 --seed=$$seed \ 77 "$@"; \ 78 done 79 $(JULIA) $(PLOTSCRIPT) \ 80 --nruns=$(NRUNS) \ 81 "$@" 82 83 $(PREFIX)_mu0.1_sigma_c0kPa.pdf: $(SCRIPT) 84 for seed in `seq 1 $(NRUNS)`; do \ 85 $(JULIA) $< \ 86 --mu_s=0.1 --mu_d=0.1 \ 87 --mu_s_wall=0.1 --mu_d_wall=0.1 \ 88 --tensile_strength=0000e3 \ 89 --seed=$$seed \ 90 "$@"; \ 91 done 92 $(JULIA) $(PLOTSCRIPT) \ 93 --nruns=$(NRUNS) \ 94 "$@" 95 96 $(PREFIX)_mu0.15_sigma_c0kPa.pdf: $(SCRIPT) 97 for seed in `seq 1 $(NRUNS)`; do \ 98 $(JULIA) $< \ 99 --mu_s=0.15 --mu_d=0.15 \ 100 --mu_s_wall=0.15 --mu_d_wall=0.15 \ 101 --tensile_strength=0000e3 \ 102 --seed=$$seed \ 103 "$@"; \ 104 done 105 $(JULIA) $(PLOTSCRIPT) \ 106 --nruns=$(NRUNS) \ 107 "$@" 108 109 $(PREFIX)_mu0.2_sigma_c0kPa.pdf: $(SCRIPT) 110 for seed in `seq 1 $(NRUNS)`; do \ 111 $(JULIA) $< \ 112 --mu_s=0.2 --mu_d=0.2 \ 113 --mu_s_wall=0.2 --mu_d_wall=0.2 \ 114 --tensile_strength=0000e3 \ 115 --seed=$$seed \ 116 "$@"; \ 117 done 118 $(JULIA) $(PLOTSCRIPT) \ 119 --nruns=$(NRUNS) \ 120 "$@" 121 122 $(PREFIX)_mu0.25_sigma_c0kPa.pdf: $(SCRIPT) 123 for seed in `seq 1 $(NRUNS)`; do \ 124 $(JULIA) $< \ 125 --mu_s=0.25 --mu_d=0.25 \ 126 --mu_s_wall=0.25 --mu_d_wall=0.25 \ 127 --tensile_strength=0000e3 \ 128 --seed=$$seed \ 129 "$@"; \ 130 done 131 $(JULIA) $(PLOTSCRIPT) \ 132 --nruns=$(NRUNS) \ 133 "$@" 134 135 $(PREFIX)_mu0.3_sigma_c0kPa.pdf: $(SCRIPT) 136 for seed in `seq 1 $(NRUNS)`; do \ 137 $(JULIA) $< \ 138 --mu_s=0.3 --mu_d=0.3 \ 139 --mu_s_wall=0.3 --mu_d_wall=0.3 \ 140 --tensile_strength=0000e3 \ 141 --seed=$$seed \ 142 "$@"; \ 143 done 144 $(JULIA) $(PLOTSCRIPT) \ 145 --nruns=$(NRUNS) \ 146 "$@" 147 148 $(PREFIX)_mu0.35_sigma_c0kPa.pdf: $(SCRIPT) 149 for seed in `seq 1 $(NRUNS)`; do \ 150 $(JULIA) $< \ 151 --mu_s=0.35 --mu_d=0.35 \ 152 --mu_s_wall=0.35 --mu_d_wall=0.35 \ 153 --tensile_strength=0000e3 \ 154 --seed=$$seed \ 155 "$@"; \ 156 done 157 $(JULIA) $(PLOTSCRIPT) \ 158 --nruns=$(NRUNS) \ 159 "$@" 160 161 $(PREFIX)_mu0.4_sigma_c0kPa.pdf: $(SCRIPT) 162 for seed in `seq 1 $(NRUNS)`; do \ 163 $(JULIA) $< \ 164 --mu_s=0.4 --mu_d=0.4 \ 165 --mu_s_wall=0.4 --mu_d_wall=0.4 \ 166 --tensile_strength=0000e3 \ 167 --seed=$$seed \ 168 "$@"; \ 169 done 170 $(JULIA) $(PLOTSCRIPT) \ 171 --nruns=$(NRUNS) \ 172 "$@" 173 174 $(PREFIX)_mu0.45_sigma_c0kPa.pdf: $(SCRIPT) 175 for seed in `seq 1 $(NRUNS)`; do \ 176 $(JULIA) $< \ 177 --mu_s=0.45 --mu_d=0.45 \ 178 --mu_s_wall=0.45 --mu_d_wall=0.45 \ 179 --tensile_strength=0000e3 \ 180 --seed=$$seed \ 181 "$@"; \ 182 done 183 $(JULIA) $(PLOTSCRIPT) \ 184 --nruns=$(NRUNS) \ 185 "$@" 186 187 $(PREFIX)_mu0.0_sigma_c010kPa.pdf: $(SCRIPT) 188 for seed in `seq 1 $(NRUNS)`; do \ 189 $(JULIA) $< \ 190 --mu_s=0.0 --mu_d=0.0 \ 191 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 192 --tensile_strength=0010e3 \ 193 --seed=$$seed \ 194 "$@"; \ 195 done 196 $(JULIA) $(PLOTSCRIPT) \ 197 --nruns=$(NRUNS) \ 198 "$@" 199 200 $(PREFIX)_mu0.0_sigma_c100kPa.pdf: $(SCRIPT) 201 for seed in `seq 1 $(NRUNS)`; do \ 202 $(JULIA) $< \ 203 --mu_s=0.0 --mu_d=0.0 \ 204 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 205 --tensile_strength=0100e3 \ 206 --seed=$$seed \ 207 "$@"; \ 208 done 209 $(JULIA) $(PLOTSCRIPT) \ 210 --nruns=$(NRUNS) \ 211 "$@" 212 213 $(PREFIX)_mu0.0_sigma_c200kPa.pdf: $(SCRIPT) 214 for seed in `seq 1 $(NRUNS)`; do \ 215 $(JULIA) $< \ 216 --mu_s=0.0 --mu_d=0.0 \ 217 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 218 --tensile_strength=0200e3 \ 219 --seed=$$seed \ 220 "$@"; \ 221 done 222 $(JULIA) $(PLOTSCRIPT) \ 223 --nruns=$(NRUNS) \ 224 "$@" 225 226 $(PREFIX)_mu0.0_sigma_c300kPa.pdf: $(SCRIPT) 227 for seed in `seq 1 $(NRUNS)`; do \ 228 $(JULIA) $< \ 229 --mu_s=0.0 --mu_d=0.0 \ 230 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 231 --tensile_strength=0300e3 \ 232 --seed=$$seed \ 233 "$@"; \ 234 done 235 $(JULIA) $(PLOTSCRIPT) \ 236 --nruns=$(NRUNS) \ 237 "$@" 238 239 $(PREFIX)_mu0.0_sigma_c400kPa.pdf: $(SCRIPT) 240 for seed in `seq 1 $(NRUNS)`; do \ 241 $(JULIA) $< \ 242 --mu_s=0.0 --mu_d=0.0 \ 243 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 244 --tensile_strength=0400e3 \ 245 --seed=$$seed \ 246 "$@"; \ 247 done 248 $(JULIA) $(PLOTSCRIPT) \ 249 --nruns=$(NRUNS) \ 250 "$@" 251 252 $(PREFIX)_mu0.0_sigma_c500kPa.pdf: $(SCRIPT) 253 for seed in `seq 1 $(NRUNS)`; do \ 254 $(JULIA) $< \ 255 --mu_s=0.0 --mu_d=0.0 \ 256 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 257 --tensile_strength=0500e3 \ 258 --seed=$$seed \ 259 "$@"; \ 260 done 261 $(JULIA) $(PLOTSCRIPT) \ 262 --nruns=$(NRUNS) \ 263 "$@" 264 265 $(PREFIX)_mu0.0_sigma_c600kPa.pdf: $(SCRIPT) 266 for seed in `seq 1 $(NRUNS)`; do \ 267 $(JULIA) $< \ 268 --mu_s=0.0 --mu_d=0.0 \ 269 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 270 --tensile_strength=0600e3 \ 271 --seed=$$seed \ 272 "$@"; \ 273 done 274 $(JULIA) $(PLOTSCRIPT) \ 275 --nruns=$(NRUNS) \ 276 "$@" 277 278 $(PREFIX)_mu0.0_sigma_c700kPa.pdf: $(SCRIPT) 279 for seed in `seq 1 $(NRUNS)`; do \ 280 $(JULIA) $< \ 281 --mu_s=0.0 --mu_d=0.0 \ 282 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 283 --tensile_strength=0700e3 \ 284 --seed=$$seed \ 285 "$@"; \ 286 done 287 $(JULIA) $(PLOTSCRIPT) \ 288 --nruns=$(NRUNS) \ 289 "$@" 290 291 $(PREFIX)_mu0.0_sigma_c800kPa.pdf: $(SCRIPT) 292 for seed in `seq 1 $(NRUNS)`; do \ 293 $(JULIA) $< \ 294 --mu_s=0.0 --mu_d=0.0 \ 295 --mu_s_wall=0.0 --mu_d_wall=0.0 \ 296 --tensile_strength=0800e3 \ 297 --seed=$$seed \ 298 "$@"; \ 299 done 300 $(JULIA) $(PLOTSCRIPT) \ 301 --nruns=$(NRUNS) \ 302 "$@" 303 304 .PHONY: clean 305 clean: 306 $(RM) $(PREFIX)_*.run 307 $(RM) -rf $(PREFIX)_*pdf*/ 308 $(RM) -rf test*pdf*/ 309 for f in $(PREFIX)_*.vts; do $(RM) $$f; done 310 for f in $(PREFIX)_*.vtp; do $(RM) $$f; done 311 for f in $(PREFIX)_*.vtu; do $(RM) $$f; done 312 $(RM) $(PREFIX)_*.pdf 313 $(RM) $(PREFIX)_*.txt 314 $(RM) test*.pdf 315 $(RM) test*.txt