//
// Listing 1.1: The direct realization of the Discrete Sine Transform (DST):
//

#define M_PI 3.14159265358979323846

long bin,k;
double arg;
for (bin = 0; bin < transformLength; bin++) {
	
	transformData[bin] = 0.;
	for (k = 0; k < transformLength; k++) {
		
		arg = (float)bin * M_PI *(float)k / (float)transformLength;
		transformData[bin] += inputData[k] * sin(arg);
		
	}
	
}