blob: 788dd6d41ceb3ed53dd17d4eb3a49ce19f3b1d67 (
plain)
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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
*----
* TEST CASE TCM56
* CNS 2D CARTESIAN PROBLEM
* ONE-GROUP MACROLIB-DEFINED CROSS SECTIONS
*
* Ref: A. Calloo and A. Hebert, CNS 2019
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST GEOM MACRO TRACK SYSTEM FLUX EDIT ;
MODULE GEO: MAC: SNT: ASM: FLU: UTL: END: ;
PROCEDURE assertS ;
*
* s = region split
* sch = scheme parameter
* spa = spatial order
* and = sn order
* sca = scattering
*
INTEGER s := 6 ;
INTEGER sch := 2 ;
INTEGER spa := 3 ;
INTEGER ang := 18 ;
INTEGER sca := 2 ;
*
GEOM := GEO: :: CAR2D 5 5
X- REFL X+ VOID
Y- REFL Y+ VOID
MIX 1 1 2 2 3
1 1 2 2 3
2 2 2 2 3
2 2 2 3 3
3 3 3 3 0
MESHX 0.0 40.0 80.0 120.0 160.0 200.0
MESHY 0.0 40.0 80.0 120.0 160.0 200.0
SPLITX <<s>> <<s>> <<s>> <<s>> <<s>>
SPLITY <<s>> <<s>> <<s>> <<s>> <<s>>
;
MACRO := MAC: ::
EDIT 2 NGRO 1 NMIX 3 NIFI 1 ANIS 2
READ INPUT
MIX 1
TOTAL 0.025
NUSIGF 0.0155
CHI 1.0
SCAT 1 1 0.013
1 1 0.0
MIX 2
TOTAL 0.025
SCAT 1 1 0.024
1 1 0.006
MIX 3
TOTAL 0.075
SCAT 1 1 0.0
1 1 0.0
;
TRACK := SNT: GEOM ::
TITLE '2D ANISOTROPIC CNS DG BENCHMARK, 1 GROUP'
EDIT 2 MAXR 25000
SCHM <<sch>>
DIAM <<spa>>
SN <<ang>>
SCAT <<sca>>
MAXI 300 EPSI 1.E-5
QUAD 2 (*Level-Symmetric*) ;
SYSTEM := ASM: MACRO TRACK :: EDIT 2 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 1E-6 ;
assertS FLUX :: 'K-EFFECTIVE' 1 9.922985E-01 ;
ECHO "test TCM56 completed" ;
END: ;
QUIT "LIST" .
|
