DFT matrix generator (gcc -std=c11 -lm dftmgen.c)

dftmgen.c

#include <complex.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>

#define M_PI 3.1415926535897932384626433832795L

int main(int argc, char** argv) {
	int n = 0;
	if(argc > 1) n = atoi(argv[1]);
	for(int i = 0; i < n; i++) {
		for(int j = 0; j < n; j++) {
			complex double omega = cos(2*M_PI/n) - sin(2*M_PI/n)*I;
			omega = cpow(omega, i*j);
			printf("(%f,%f)", creal(omega), cimag(omega));
		}
		printf("\n");
	}
}

Line 13 should read cos(2*M_PI/n) - sin(2*M_PI/n)*I, as per the trigonometric identity for e^(-ix). Otherwise, thanks for the snippet!