Initial commit from Polytechnique Montreal

1 files changed, 173 insertions, 0 deletions

diff --git a/Trivac/data/DIFtst_proc/vv1k3d.c2m b/Trivac/data/DIFtst_proc/vv1k3d.c2m
new file mode 100755
index 0000000..683ada5
--- /dev/null
+++ b/Trivac/data/DIFtst_proc/vv1k3d.c2m
@@ -0,0 +1,173 @@
+*------------------------
+* Benchmark VV1K3D
+*
+* Reference: A. Bernal, A. Hebert, J. E. Roman, R. Miro and G. Verdu,
+* "A Krylov-Schur solution of the eigenvalue problem for the neutron
+* diffusion equation discretized with the Raviart-Thomas method," J. of
+* Nuclear Science and Technology, Vol. 54, No. 10, 1085-1094 (2017).
+*------------------------
+LINKED_LIST GEOM XSEC TRACK SYSTEM FLUX RESU ; 
+MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: OUT: END: ABORT: ;
+PROCEDURE assertS ;
+
+GEOM := GEO: :: 
+                HEXZ 169 10
+                EDIT 2
+                Z- ZERO Z+ ZERO HBC COMPLETE ZERO
+                MESHZ 0. 20. 40. 60. 80. 100. 120. 140. 160. 180. 200.
+                SIDE 13.6255 
+                SPLITL 1
+                MIX 
+                ! PLANE 1
+                2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 2
+                2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 3
+                2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 4
+                2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 5
+                2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 6
+                2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 7
+                2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 8
+                2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 9
+                2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ! PLANE 10
+                2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1
+                3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1
+                1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4
+                1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4
+                5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3
+                5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
+                5 5 5 5 5 5 5 5
+                ;
+
+XSEC := MAC: :: 
+                EDIT 2
+                NGRO 2 NIFI 1 NMIX 5
+                READ INPUT 
+                MIX 1
+                   TOTAL 2.48836E-02  6.73049E-2  
+                   NUSIGF 4.81619E-3 8.46154E-2  
+                   DIFF   1.38320 3.86277E-1  
+                   CHI 1. 0.
+                   SCAT 1 1 0.
+                        1 1 1.64977E-2
+                MIX 2
+                   TOTAL 2.62865E-02  8.10328E-2 
+                   NUSIGF  4.66953E-3 8.52264E-2 
+                   DIFF    1.38299  3.89403E-1 
+                   CHI 1. 0.
+                   SCAT 1 1 0.
+                        1 1 1.47315E-2    
+                MIX 3
+                   TOTAL  2.45662E-02 8.44801E-2 
+                   NUSIGF  6.04889E-3  1.19428E-1 
+                   DIFF    1.39522  3.86225E-1  
+                   CHI 1. 0.
+                   SCAT 1 1 0.
+                        1 1 1.56219E-2 
+                MIX 4
+                   TOTAL 2.60117E-02   9.89671E-2 
+                   NUSIGF 5.91507E-3 1.20497E-1  
+                   DIFF   1.39446  3.87723E-1   
+                   CHI 1. 0.
+                   SCAT 1 1 0.
+                        1 1 1.40185E-2 
+                MIX 5
+                   TOTAL 2.46141E-02  8.93878E-2
+                   NUSIGF 6.40256E-3 1.29281E-1 
+                   DIFF   1.39506 3.84492E-1  
+                   CHI 1. 0.
+                   SCAT 1 1 0.
+                        1 1 1.54981E-2
+                   ;
+
+TRACK := TRIVAT: GEOM ::
+               EDIT 2
+               TITL 'VV1K3D'
+               MAXR 10000 DUAL 2 1
+               ;
+
+SYSTEM := TRIVAA: XSEC TRACK :: EDIT 3 ;
+
+FLUX := FLUD: SYSTEM TRACK :: 
+                EDIT 2 ADJ IRAM 3 6 EXTE 30 1.0E-8
+                ;  
+
+RESU := OUT: FLUX SYSTEM TRACK XSEC GEOM :: 
+                 EDIT 2
+                 MODE 1 POWR 1.0 INTG MIX
+                 ;
+
+assertS FLUX :: 'K-EFFECTIVE' 1 1.005450 ;
+ECHO "test vv1k3d completed" ;
+END: ;