#ifndef dm_spectrum_H
#define dm_spectrum_H

/* Filtered DNA sequences */
/* Last edited on 2008-01-29 15:00:18 by hcgl */

#define dm_spectrum_H_COPYRIGHT \
  "Copyright © 2005  by the State University of Campinas (UNICAMP)" \
  " and the Federal Fluminense University (UFF)"

#include <dm_basic.h>
#include <dm_sample.h>
#include <dm_datum.h>
#include <dm_seq.h>

#include <msm_ps_tools.h>

#include <vec.h>
#include <pswr.h>

/* This interface provides tools to compute and plot
  the power spectrum of a sampled multi-channel signal
  (a {dm-seq_t}). */

double_vec_t dm_spectrum_from_seq(dm_seq_t *seqp);
  /* Computes the power spectrum {pwr[0..fMax]} of sequence {seqp},
    where {fMax = nd/2} is the max absolute frequency, and {nd} is the
    number of datums in the sequence. The the three channels are
    decoded and scaled by {seqp->sfac}, and their power spectra are
    added together. */

/* POSTSCRIPT PLOTTING */

void dm_spectrum_postscript_plot(msm_ps_tools_t *dp, dm_seq_t *seqp);
  /* Generates a Postscript plot of of the power spectrum of sequence
    {seqp}, written out to the Postscript plotting stream {dp}. The
    power spectra of all three channels are added together. The plot
    uses the whole plotting area of {dp} (with a bit of slack). */ 

void dm_spectrum_postscript_plot_named
  ( dm_seq_t *seqp, 
    double hSize, 
    double vSize,
    double mrg,
    char *name,
    char *tag
  );
  /* Generates a Postscript plot of the power spectrum of sequence
    {seqp}, written out to the Encapsulated Postscript file called
    "{name}{tag}.eps".  The power spectra of all three channels are added
    together.  The figure will have {hSize} by {vSize}
    millimeters, including a blank margin {mrg} millimeters wide all
    around.  The plot uses the whole plotting area of {dp} (with a bit
    of slack). */

#endif