Initial commit from Polytechnique Montreal

1 files changed, 156 insertions, 0 deletions

diff --git a/Dragon/data/tmacro_proc/TCM50.c2m b/Dragon/data/tmacro_proc/TCM50.c2m
new file mode 100644
index 0000000..5f5cc12
--- /dev/null
+++ b/Dragon/data/tmacro_proc/TCM50.c2m
@@ -0,0 +1,156 @@
+*----TEST CASE TCM50--------*
+*  REGRESSION TEST FOR SN 3D CARTESIAN GEOMETRY
+*  IAEA 3D BENCHMARK -- SN THEORY (ISOTROPIC SCATTERING)
+*  MACROLIB-DEFINED CROSS SECTIONS
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST IAEA3D MACRO MACROT TRACK SYSTEM FLUX EDIT ;
+MODULE GEO: MAC: T: SNT: ASM: FLU: DELETE: END: ;
+PROCEDURE assertS ;
+
+IAEA3D := GEO: :: CAR3D 9 9 4
+          EDIT 0
+          X- DIAG  X+ VOID 
+          Y- REFL  Y+ DIAG 
+          Z- VOID  Z+ VOID 
+          MESHX 10.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0  
+          MESHZ 0.0 20.0 280.0 360.0 380.0 
+          SPLITZ 1 2 1 1 
+          ! PLANE NB 1
+          MIX 4 4 4 4 4 4 4 4 4 
+                4 4 4 4 4 4 4 4 
+                  4 4 4 4 4 4 4 
+                    4 4 4 4 4 4 
+                      4 4 4 4 0 
+                        4 4 4 0 
+                          4 0 0 
+                            0 0 
+                              0 
+              ! PLANE NB 2
+              3 2 2 2 3 2 2 1 4 
+                2 2 2 2 2 2 1 4 
+                  2 2 2 2 1 1 4 
+                    2 2 2 1 4 4 
+                      3 1 1 4 0 
+                        1 4 4 0 
+                          4 0 0 
+                            0 0 
+                              0 
+              ! PLANE NB 3
+              3 2 2 2 3 2 2 1 4 
+                2 2 2 2 2 2 1 4 
+                  3 2 2 2 1 1 4 
+                    2 2 2 1 4 4 
+                      3 1 1 4 0 
+                        1 4 4 0 
+                          4 0 0 
+                            0 0 
+                              0 
+              ! PLANE NB 4
+              4 4 4 4 4 4 4 4 4 
+                4 4 4 4 4 4 4 4 
+                  4 4 4 4 4 4 4 
+                    4 4 4 4 4 4 
+                      4 4 4 4 0 
+                        4 4 4 0 
+                          4 0 0 
+                            0 0 
+                              0 
+           ;
+MACRO := MAC: ::
+ EDIT 0 NGRO 2 NMIX 4 NIFI 1 
+ READ INPUT
+ MIX     1
+     TOTAL  2.222E-01  8.3333E-01
+    NUSIGF  0.000E+00  1.3500E-01
+       CHI  1.000E+00  0.0000E+00
+  H-FACTOR  0.000E+00  1.3500E-01
+      SCAT  1 1 0.1921 2 2 0.7533 0.02
+ MIX     2
+     TOTAL  2.222E-01  8.3333E-01
+    NUSIGF  0.000E+00  1.3500E-01
+       CHI  1.000E+00  0.0000E+00
+  H-FACTOR  0.000E+00  1.3500E-01
+      SCAT  1 1 0.1921 2 2 0.7483 0.02
+ MIX     3
+     TOTAL  2.222E-01  8.3333E-01
+    NUSIGF  0.000E+00  1.35000E-01
+       CHI  1.000E+00  0.00000E+00
+  H-FACTOR  0.000E+00  1.35000E-01
+      SCAT  1 1 0.1921 2 2 0.7033 0.02
+ MIX     4
+     TOTAL  1.667E-01  1.1111
+      SCAT  1 1 0.1265 2 2 1.1011 0.04
+ ;
+!*----
+!* HIGH-ORDER DIAMOND DIFFERENCE
+!*----
+!FLAT ORDER NO DSA
+TRACK := SNT: IAEA3D ::
+TITLE 'IAEA3D BENCHMARK,ISOTROPIC SCATTERING NDSA, 2 GROUPS'
+EDIT 1 MAXR 50000 DIAM 0 SN 6 SCAT 1 QUAD 1 NDSA GMRES 10 MAXI 50 ;
+SYSTEM := ASM: MACRO TRACK ::
+ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.025034 ;
+SYSTEM FLUX := DELETE: SYSTEM FLUX ;
+!
+!FLAT ORDER (ADJOINT CALCULATION)
+MACROT := T: MACRO ;
+SYSTEM := ASM: MACROT TRACK ::
+ARM ;
+FLUX :=  FLU: MACROT TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.025038 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!LINEAR ORDER NO DSA
+TRACK := SNT: IAEA3D ::
+TITLE 'IAEA3D BENCHMARK,ISOTROPIC HODD 1 NDSA, 2 GROUPS'
+MAXR 50000 DIAM 1 SN 6 SCAT 1 QUAD 1 NDSA GMRES 10 MAXI 50 ;
+SYSTEM := ASM: MACRO TRACK ::
+ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.025936 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!LINEAR ORDER WITH DSA
+TRACK := SNT: IAEA3D ::
+TITLE 'IAEA3D BENCHMARK,ISOTROPIC HODD 1 DSA, 2 GROUPS'
+MAXR 50000 DIAM 1 SN 6 SCAT 1 QUAD 1 DSA 10 0 2 GMRES 10 MAXI 50 ;
+SYSTEM := ASM: MACRO TRACK ::
+ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.025935 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+*----
+* DISCONTINUOUS GALERKIN
+*----
+!LINEAR ORDER NO DSA
+TRACK := SNT: IAEA3D ::
+TITLE 'IAEA3D BENCHMARK,ISOTROPIC DG NDSA 1, 2 GROUPS'
+MAXR 50000 SCHM 2 DIAM 1 SN 6 SCAT 1 QUAD 1 NDSA GMRES 10 MAXI 50 ;
+SYSTEM := ASM: MACRO TRACK ::
+ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.027354 ;
+TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
+!
+!LINEAR ORDER WITH DSA
+TRACK := SNT: IAEA3D ::
+TITLE 'IAEA3D BENCHMARK,ISOTROPIC DG 1 DSA, 2 GROUPS'
+MAXR 50000 SCHM 2 DIAM 1 SN 6 SCAT 1 QUAD 1 DSA 10 0 2 GMRES 10 MAXI 50 ;
+SYSTEM := ASM: MACRO TRACK ::
+ARM ;
+FLUX :=  FLU: MACRO TRACK SYSTEM ::
+EDIT 2 TYPE K EXTE 100 5.0E-8 ;
+assertS FLUX :: K-EFFECTIVE 1 1.027353 ;
+
+ECHO "test TCM50 completed" ;
+END: ;
+QUIT "LIST" .