source: buchla-68k/orig/BUCHLA/FRQPROM.C

/*
   =============================================================================
        frqprom.c -- generate the frequency / pitch PROM for the Buchla 700
        Written by:  D.N. Lynx Crowe -- see VER below for version
   =============================================================================
*/

#define VER     "Version 12 - 1987-09-17"       /* version ID string */

#define DEBUGIT         0

#include "stdio.h"
#include "math.h"
#include "stddefs.h"

#define N       (double)192.0           /* number of steps in the oscillator */
#define R       (double)524288.0        /* number of steps in the phase angle */
#define C       (double)10066000.0      /* clock frequency */

#define RNC     ((R*N)/C)               /* frequency multiplier */

#if DEBUGIT
#define FILEOUT         "FRQPROM.PRN"
#else
#define FILEOUT         "FRQPROM.MOT"
#endif

#define NOTIFY1 10      /* frquency notification interval */
#define NOTIFY2 600     /* pitch notification interval */

#define HILIMIT (double)16383.5

extern  int     msrec();

double  p, q, s, t, inc, f;

unsigned short tab[32768];
unsigned char otab[65536];

int count, notify;

FILE *fp;

/* 

*/

main()
{
        register long i;

        printf("Buchla 700 Frequency / Pitch PROM Generator %s\n\n", VER);

        q = RNC;

        if (NULL EQ (fp = fopen(FILEOUT, "wa"))) {      /* open output file */

                printf("ERROR:  Unable to open [%s] for output.\n", FILEOUT);
                exit(1);
        }

        printf("Output on [%s]\n", FILEOUT);
        printf("Using:  R = %f, N = %f\n", R, N);
        printf("        C = %f, (R*N)/C = %f\n\n", C, q);

#if DEBUGIT
        fprintf(fp, "Buchla 700 Frequency / Pitch PROM Generator %s\n\n", VER);
        fprintf(fp, "Using:  R = %f, N = %f\n", R, N);
        fprintf(fp, "        C = %f, (R*N)/C = %f\n\n", C, q);
        fprintf(fp, "Frequency portion of PROM\n");
#endif

        for (i = 0; i < 65536; i++)     /* preset table to saturation value */
                tab[i] = 65535;

        /* generate the frequency portion of the table from 0.0 to 15.95 Hz */
        /* using the formula:  inc = f * ((R*N)/C)  */

        printf("Generating frequency portion of PROM file\n");
        notify = NOTIFY1;
        count = notify;

#if DEBUGIT
        for (f = (double)0.0; f < (double)16.0; f += ((double)0.05 * (double)NOTIFY1)) {
#else
        for (f = (double)0.0; f < (double)16.0; f += (double)0.05) {
#endif
                inc = f * q;
                i = (long)floor(f * 20.0);
                tab[i] = floor(inc);

#if DEBUGIT
                fprintf(fp, "%5d 0x%04x %14.4f %14.4f\n",
                        i, tab[i], f, inc);
#endif

                if (++count GE notify) {

                        count = 0;
                        printf(".");
                }
        }

        printf("\n");

#if DEBUGIT
        fprintf(fp, "\n\nPitch portion of PROM\n");
#endif
/* 

*/
        /* generate the pitch portion of the table from 160.0 to HILIMIT cents */
        /* according to the formula:  inc = ((R*K*N)/C)*(2**(p/1200)) */

        printf("Generating pitch portion of PROM file\n");
        notify = NOTIFY2;
        count = notify;

#if DEBUGIT
        for (p = (double)160.0; p < HILIMIT; p += ((double)0.5 * (double)NOTIFY2)) {
#else
        for (p = (double)160.0; p < HILIMIT; p += (double)0.5) {
#endif

                s = p / (double)1200.0;         /* calculate the increment */
                t = pow((double)2.0, s);
                inc = (double)149.082 * t;

                if (inc GE 65536.0)             /* saturate */
                        break;

                i = (long)floor((double)2.0 * p);       /* stuff the table */
                tab[i] = floor(inc);

#if DEBUGIT
                fprintf(fp, "%5d 0x%04x %14.4f %14.4f %14.4f %14.4f\n",
                        i, tab[i], p, s, t, inc);
#endif

                if (++count GE notify) {

                        count = 0;
                        printf(".");
                }
        }

        printf("\n");

/* 

*/
        /* output the table */

#if !DEBUGIT

        printf("Splitting table into LS and MS bytes for output\n");

        for (i = 0; i < 32768; i++) {   /* split into high and low bytes */

                otab[i] = tab[i] & 0x00FF;      /* LS byte in low half */
                otab[i + 32768L] = tab[i] >> 8; /* MS byte in high half */
        }

        printf("Writing Motorola S-Records to [%s]\n", FILEOUT);
        msrec(fp, 0L, 65536L, otab);    /* output the table as S-Records */

#endif

        fclose(fp);                     /* close the file */
        printf("Program complete.\n");  /* say we're done */
        exit(0);                        /* go away */
}