#include #include #include #include int readPaulicoeffs(int* alpha, int* beta, int* gamma, int* delta, int numqubits); complex double Gausssum1d(int quadraticcoeff, int linearcoeff); complex double Kroneck(int arg); // order of matrix elements is [row][column]!!! int main() { int i, j; int x, y, xp, yp, xpp, ypp, xppp, yppp; int N; // number of qubits scanf("%d", &N); if(N%12!=0) { printf("Error: N must be a multiple of 12!\n"); return 1; } int alpha[N], beta[N], gamma[N], delta[N]; double complex summand, sum; while(readPaulicoeffs(alpha, beta, gamma, delta, N)) { /* for(i=0; i