#include "SU21-base-defs.c"
#include "a7p2-base-defs.c"

// Generators for (trivial) finite subgroup $K$

// The group $K$ itself
#define GpKSiz 1
Mtx_t GpK[GpKSiz];

/* Set up the matrices for the finite subgroup $K$ */
void MakeGpK() {
  int Ix,ir,ic;

  memcpy(GpK[0],Id3,sizeof(Mtx_t));
  return;
}
  
// Other generators of our group $\Gamma$
#define NoGens 3
#define MaxNoGens 3

Mtx_t Gens[MaxNoGens];   // These will be the conjugated generators
Mtx_t Gens0[MaxNoGens]={ // These are the given generators
  {{1/18.0*(ZZ*ZZ+3*ZZ+3)*sqm7+1/6.0*(ZZ*ZZ+ZZ-3),  // This is U
      1/126.0*(11*ZZ*ZZ+18.0*ZZ-36)*sqm7+1/6.0*(-ZZ*ZZ-4*ZZ-2),
      1/126.0*(-2*ZZ*ZZ-9*ZZ-3)*sqm7+1/6.0*(ZZ+1)
    },
    {1/63.0*(-2*ZZ*ZZ-9*ZZ-3)*sqm7+1/3.0*(-ZZ-1),
     1/18.0*(ZZ*ZZ+3*ZZ+3)*sqm7+1/6.0*(-ZZ*ZZ-ZZ+1),
     -1/18.0*ZZ*ZZ*sqm7+1/6.0*(-ZZ*ZZ-2*ZZ)
    },
    {1/63.0*(-4*ZZ*ZZ-18.0*ZZ-27)*sqm7+1/3.0*(2*ZZ*ZZ+6*ZZ-1),
     1/63.0*(4*ZZ*ZZ+18.0*ZZ+6)*sqm7,
     1/18.0*(-2*ZZ*ZZ-6*ZZ+3)*sqm7-1/6.0
    }
  },
  {{1/126.0*(8*ZZ*ZZ+15*ZZ-9)*sqm7+1/6.0*(ZZ+5), // This is B
    1/63.0*(2*ZZ*ZZ+9*ZZ+3)*sqm7,
    1/252.0*(-11*ZZ*ZZ-18*ZZ-6)*sqm7+1/12.0*(-ZZ*ZZ-2*ZZ+2)
  },
   {1/126.0*(-2*ZZ*ZZ-9*ZZ-24)*sqm7+1/6.0*(2*ZZ*ZZ+5*ZZ),
    1/18.0*(-ZZ*ZZ-3*ZZ)*sqm7+1/6.0*(-ZZ*ZZ-3*ZZ+4),
    1/126.0*(-2*ZZ*ZZ-9*ZZ-3)*sqm7+1/6.0*(ZZ+1)
   },
   {1/63.0*(-2*ZZ*ZZ-9*ZZ-3)*sqm7+1/3.0*(-ZZ-1),
    1/126.0*(13*ZZ*ZZ+27*ZZ-33)*sqm7+1/6.0*(-ZZ*ZZ-3*ZZ+1),
    1/126.0*(-ZZ*ZZ+6*ZZ+9)*sqm7+1/6.0*(ZZ*ZZ+2*ZZ+3)
   }
  },
  {{1/126.0*(2*ZZ*ZZ+9*ZZ+24)*sqm7+1/6.0*(ZZ-2), // This is A
    1/126.0*(-8*ZZ*ZZ-15*ZZ+9)*sqm7+1/6.0*(-ZZ+1),
    1/126.0*(2*ZZ*ZZ+9*ZZ+3)*sqm7+1/6.0*(ZZ+3)
  },
   {1/63.0*(-4*ZZ*ZZ-18*ZZ-6)*sqm7+2/3.0,
    1/126.0*(-4*ZZ*ZZ-18*ZZ+15)*sqm7-1/2.0,
    1/18.0*(ZZ*ZZ+3*ZZ)*sqm7+1/6.0*(ZZ*ZZ+3*ZZ+2)
   },
   {1/63.0*(ZZ*ZZ-6*ZZ-9)*sqm7+1/3.0*(-ZZ*ZZ-2*ZZ+3),
    1/63.0*(2*ZZ*ZZ+9*ZZ+3)*sqm7+1/3.0*(-ZZ+1),
    1/126.0*(2*ZZ*ZZ+9*ZZ+24)*sqm7+1/6.0*(-ZZ-4)
   }
  }
};

#include "a7p2-base.c"