bytemedia

kagome_kagome_clarinet.c (3890B)

//#define _POSIX_C_SOURCE 200809L
#define _XOPEN_SOURCE 700
#include <stdio.h>
#include <inttypes.h> // uint32_t
#include <math.h> // sinf(), M_PI, pow()
#include <stdlib.h> // abort()
#include "au.h"
// https://pages.mtu.edu/~suits/NoteFreqCalcs.html
#define TWELTH_ROOT_2 1.059463094359
double
note_to_hertz(double ref_tune, int steps) {
	return ref_tune*pow(TWELTH_ROOT_2,steps);
}
void
sin_note(int dur, struct au_header *header, int hertz, float volume, float delay) {
	char buf[4];
	for(int t = 0;t<dur;t++) {
		//putchar((uint8_t)(sinf((delay + t) * 2*M_PI*hertz/header->samplerate) * volume) % 210);
		float w1 = 2*M_PI*hertz/header->samplerate;
		// https://pages.mtu.edu/~suits/clarinet.html
		float tune = sin(w1*t) + 0.75*sin(3*w1*t) + 0.5*sin(5*w1*t) + 0.14*sin(7*w1*t) + 0.5*sin(9*w1*t) + 0.12*sin(11*w1*t) + 0.17*sin(13*w1*t);
		putchar((uint8_t)(tune * volume));
	}
}
#if 0
void
note(int dur, int pitch) {
	for(int t = 0;t<dur;t++) {
		putchar( t*pitch );
	}
}
#else
#define REF_TUNE 440/4
#define VOLUME 60
#define note(dur,pitch) sin_note(dur, header, note_to_hertz(REF_TUNE, pitch-8), VOLUME, 0)
#endif
void
silence(int dur) {
	for(int t = 0;t<dur;t++)putchar( 0x00 );
}
void
kagome_kagome(struct au_header *header) {
	int tick = header->samplerate / 1.9;
	//int note = tick * 0.99;
	//int gap  = tick * 0.01;
	int note = tick;
	int gap = 0;
	note(note*2, 8); silence(gap*2);
	note(note*1, 8); silence(gap);
	note(note*1, 9); silence(gap);
	// ----
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	silence(tick*1);
	// ----
	note(note*1, 8); silence(gap);
	note(note/2, 8); silence(gap/2);
	note(note/2, 8); silence(gap/2);
	note(note*1, 8); silence(gap);
	note(note/2, 7); silence(gap/2);
	note(note/2, 7); silence(gap/2);
	// ----
	note(note*1, 8); silence(gap);
	note(note/2, 8); silence(gap/2);
	note(note/2, 7); silence(gap/2);
	note(note*1, 5); silence(gap);
	silence(tick*1);
	// ----
	note(note, 8); silence(gap);
	note(note, 7); silence(gap);
	note(note, 8); silence(gap);
	note(note, 7); silence(gap);
	// ----
	note(note*1, 8); silence(gap);
	note(note/2, 8); silence(gap/2);
	note(note/2, 7); silence(gap/2);
	note(note*1, 5); silence(gap);
	silence(tick*1);
	// ----
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	note(note*1, 9); silence(gap);
	// ----
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	silence(tick*1);
	// ----
	note(note*1, 8); silence(gap);
	note(note/2, 7); silence(gap/2);
	note(note/2, 7); silence(gap/2);
	note(note*1, 8); silence(gap);
	note(note/2, 7); silence(gap/2);
	note(note/2, 7); silence(gap/2);
	// ----
	note(note*1, 8); silence(gap);
	note(note*1, 8); silence(gap);
	note(note*1, 5); silence(gap);
	silence(tick*1);
	// ----
	note(note/2, 8); silence(gap/2);
	note(note/2, 8); silence(gap/2);
	note(note/2, 8); silence(gap/2);
	note(note/2, 8); silence(gap/2);
	note(note*1, 8); silence(gap);
	note(note*1, 9); silence(gap);
	// ----
	note(tick*1, 8); // silence(gap);
	note(note*1, 7); silence(gap);
	note(note*1, 8); silence(gap);
	silence(tick*1);
}
int
main(void)
{
	// Use AU defaults
	struct au_header header = {
		.offset = 24, // no-annotation, in octets
		.length = 0xFFFFFFFF, // unknown data size
		//.encoding = AU_ENCODING_32B_FPCM,
		.encoding = 02,
		.samplerate = 48000, // Hz
		.channels = 1
	};
	// fd 1 is stdout
	write_au_header(1, &header);
	kagome_kagome(&header);
	//       dur    pt  vol del
	//sin_note(header.samplerate, &header, 440, 150, 0);
	//sin_note(header.samplerate/2, &header, note_to_hertz(440, 0), 150, 0);
	//sin_note(header.samplerate/4, &header, note_to_hertz(440, 1), 150, 0);
	//sin_note(header.samplerate/4, &header, note_to_hertz(440, 1), 150, 0);
	return 0;
}