Initial commit from Polytechnique Montreal

1 files changed, 243 insertions, 0 deletions

diff --git a/Trivac/data/DIFtst_proc/pertdiff.c2m b/Trivac/data/DIFtst_proc/pertdiff.c2m
new file mode 100755
index 0000000..2196541
--- /dev/null
+++ b/Trivac/data/DIFtst_proc/pertdiff.c2m
@@ -0,0 +1,243 @@
+*-----
+*  GPT TEST pertdiff
+*-----
+LINKED_LIST IAEA MACRO_REF MACRO_0 DMACRO TRACK SYSTEM_REF FLUX DSYSTEM
+            DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 EDIT ADFLUX MACRO ;
+MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: DELETE: DELTA: GPTFLU: OUT:
+       END: ADD: UTL: STAT: ;
+PROCEDURE assertS ;
+*
+IAEA := GEO: :: CAR2D 9 9
+           EDIT 2
+           X- DIAG X+ VOID
+           Y- SYME Y+ DIAG
+           MIX  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
+           MESHX 0.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0
+      ;
+*
+* REFERENCE CASE:
+MACRO_REF := MAC: ::
+ EDIT 2 NGRO 2 NMIX 4 NIFI 1
+ READ INPUT
+ MIX     1
+      DIFF  1.900000E+00  4.400000E-01
+     TOTAL  3.012000E-02  8.303201E-02
+    NUSIGF  0.000000E+00  1.650000E-01
+  H-FACTOR  0.000000E+00  1.650000E-01
+      SCAT  1 1 0.0 2 2 0.0   0.200000E-01
+ MIX     2
+      DIFF  1.500000E+00  4.000000E-01
+     TOTAL  3.012000E-02  8.503199E-02
+    NUSIGF  0.000000E+00  1.350000E-01
+  H-FACTOR  0.000000E+00  1.350000E-01
+      SCAT  1 1 0.0 2 2 0.0   0.200000E-01
+ MIX     3
+      DIFF  1.500000E+00  4.000000E-01
+     TOTAL  2.912000E-02  1.260320E-01
+    NUSIGF  0.000000E+00  1.350000E-01
+  H-FACTOR  0.000000E+00  1.350000E-01
+      SCAT  1 1 0.0 2 2 0.0   0.200000E-01
+ MIX     4
+      DIFF  2.000000E+00  3.000000E-01
+     TOTAL  4.016000E-02  1.002400E-02
+      SCAT  1 1 0.0 2 2 0.0   0.400000E-01
+ ;
+*
+* UNPERTURBED CASE:
+MACRO_0 := MAC: ::
+ EDIT 2 NGRO 2 NMIX 4 NIFI 1
+ READ INPUT
+ MIX     1
+      DIFF  1.500000E+00  4.000000E-01
+     TOTAL  3.012000E-02  8.003199E-02
+    NUSIGF  0.000000E+00  1.350000E-01
+  H-FACTOR  0.000000E+00  1.350000E-01
+      SCAT  1 1 0.0 2 2 0.0  0.200000E-01
+ MIX     2
+      DIFF  1.500000E+00  4.000000E-01
+     TOTAL  3.012000E-02  8.503199E-02
+    NUSIGF  0.000000E+00  1.350000E-01
+  H-FACTOR  0.000000E+00  1.350000E-01
+      SCAT  1 1 0.0 2 2 0.0  0.200000E-01
+ MIX     3
+      DIFF  1.500000E+00  4.000000E-01
+     TOTAL  3.012000E-02  1.300320E-01
+    NUSIGF  0.000000E+00  1.350000E-01
+  H-FACTOR  0.000000E+00  1.350000E-01
+      SCAT  1 1 0.0 2 2 0.0  0.200000E-01
+ MIX     4
+      DIFF  2.000000E+00  3.000000E-01
+     TOTAL  4.016000E-02  1.002400E-02
+      SCAT  1 1 0.0 2 2 0.0  0.400000E-01
+ ;
+*
+* PERTURBATION:
+DMACRO := MAC: ::
+ EDIT 2 NGRO 2 NMIX 4 NIFI 1
+ READ INPUT
+ MIX     1
+      DIFF  4.000000E-01  4.000000E-02
+     TOTAL  0.000000E+00  3.000000E-03
+    NUSIGF  0.000000E+00  3.000000E-02
+  H-FACTOR  0.000000E+00  3.000000E-02
+      SCAT  1 1 0.0 2 2 0.0  1.0E-10
+ MIX     2
+      SCAT  1 1 0.0 2 2 0.0  1.0E-10
+ MIX     3
+     TOTAL -9.999999E-04 -4.000001E-03
+      SCAT  1 1 0.0 2 2 0.0  1.0E-10
+ MIX     4
+      SCAT  1 1 0.0 2 2 0.0  1.0E-10
+ ;
+*----
+* MCFD1 case, Livolant iteration
+*----
+TRACK := TRIVAT: IAEA ::
+      TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)'
+      EDIT 5 MAXR 1156
+      MCFD 1 ;
+SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_REF TRACK :: EDIT 2 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ;
+FLUX := DELETE: FLUX ;
+*
+SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ;
+DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ;
+DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ;
+DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 EXPLICIT
+         FROM-TO 1 1 ;
+assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ;
+ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 IMPLICIT
+          FROM-TO 1 1 ;
+assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+* Reset the perturbation flag of DSYSTEM to 0 so that it can be added.
+DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ;
+SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ;
+MACRO := ADD: MACRO_0 DMACRO ;
+*
+STAT: SYSTEM_REF SYSTEM_0 ;
+STAT: MACRO_REF MACRO ;
+FLUX2 :=  FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ;
+assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ;
+EDIT := OUT: FLUX2 TRACK MACRO IAEA ::
+       EDIT 2 INTG
+      1  2  3  4  5  6  7  8  0
+         9 10 11 12 13 14 15  0
+           16 17 18 19 20 21  0
+              22 23 24 25  0  0
+                 26 27 28  0  0
+                    29  0  0  0
+                        0  0  0
+                           0  0
+                              0
+      ;
+TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2
+EDIT ADFLUX MACRO
+:= DELETE:
+TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2
+EDIT ADFLUX MACRO ;
+*----
+* RAVIART-THOMAS DUAL 1 2 case, Livolant iteration
+*----
+TRACK := TRIVAT: IAEA ::
+      TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)'
+      EDIT 5 MAXR 1156
+      DUAL 1 2 ;
+SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_REF TRACK :: EDIT 2 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ;
+FLUX := DELETE: FLUX ;
+*
+SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ;
+DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ;
+DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ;
+DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 EXPLICIT
+         FROM-TO 1 1 ;
+assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ;
+ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 IMPLICIT
+          FROM-TO 1 1 ;
+assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+* Reset the perturbation flag of DSYSTEM to 0 so that it can be added.
+DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ;
+SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ;
+MACRO := ADD: MACRO_0 DMACRO ;
+*
+STAT: SYSTEM_REF SYSTEM_0 ;
+STAT: MACRO_REF MACRO ;
+FLUX2 :=  FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ;
+assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ;
+EDIT := OUT: FLUX2 TRACK MACRO IAEA ::
+       EDIT 2 INTG
+      1  2  3  4  5  6  7  8  0
+         9 10 11 12 13 14 15  0
+           16 17 18 19 20 21  0
+              22 23 24 25  0  0
+                 26 27 28  0  0
+                    29  0  0  0
+                        0  0  0
+                           0  0
+                              0
+      ;
+TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2
+EDIT ADFLUX MACRO
+:= DELETE:
+TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2
+EDIT ADFLUX MACRO ;
+*----
+* RAVIART-THOMAS DUAL 1 2 case, GMRES iteration
+*----
+TRACK := TRIVAT: IAEA ::
+      TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)'
+      EDIT 5 MAXR 1156
+      DUAL 1 2 ;
+SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_REF TRACK :: EDIT 2 ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ;
+FLUX := DELETE: FLUX ;
+*
+SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ;
+FLUX :=  FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ;
+DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ;
+DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ;
+DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 GMRES 10 EXPLICIT
+         FROM-TO 1 1 ;
+assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ;
+ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 GMRES 10 IMPLICIT
+          FROM-TO 1 1 ;
+assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ;
+* Reset the perturbation flag of DSYSTEM to 0 so that it can be added.
+DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ;
+SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ;
+MACRO := ADD: MACRO_0 DMACRO ;
+*
+STAT: SYSTEM_REF SYSTEM_0 ;
+STAT: MACRO_REF MACRO ;
+FLUX2 :=  FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ;
+assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ;
+EDIT := OUT: FLUX2 TRACK MACRO IAEA ::
+       EDIT 2 INTG
+      1  2  3  4  5  6  7  8  0
+         9 10 11 12 13 14 15  0
+           16 17 18 19 20 21  0
+              22 23 24 25  0  0
+                 26 27 28  0  0
+                    29  0  0  0
+                        0  0  0
+                           0  0
+                              0
+      ;
+ECHO "test pertdiff completed" ;
+END: ;