Commit bde18b31 authored by pia.stammer's avatar pia.stammer
Browse files

Commented icru.h and io.h + added config and mesh file for 2D phantom

parent 7f80e4a6
Pipeline #119904 canceled with stages
/*!
* @file io.h
* @brief Creates required input from given config file
*/
#ifndef IO_H
#define IO_H
......@@ -9,17 +13,41 @@
class Config;
class Mesh;
/*!
* @brief
* @param[in] fileName -name of the file that results shall be exported to
* @param[in] results - results of computation
* @param[in] fieldNames - specify different output values
* @param[in] mesh - mesh results were computed on
*/
void ExportVTK( const std::string fileName,
const std::vector<std::vector<std::vector<double>>>& results,
const std::vector<std::string> fieldNames,
const Mesh* mesh );
/*!
* @brief Loads mesh from an SU2 file, returns Mesh element
*/
Mesh* LoadSU2MeshFromFile( const Config* settings );
/*!
* @brief Parses arguments given when calling program from command line
* @param[in] argc - number arguments
* @param[in] argv - string with arguments given
*/
std::string ParseArguments( int argc, char* argv[] );
/*!
* @brief prints configurations set in inputFile to console
*/
void PrintLogHeader( std::string inputFile );
/*!
* @brief creates a mesh from image given in imageName and writes it to SU2Filename
* @param[in] imageName - name of image file
* @param[in] SU2Filename - name of SU2 file where output is saved
*/
Matrix createSU2MeshFromImage( std::string imageName, std::string SU2Filename );
#endif // IO_H
/*!
* @file icru.h
* @brief Extracts physical values from ICRU 77 database
*
* Disclaimer: This is copied from the code distributed with the ICRU 77 report
* Link: https://icru.org/home/reports/elastic-scattering-of-electrons-and-positrons-icru-report-77
*/
#ifndef ICRU_H
#define ICRU_H
......
NDIME= 2
NELEM= 162
5 67 78 37 0
5 37 71 67 1
5 34 68 61 2
5 39 62 58 3
5 61 69 34 4
5 39 85 62 5
5 14 70 13 6
5 40 66 64 7
5 36 75 66 8
5 64 87 40 9
5 13 70 37 10
5 34 70 14 11
5 35 64 63 12
5 35 63 62 13
5 48 77 65 14
5 66 75 65 15
5 41 78 67 16
5 69 70 34 17
5 65 77 63 18
5 65 75 48 19
5 69 92 70 20
5 70 92 71 21
5 68 82 61 22
5 58 83 39 23
5 64 65 63 24
5 80 83 42 25
5 83 96 42 26
5 71 84 67 27
5 70 71 37 28
5 64 66 65 29
5 66 81 36 30
5 15 68 34 31
5 22 72 33 32
5 44 45 32 33
5 44 67 45 34
5 8 44 32 35
5 9 44 8 36
5 16 68 15 37
5 41 44 9 38
5 6 73 36 39
5 23 72 22 40
5 8 32 7 41
5 43 68 16 42
5 41 67 44 43
5 32 73 7 44
5 15 34 14 45
5 33 74 21 46
5 13 37 12 47
5 55 56 52 48
5 7 73 6 49
5 24 42 23 50
5 42 72 23 51
5 21 74 20 52
5 20 74 38 53
5 29 35 28 54
5 48 75 46 55
5 46 75 73 56
5 31 40 30 57
5 45 46 32 58
5 53 55 52 59
5 6 36 5 60
5 22 33 21 61
5 73 75 36 62
5 45 47 46 63
5 47 49 48 64
5 50 51 49 65
5 55 58 56 66
5 17 43 16 67
5 46 73 32 68
5 10 41 9 69
5 47 50 49 70
5 54 57 53 71
5 63 77 56 72
5 20 38 19 73
5 47 48 46 74
5 51 54 53 75
5 62 63 56 76
5 27 39 26 77
5 51 53 52 78
5 54 60 57 79
5 58 62 56 80
5 56 77 52 81
5 55 59 58 82
5 57 76 53 83
5 52 77 49 84
5 49 77 48 85
5 54 61 60 86
5 55 76 59 87
5 51 52 49 88
5 53 76 55 89
5 61 86 69 90
5 79 97 43 91
5 43 97 82 92
5 19 79 18 93
5 26 80 25 94
5 12 78 11 95
5 62 85 35 96
5 5 81 4 97
5 36 81 5 98
5 35 87 64 99
5 37 78 12 100
5 39 80 26 101
5 38 79 19 102
5 50 92 86 103
5 40 81 66 104
5 67 84 45 105
5 47 84 50 106
5 35 85 28 107
5 27 85 39 108
5 50 86 51 109
5 54 86 61 110
5 40 87 30 111
5 29 87 35 112
5 0 88 31 113
5 3 90 24 114
5 2 91 17 115
5 4 88 0 116
5 1 89 10 117
5 11 89 1 118
5 25 90 3 119
5 18 91 2 120
5 84 92 50 121
5 43 82 68 122
5 61 82 60 123
5 59 83 58 124
5 86 92 69 125
5 45 84 47 126
5 39 83 80 127
5 28 85 27 128
5 51 86 54 129
5 71 92 84 130
5 78 89 11 131
5 79 91 18 132
5 80 90 25 133
5 30 87 29 134
5 24 90 42 135
5 10 89 41 136
5 17 91 43 137
5 31 88 40 138
5 81 88 4 139
5 94 96 59 140
5 57 95 93 141
5 93 95 74 142
5 72 94 33 143
5 33 93 74 144
5 60 95 57 145
5 33 94 93 146
5 74 95 38 147
5 57 93 76 148
5 93 94 76 149
5 76 94 59 150
5 72 96 94 151
5 42 96 72 152
5 41 89 78 153
5 43 91 79 154
5 42 90 80 155
5 60 97 95 156
5 95 97 38 157
5 40 88 81 158
5 38 97 79 159
5 59 96 83 160
5 82 97 60 161
NPOIN= 98
0 0 0
1 0 1
1 1 2
0 1 3
0.1249999999999998 0 4
0.2499999999999998 0 5
0.3749999999999994 0 6
0.4999999999999988 0 7
0.6249999999999988 0 8
0.7499999999999991 0 9
0.8749999999999994 0 10
1 0.1249999999999998 11
1 0.2499999999999998 12
1 0.3749999999999994 13
1 0.4999999999999988 14
1 0.6249999999999988 15
1 0.7499999999999991 16
1 0.8749999999999994 17
0.8750000000000001 1 18
0.7500000000000002 1 19
0.6250000000000007 1 20
0.5000000000000011 1 21
0.3750000000000012 1 22
0.2500000000000009 1 23
0.1250000000000006 1 24
0 0.8750000000000001 25
0 0.7500000000000002 26
0 0.6250000000000007 27
0 0.5000000000000011 28
0 0.3750000000000012 29
0 0.2500000000000009 30
0 0.1250000000000006 31
0.5624999999999989 0.1082531754730547 32
0.4375000000000012 0.8917468245269454 33
0.8787408590474762 0.5688788296333867 34
0.1107249391757801 0.4359051368937709 35
0.3098918535636471 0.1070150426096389 36
0.8917468245269455 0.3124999999999996 37
0.6871891340151978 0.8955068139135626 38
0.1182626851017918 0.6692259739230109 39
0.1094851921887065 0.2033420655976694 40
0.8119703288398898 0.1058465056319792 41
0.1875000000000007 0.8917468245269449 42
0.879230463236312 0.8130906423764597 43
0.687411721473314 0.1078520638328755 44
0.6222382678779044 0.2148535068271494 45
0.4995397113129832 0.2162308769262824 46
0.5592159774975001 0.3228521478837691 47
0.4350399505186865 0.3232496269229517 48
0.4995115585047765 0.4327221170255846 49
0.6130571419252585 0.4327008212914639 50
0.55698961825214 0.5431506844806554 51
0.4365001961261518 0.5415315810733969 52
0.4989149690630472 0.6498774713089112 53
0.6204622660663317 0.6518781417959357 54
0.3748733700916652 0.6495569611339232 55
0.3123122610363024 0.5413164793658919 56
0.5624999999999974 0.757772228311381 57
0.2499999999999988 0.6495190528383259 58
0.3133176815799771 0.7523710733550596 59
0.685881360403847 0.7503445567222474 60
0.7558138451570435 0.6422638761822285 61
0.1875710532781786 0.5449975168824133 62
0.2499999999999996 0.433012701892217 63
0.1778732193721211 0.3232140480521861 64
0.3014840143016392 0.3178829819971145 65
0.2375681941678718 0.2131232466900193 66
0.7646672502010369 0.2290660061068141 67
0.891746824526945 0.6874999999999991 68
0.7774226411965255 0.5265527835888375 69
0.8628769461960453 0.4382017528180521 70
0.7819074654780986 0.3451910717748129 71
0.3125000000000009 0.891746824526945 72
0.4369885941461045 0.1080009076591809 73
0.5625000000000008 0.8917468245269453 74
0.3714895616679679 0.2148949910804582 75
0.4374999999999979 0.7577722283113808 76
0.3724746634145927 0.4316192480461928 77
0.8959004857786097 0.1879321779033213 78
0.8081462788186125 0.9070743513265296 79
0.103415783057214 0.8066379931605242 80
0.1911333561284868 0.1028269135594176 81
0.7923273839860845 0.74371081557924 82
0.2010889394410991 0.7499266504670594 83
0.6739266620600561 0.3223988573944038 84
0.08453570123751097 0.5625000000000009 85
0.671840423222428 0.54227566419079 86
0.07733755452579338 0.3125000000000011 87
0.09150635094610993 0.09150635094611009 88
0.90849364905389 0.09150635094610993 89
0.09150635094611009 0.9084936490538901 90
0.9084936490538901 0.90849364905389 91
0.7301718800130687 0.43455349184306 92
0.4999999999999993 0.8247595264191632 93
0.374999999999998 0.8247595264191614 94
0.6250000000000007 0.8247595264191651 95
0.260160832435939 0.8246932292153024 96
0.7489644266062752 0.8253550026152655 97
NMARK= 1
MARKER_TAG= void
MARKER_ELEMS= 32
3 0 4
3 4 5
3 5 6
3 6 7
3 7 8
3 8 9
3 9 10
3 10 1
3 1 11
3 11 12
3 12 13
3 13 14
3 14 15
3 15 16
3 16 17
3 17 2
3 2 18
3 18 19
3 19 20
3 20 21
3 21 22
3 22 23
3 23 24
3 24 3
3 3 25
3 25 26
3 26 27
3 27 28
3 28 29
3 29 30
3 30 31
3 31 0
OUTPUT_DIR = ../result
OUTPUT_FILE = phantom2d_csd_10MeV_fine
LOG_DIR = ../result/logs
MESH_FILE = phantom2D.su2
PROBLEM = PHANTOM2D
SOLVER = CSD_SN_FOKKERPLANCK_SOLVER
CONTINUOUS_SLOWING_DOWN = YES
HYDROGEN_FILE = ENDL_H.txt
OXYGEN_FILE = ENDL_O.txt
KERNEL = ISOTROPIC
CFL_NUMBER = 0.01
TIME_FINAL = 1.0
CLEAN_FLUX_MATRICES = NO
BC_DIRICHLET = ( dirichlet )
BC_NEUMANN = ( wall_low, wall_up )
QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
QUAD_ORDER = 10
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment