|
| 1 | + |
| 2 | +#!/usr/bin/env bash |
| 3 | +# |
| 4 | +# Copyright 2022 Luca Fedeli |
| 5 | +# |
| 6 | +# License: BSD-3-Clause-LBNL |
| 7 | + |
| 8 | +set -eu -o pipefail |
| 9 | + |
| 10 | +export OMP_NUM_THREADS=2 |
| 11 | + |
| 12 | +# |
| 13 | +# Generate QED lookup tables using external tool |
| 14 | +# |
| 15 | +./build/bin/qed_table_generator \ |
| 16 | + --table BW --mode DP --dndt_chi_min 0.01 \ |
| 17 | + --dndt_chi_max 100 --dndt_how_many 64 \ |
| 18 | + --pair_chi_min 0.01 --pair_chi_max 100 \ |
| 19 | + --pair_chi_how_many 64 --pair_frac_how_many 64 \ |
| 20 | + -o bw_table_tool |
| 21 | +./build/bin/qed_table_generator \ |
| 22 | + --table QS --mode DP --dndt_chi_min 0.001 \ |
| 23 | + --dndt_chi_max 100 --dndt_how_many 64 \ |
| 24 | + --em_chi_min 0.001 --em_chi_max 100 \ |
| 25 | + --em_chi_how_many 64 --em_frac_how_many 64 \ |
| 26 | + --em_frac_min 1e-12 -o qs_table_tool |
| 27 | + |
| 28 | +# |
| 29 | +# Generate QED lookup tables using WarpX |
| 30 | +# |
| 31 | +./build/bin/warpx.2d \ |
| 32 | + ./Examples/Tests/qed/quantum_synchrotron/inputs_2d \ |
| 33 | + qed_bw.lookup_table_mode = "generate" \ |
| 34 | + qed_bw.tab_dndt_chi_min = 0.01 \ |
| 35 | + qed_bw.tab_dndt_chi_max = 100.0 \ |
| 36 | + qed_bw.tab_dndt_how_many = 64 \ |
| 37 | + qed_bw.tab_pair_chi_min = 0.01 \ |
| 38 | + qed_bw.tab_pair_chi_max = 100.0 \ |
| 39 | + qed_bw.tab_pair_chi_how_many = 64 \ |
| 40 | + qed_bw.tab_pair_frac_how_many = 64 \ |
| 41 | + qed_bw.save_table_in = "bw_table" \ |
| 42 | + qed_qs.lookup_table_mode = "generate" \ |
| 43 | + qed_qs.tab_dndt_chi_min = 0.001 \ |
| 44 | + qed_qs.tab_dndt_chi_max = 100.0 \ |
| 45 | + qed_qs.tab_dndt_how_many = 64 \ |
| 46 | + qed_qs.tab_em_chi_min = 0.001 \ |
| 47 | + qed_qs.tab_em_frac_min = 1.0e-12 \ |
| 48 | + qed_qs.tab_em_chi_max = 100.0 \ |
| 49 | + qed_qs.tab_em_chi_how_many = 64 \ |
| 50 | + qed_qs.tab_em_frac_how_many = 64 \ |
| 51 | + qed_qs.save_table_in = "qs_table" |
| 52 | + |
| 53 | +# |
| 54 | +# Convert lookup tables (generated with WarpX and with the external tool) in human-readable format |
| 55 | +# |
| 56 | +./build/bin/qed_table_reader -i qs_table --table QS --mode DP -o qs_table |
| 57 | +./build/bin/qed_table_reader -i qs_table_tool --table QS --mode DP -o qs_table_tool |
| 58 | +./build/bin/qed_table_reader -i bw_table --table BW --mode DP -o bw_table |
| 59 | +./build/bin/qed_table_reader -i bw_table_tool --table BW --mode DP -o bw_table_tool |
| 60 | + |
| 61 | +# |
| 62 | +# Compare the generated lookup tables |
| 63 | +# |
| 64 | +diff bw_table_dndt bw_table_tool_dndt |
| 65 | +diff bw_table_pair bw_table_tool_pair |
| 66 | +diff qs_table_phot_em qs_table_tool_phot_em |
| 67 | +diff qs_table_dndt qs_table_tool_dndt |
| 68 | + |
| 69 | +# |
| 70 | +# Run a WarpX simulation using the lookup tables generated by the external tool |
| 71 | +# |
| 72 | +./build/bin/warpx.2d \ |
| 73 | + ./Examples/Tests/qed/quantum_synchrotron/inputs_2d \ |
| 74 | + qed_bw.lookup_table_mode = "load" \ |
| 75 | + qed_bw.load_table_from = "bw_table_tool" \ |
| 76 | + qed_qs.lookup_table_mode = "load" \ |
| 77 | + qed_qs.load_table_from = "qs_table_tool" |
0 commit comments