--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/spandsp-0.0.3/spandsp-0.0.3/tests/dds_tests.c	Fri Jun 25 16:00:21 2010 +0200
@@ -0,0 +1,257 @@
+/*
+ * SpanDSP - a series of DSP components for telephony
+ *
+ * dds_tests.c
+ *
+ * Written by Steve Underwood <steveu@coppice.org>
+ *
+ * Copyright (C) 2003 Steve Underwood
+ *
+ * All rights reserved.
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License version 2, as
+ * published by the Free Software Foundation.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+ *
+ * $Id: dds_tests.c,v 1.17 2006/11/19 14:07:26 steveu Exp $
+ */
+
+/*! \file */
+
+/*! \page dds_tests_page Direct digital synthesis tests
+\section dds_tests_page_sec_1 What does it do?
+???.
+
+\section dds_tests_page_sec_2 How does it work?
+???.
+*/
+
+#ifdef HAVE_CONFIG_H
+#include "config.h"
+#endif
+
+#include <inttypes.h>
+#include <stdlib.h>
+#include <stdio.h>
+#include <fcntl.h>
+#include <memory.h>
+#if defined(HAVE_TGMATH_H)
+#include <tgmath.h>
+#endif
+#if defined(HAVE_MATH_H)
+#include <math.h>
+#endif
+#include <audiofile.h>
+#include <tiffio.h>
+
+#include "spandsp.h"
+
+#define OUTPUT_FILE_NAME            "dds.wav"
+#define OUTPUT_FILE_NAME_COMPLEX    "complex_dds.wav"
+
+#define SAMPLES_PER_CHUNK           8000
+
+int main(int argc, char *argv[])
+{
+    int i;
+    uint32_t phase;
+    int32_t phase_inc;
+    int outframes;
+    complexf_t camp;
+    int16_t buf[2*SAMPLES_PER_CHUNK];
+    AFfilehandle outhandle;
+    AFfilesetup filesetup;
+    power_meter_t meter;
+    power_meter_t meter_i;
+    power_meter_t meter_q;
+
+    power_meter_init(&meter, 10);
+
+    printf("Non-complex DDS tests.\n");
+    if ((filesetup = afNewFileSetup()) == AF_NULL_FILESETUP)
+    {
+        fprintf(stderr, "    Failed to create file setup\n");
+        exit(2);
+    }
+    afInitSampleFormat(filesetup, AF_DEFAULT_TRACK, AF_SAMPFMT_TWOSCOMP, 16);
+    afInitRate(filesetup, AF_DEFAULT_TRACK, (float) SAMPLE_RATE);
+    afInitFileFormat(filesetup, AF_FILE_WAVE);
+    afInitChannels(filesetup, AF_DEFAULT_TRACK, 1);
+
+    if ((outhandle = afOpenFile(OUTPUT_FILE_NAME, "w", filesetup)) == AF_NULL_FILEHANDLE)
+    {
+        fprintf(stderr, "    Cannot create wave file '%s'\n", OUTPUT_FILE_NAME);
+        exit(2);
+    }
+
+    phase = 0;
+
+    printf("Test with 123.456789Hz.\n");
+    phase_inc = dds_phase_rate(123.456789);
+    for (i = 0;  i < SAMPLES_PER_CHUNK;  i++)
+    {
+        buf[i] = alaw_to_linear(linear_to_alaw(dds(&phase, phase_inc)));
+        power_meter_update(&meter, buf[i]);
+    }
+    outframes = afWriteFrames(outhandle,
+                              AF_DEFAULT_TRACK,
+                              buf,
+                              SAMPLES_PER_CHUNK);
+    if (outframes != SAMPLES_PER_CHUNK)
+    {
+        fprintf(stderr, "    Error writing wave file\n");
+        exit(2);
+    }
+    printf("Level is %fdBOv\n", power_meter_dbov(&meter));
+    if (fabs(power_meter_dbov(&meter) + 3.02) > 0.05)
+    {
+        printf("Test failed.\n");
+        exit(2);
+    }
+
+    printf("Test with 12.3456789Hz.\n");
+    phase_inc = dds_phase_rate(12.3456789);
+    for (i = 0;  i < SAMPLES_PER_CHUNK;  i++)
+    {
+        buf[i] = alaw_to_linear(linear_to_alaw(dds(&phase, phase_inc)));
+        power_meter_update(&meter, buf[i]);
+    }
+    outframes = afWriteFrames(outhandle,
+                              AF_DEFAULT_TRACK,
+                              buf,
+                              SAMPLES_PER_CHUNK);
+    if (outframes != SAMPLES_PER_CHUNK)
+    {
+        fprintf(stderr, "    Error writing wave file\n");
+        exit(2);
+    }
+    printf("Level is %fdBOv\n", power_meter_dbov(&meter));
+    /* Use a wider tolerance for this very low frequency - the power meter will ripple */
+    if (fabs(power_meter_dbov(&meter) + 3.02) > 0.2)
+    {
+        printf("Test failed.\n");
+        exit(2);
+    }
+
+    printf("Test with 2345.6789Hz.\n");
+    phase_inc = dds_phase_rate(2345.6789);
+    for (i = 0;  i < SAMPLES_PER_CHUNK;  i++)
+    {
+        buf[i] = alaw_to_linear(linear_to_alaw(dds(&phase, phase_inc)));
+        power_meter_update(&meter, buf[i]);
+    }
+    outframes = afWriteFrames(outhandle,
+                              AF_DEFAULT_TRACK,
+                              buf,
+                              SAMPLES_PER_CHUNK);
+    if (outframes != SAMPLES_PER_CHUNK)
+    {
+        fprintf(stderr, "    Error writing wave file\n");
+        exit(2);
+    }
+    printf("Level is %fdBOv\n", power_meter_dbov(&meter));
+    if (fabs(power_meter_dbov(&meter) + 3.02) > 0.05)
+    {
+        printf("Test failed.\n");
+        exit(2);
+    }
+
+    printf("Test with 3456.789Hz.\n");
+    phase_inc = dds_phase_rate(3456.789);
+    for (i = 0;  i < SAMPLES_PER_CHUNK;  i++)
+    {
+        buf[i] = alaw_to_linear(linear_to_alaw(dds(&phase, phase_inc)));
+        power_meter_update(&meter, buf[i]);
+    }
+    outframes = afWriteFrames(outhandle,
+                              AF_DEFAULT_TRACK,
+                              buf,
+                              SAMPLES_PER_CHUNK);
+    if (outframes != SAMPLES_PER_CHUNK)
+    {
+        fprintf(stderr, "    Error writing wave file\n");
+        exit(2);
+    }
+    printf("Level is %fdBOv\n", power_meter_dbov(&meter));
+    if (fabs(power_meter_dbov(&meter) + 3.02) > 0.05)
+    {
+        printf("Test failed.\n");
+        exit(2);
+    }
+
+    if (afCloseFile(outhandle) != 0)
+    {
+        fprintf(stderr, "    Cannot close wave file '%s'\n", OUTPUT_FILE_NAME);
+        exit(2);
+    }
+    afFreeFileSetup(filesetup);
+
+    printf("Complex DDS tests,\n");
+    if ((filesetup = afNewFileSetup()) == AF_NULL_FILESETUP)
+    {
+        fprintf(stderr, "    Failed to create file setup\n");
+        exit(2);
+    }
+    afInitSampleFormat(filesetup, AF_DEFAULT_TRACK, AF_SAMPFMT_TWOSCOMP, 16);
+    afInitRate(filesetup, AF_DEFAULT_TRACK, (float) SAMPLE_RATE);
+    afInitFileFormat(filesetup, AF_FILE_WAVE);
+    afInitChannels(filesetup, AF_DEFAULT_TRACK, 2);
+
+    if ((outhandle = afOpenFile(OUTPUT_FILE_NAME_COMPLEX, "w", filesetup)) == AF_NULL_FILEHANDLE)
+    {
+        fprintf(stderr, "    Cannot create wave file '%s'\n", OUTPUT_FILE_NAME_COMPLEX);
+        exit(2);
+    }
+
+    power_meter_init(&meter_i, 7);
+    power_meter_init(&meter_q, 7);
+
+    phase = 0;
+    phase_inc = dds_phase_ratef(123.456789);
+    for (i = 0;  i < SAMPLES_PER_CHUNK;  i++)
+    {
+        camp = dds_complexf(&phase, phase_inc);
+        buf[2*i] = camp.re*10000.0;
+        buf[2*i + 1] = camp.im*10000.0;
+        power_meter_update(&meter_i, buf[2*i]);
+        power_meter_update(&meter_q, buf[2*i]);
+    }
+    outframes = afWriteFrames(outhandle,
+                              AF_DEFAULT_TRACK,
+                              buf,
+                              SAMPLES_PER_CHUNK);
+    if (outframes != SAMPLES_PER_CHUNK)
+    {
+        fprintf(stderr, "    Error writing wave file\n");
+        exit(2);
+    }
+    printf("Level is %fdBOv/%fdBOv\n", power_meter_dbov(&meter_i), power_meter_dbov(&meter_q));
+    if (fabs(power_meter_dbov(&meter_i) + 13.42) > 0.05
+        ||
+        fabs(power_meter_dbov(&meter_q) + 13.42) > 0.05)
+    {
+        printf("Test failed.\n");
+        exit(2);
+    }
+
+    if (afCloseFile(outhandle) != 0)
+    {
+        fprintf(stderr, "    Cannot close wave file '%s'\n", OUTPUT_FILE_NAME_COMPLEX);
+        exit(2);
+    }
+    afFreeFileSetup(filesetup);
+
+    printf("Tests passed.\n");
+    return  0;
+}
+/*- End of function --------------------------------------------------------*/
+/*- End of file ------------------------------------------------------------*/