Initial commit from Polytechnique Montreal

1 files changed, 133 insertions, 0 deletions

diff --git a/Dragon/data/tmacro_proc/TCM49.c2m b/Dragon/data/tmacro_proc/TCM49.c2m
new file mode 100644
index 0000000..3a0afa6
--- /dev/null
+++ b/Dragon/data/tmacro_proc/TCM49.c2m
@@ -0,0 +1,133 @@
+*----
+*  TEST CASE TCM49
+*  NSE 2D CARTESIAN PROBLEM MODIFIED FOR 3D SN REGRESSION TEST
+*  ONE-GROUP MACROLIB-DEFINED CROSS SECTIONS
+*
+*  Ref: A. Hebert, "The Search for Superconvergence in Spherical
+*  Harmonics Approximations", Nucl. Sci. Eng, 2006.
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM MACRO TRACK SYSTEM FLUX EDIT ;
+MODULE GEO: MAC: SNT: ASM: FLU: DELETE: END: ;
+INTEGER s := 2 ;
+PROCEDURE assertS ;
+*
+*--------
+*GEOMETRY
+*--------
+GEOM := GEO: :: CAR3D  5 5 1
+           X- REFL X+ VOID
+           Y- REFL Y+ VOID
+           Z- VOID Z+ 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
+           MESHZ 0.0 1000.0
+           SPLITX <<s>> <<s>> <<s>> <<s>> <<s>>
+           SPLITY <<s>> <<s>> <<s>> <<s>> <<s>>
+           ;
+*--------
+*MACROLIB
+*--------
+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
+ ;
+*--------
+*TRACKING
+*--------
+*----
+* HIGH-ORDER DIAMOND DIFFERENCE
+*----
+!
+!FLAT ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 DIAM 0 SN 6 SCAT 2 QUAD 1 DSA 10 0 2 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9178295 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!LINEAR ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 DIAM 1 SN 6 SCAT 2 QUAD 1 DSA 10 1 2 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9817208 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!PARABOLIC ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 DIAM 2 SN 6 SCAT 2 QUAD 1 DSA 10 2 2 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9829153 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+*----
+* DISCONTINUOUS GALERKIN
+*----
+!
+!LINEAR ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC DG LINEAR BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 SCHM 2 DIAM 1 SN 6 SCAT 2 QUAD 1 DSA 10 1 2 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9509456 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!PARABOLIC ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC DG PARABOLIC BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 SCHM 2 DIAM 2 SN 6 SCAT 2 QUAD 1 DSA 10 2 2 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9830175 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!CUBIC ORDER
+TRACK := SNT: GEOM ::
+      TITLE '3D ANISOTROPIC DG CUBIC BENCHMARK, 1 GROUP'
+      EDIT 2 MAXR 5000 SCHM 2 DIAM 3 SN 4 SCAT 2 QUAD 1 ;
+SYSTEM := ASM: MACRO TRACK ::
+      EDIT 1 ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+      EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9827218 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+
+ECHO "test TCM49 completed" ;
+END: ;
+QUIT "LIST" .