Quelle  mat2   Sprache: C

#############################################################################
##
#A  testmat2               EDIM-mini-package                     Frank Lübeck
##
##  
#Y  Copyright (C) 1999  Lehrstuhl D f\"ur Mathematik, RWTH Aachen
##  
##  Smaller example matrix for the EDIM package
##  

Print("Reading 34x34 integer matrix 'mat2' with elementary divisors 'eldiv2'.\n");

mat2 :=
[ [ 10, 2, 2, 2, -4, -2, -2, 0, 0, 8, -4, 2, 2, 2, 0, 0, -10, 4, 4, 2, 10, 
      -4, -4, -2, -2, 0, 0, 0, 0, 0, 2, -4, -2, 0 ], 
  [ 2, 16, -8, 4, -2, -4, 0, 4, -2, -4, -2, 4, 4, 0, -4, 2, -4, 12, 0, 0, 4, 
      -12, 6, 0, 0, -8, 4, -4, 2, 4, 0, 0, 0, -4 ], 
  [ 2, -8, 16, -2, -2, 2, 0, -2, 4, 8, -2, -2, -2, 0, 2, -4, -4, -6, 0, 0, 4, 
      6, -12, 0, 0, 4, -8, 2, 2, -2, 0, 0, 0, 2 ], 
  [ 2, 4, -2, 16, -8, 0, -4, 2, 2, -4, 0, 0, 0, 4, -2, -2, -4, 8, 2, 4, 0, 0, 
      0, 0, 0, -4, 2, -8, 4, 2, 0, -12, 0, 0 ], 
  [ -4, -2, -2, -8, 16, 0, 2, -4, 2, -4, 0, 0, 0, -2, 4, -2, 8, -4, -4, -2, 
      0, 0, 0, 0, 0, 2, 2, 4, -8, 2, 0, 6, 0, 0 ], 
  [ -2, -4, 2, 0, 0, 18, 6, -12, 6, -12, 0, 0, -6, -2, 4, -2, 2, -4, 4, 2, 
      -2, 4, -2, 2, 2, -8, 4, -4, 2, 4, -6, 0, 0, 12 ], 
  [ -2, 0, 0, -4, 2, 6, 18, -6, -6, -12, -4, 2, -2, -6, 2, 2, 2, 0, 0, 0, -2, 
      0, 0, 0, 6, -4, 2, -8, 4, 2, -2, 4, 2, 4 ], 
  [ 0, 4, -2, 2, -4, -12, -6, 24, -12, 0, 2, -4, 4, 2, -8, 4, 0, 4, -4, -2, 
      0, -4, 2, 4, -2, 4, -2, 2, 2, -8, 4, -2, 2, -8 ], 
  [ 0, -2, 4, 2, 2, 6, -6, -12, 24, 0, 2, 2, -2, 2, 4, -8, 0, -2, 2, -2, 0, 
      2, -4, -2, -2, -2, -2, 2, -4, 4, -2, -2, 2, 4 ], 
  [ 8, -4, 8, -4, -4, -12, -12, 0, 0, 48, -4, -4, 4, 4, 0, 0, -8, -4, 4, -4, 
      8, 4, -8, -4, -4, 12, -12, 12, 0, 0, 4, 4, -4, -8 ], 
  [ -4, -2, -2, 0, 0, 0, -4, 2, 2, -4, 16, -8, 0, 0, 0, 0, 0, 0, -4, -2, -8, 
      4, 4, 2, -4, 2, 2, 4, -8, 2, 0, 0, 0, 0 ], 
  [ 2, 4, -2, 0, 0, 0, 2, -4, 2, -4, -8, 16, 0, 0, 0, 0, 0, 0, 2, 4, 4, -8, 
      4, -4, 2, -4, 2, -2, 4, -4, 0, 0, 0, 0 ], 
  [ 2, 4, -2, 0, 0, -6, -2, 4, -2, 4, 0, 0, 30, 10, -20, 10, -10, 20, 0, 0, 
      10, -20, 10, -10, -10, 0, 0, 0, 0, 0, 2, -4, -2, -4 ], 
  [ 2, 0, 0, 4, -2, -2, -6, 2, 2, 4, 0, 0, 10, 30, -10, -10, -10, 0, 0, 0, 
      -6, -4, 8, -2, -2, 0, 0, 0, 0, 0, 10, -20, -10, -4 ], 
  [ 0, -4, 2, -2, 4, 4, 2, -8, 4, 0, 0, 0, -20, -10, 40, -20, 0, -20, 0, 0, 
      -4, 4, -2, 8, 8, 0, 0, 0, 0, 0, -4, 2, 4, 12 ], 
  [ 0, 2, -4, -2, -2, -2, 2, 4, -8, 0, 0, 0, 10, -10, -20, 40, 0, 10, 0, 0, 
      8, -2, -8, -4, -4, 0, 0, 0, 0, 0, -4, 2, 4, -6 ], 
  [ -10, -4, -4, -4, 8, 2, 2, 0, 0, -8, 0, 0, -10, -10, 0, 0, 50, -20, 0, 0, 
      -10, 12, 0, -6, 6, 0, 0, 0, 0, 0, -6, 12, -6, 8 ], 
  [ 4, 12, -6, 8, -4, -4, 0, 4, -2, -4, 0, 0, 20, 0, -20, 10, -20, 60, 0, 0, 
      12, -20, 10, -4, -16, 0, 0, 0, 0, 0, 8, -4, 4, -8 ], 
  [ 4, 0, 0, 2, -4, 4, 0, -4, 2, 4, -4, 2, 0, 0, 0, 0, 0, 0, 16, 8, 0, 0, 0, 
      0, 0, -4, -4, -2, 4, 2, -4, -4, -2, 0 ], 
  [ 2, 0, 0, 4, -2, 2, 0, -2, -2, -4, -2, 4, 0, 0, 0, 0, 0, 0, 8, 16, 0, 0, 
      0, 0, 0, -2, -2, -4, 2, -2, -2, -8, -4, 0 ], 
  [ 10, 4, 4, 0, 0, -2, -2, 0, 0, 8, -8, 4, 10, -6, -4, 8, -10, 12, 0, 0, 50, 
      -20, -20, -10, -10, 0, 0, 0, 0, 0, -6, 12, 6, 4 ], 
  [ -4, -12, 6, 0, 0, 4, 0, -4, 2, 4, 4, -8, -20, -4, 4, -2, 12, -20, 0, 0, 
      -20, 60, -30, 0, 0, 0, 0, 0, 0, 0, 4, 16, -4, -4 ], 
  [ -4, 6, -12, 0, 0, -2, 0, 2, -4, -8, 4, 4, 10, 8, -2, -8, 0, 10, 0, 0, 
      -20, -30, 60, 0, 0, 0, 0, 0, 0, 0, 4, -20, -4, 2 ], 
  [ -2, 0, 0, 0, 0, 2, 0, 4, -2, -4, 2, -4, -10, -2, 8, -4, -6, -4, 0, 0, 
      -10, 0, 0, 30, 0, 0, 0, 0, 0, 0, 2, -4, 10, -8 ], 
  [ -2, 0, 0, 0, 0, 2, 6, -2, -2, -4, -4, 2, -10, -2, 8, -4, 6, -16, 0, 0, 
      -10, 0, 0, 0, 30, 0, 0, 0, 0, 0, -10, -4, -2, 4 ], 
  [ 0, -8, 4, -4, 2, -8, -4, 4, -2, 12, 2, -4, 0, 0, 0, 0, 0, 0, -4, -2, 0, 
      0, 0, 0, 0, 24, -12, 12, -6, -12, 0, 0, 0, -8 ], 
  [ 0, 4, -8, 2, 2, 4, 2, -2, -2, -12, 2, 2, 0, 0, 0, 0, 0, 0, -4, -2, 0, 0, 
      0, 0, 0, -12, 24, -6, -6, 6, 0, 0, 0, 4 ], 
  [ 0, -4, 2, -8, 4, -4, -8, 2, 2, 12, 4, -2, 0, 0, 0, 0, 0, 0, -2, -4, 0, 0, 
      0, 0, 0, 12, -6, 24, -12, -6, 0, 0, 0, -4 ], 
  [ 0, 2, 2, 4, -8, 2, 4, 2, -4, 0, -8, 4, 0, 0, 0, 0, 0, 0, 4, 2, 0, 0, 0, 
      0, 0, -6, -6, -12, 24, -6, 0, 0, 0, 2 ], 
  [ 0, 4, -2, 2, 2, 4, 2, -8, 4, 0, 2, -4, 0, 0, 0, 0, 0, 0, 2, -2, 0, 0, 0, 
      0, 0, -12, 6, -6, -6, 24, 0, 0, 0, 4 ], 
  [ 2, 0, 0, 0, 0, -6, -2, 4, -2, 4, 0, 0, 2, 10, -4, -4, -6, 8, -4, -2, -6, 
      4, 4, 2, -10, 0, 0, 0, 0, 0, 30, 0, 0, -12 ], 
  [ -4, 0, 0, -12, 6, 0, 4, -2, -2, 4, 0, 0, -4, -20, 2, 2, 12, -4, -4, -8, 
      12, 16, -20, -4, -4, 0, 0, 0, 0, 0, 0, 60, 0, 0 ], 
  [ -2, 0, 0, 0, 0, 0, 2, 2, 2, -4, 0, 0, -2, -10, 4, 4, -6, 4, -2, -4, 6, 
      -4, -4, 10, -2, 0, 0, 0, 0, 0, 0, 0, 30, 0 ], 
  [ 0, -4, 2, 0, 0, 12, 4, -8, 4, -8, 0, 0, -4, -4, 12, -6, 8, -8, 0, 0, 4, 
      -4, 2, -8, 4, -8, 4, -4, 2, 4, -12, 0, 0, 36 ] ];

eldiv2 := 
[ 2, 2, 2, 2, 2, 2, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 12, 12, 12, 
  48, 48, 48, 48, 48, 48, 144, 1008, 1008, 15120 ];