Computational routine
eng

fftcmplx

File content

/* Modnumlib Scicos interfacing function
 * Copyright (C) 2009 Alan Layec
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 *
 * This library 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 library; if not, write to the
 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 * Boston, MA 02111-1307, USA.
 */
 
/* fftcmplx Scicos fft complexe
 * Type 4 simulation function ver 1.1b - scilab-3.0
 * 15 décembre 2004 - IRCOM GROUP - Author : A.Layec
 * 12 Mars 2009 - INRIA - Author : A.Layec
 */

/* REVISION HISTORY :
 * $Log$
 */

#include <scicos/scicos_block.h>
#include "modnum_lib.h"

/* Cette fonction réalise la fft complexe du vecteur {u1[];u2[]}.
 * u1 et u2 sont places dans y1 et y2 respectivement par dcopy puis la fft est réalisée.
 * Suivant la valeur de ipar[0] on réalise une fft directe (0 où -1) ou indirecte (1).
 * Suivant la valeur de ipar[1] on appele fft842(ipar[1]=0) où dfftmx(ipar[1]=1) via dfft2.
 * Si dfftmx est choisie alors on place la taille du mot de travail dans ipar[2].
 *
 * entrées réguliéres : u1[0..nu-1] : vecteur des réels du mot complexe d'entrée
 *                      u2[0..nu-1] : vecteur des imaginaires du mot complexe d'entrée
 *
 * sorties régulières : y1[0..nu-1] : vecteur des réels du mot complexe de sortie
 *                      y2[0..nu-1] : vecteur des imaginaires du mot complexe de sortie
 *
 * paramètres entiers : ipar[0] : signe de l'exponentiel (-1 où 0 : fft directe; 1 : fft indirecte)
 *                      ipar[1] : flag choix fft (0:fft842/1:dfftmx)
 *                      ipar[2] : taille du mot de travail pour dfftmx
 */

/*prototype*/
void fftcmplx(scicos_block *block,int flag)
{
 /*Déclaration des variables*/
 double *y1,*y2;
 double *u1,*u2;
 int nu,i,j,k,ierr;
 double *z__;

 /*fft variables*/
 fft_pr_struct fft_pr;

 /*Récupération des adresses des ports réguliers*/
 y1=(double *)block->outptr[0];
 y2=(double *)block->outptr[1];
 u1=(double *)block->inptr[0];
 u2=(double *)block->inptr[1];

 /*Récupération de la taille du port d'entrée*/
 nu=block->insz[0];

 /* init : fait l'allocation de work */
 if (flag==4) {
  if (block->ipar[1]==1) { /*cas fftmx*/
   /*allocation dynamique*/
   if ((*block->work=scicos_malloc(sizeof(double)*block->ipar[2]))== NULL) {
     set_block_error(-16);
     return;
   }
  }
 }
 /* calcule sortie */
 else if ((flag==1)||(flag==6)) {
  /* initialize fft_pr*/
  fft_pr.lfft   = nu;
  fft_pr.isn    = block->ipar[0];
  fft_pr.ffttyp = block->ipar[1];
  fft_pr.nwork  = block->ipar[2];
  fft_pr.work   = *block->work;

  /*Recopie u dans y*/
  cmplxcopy_c(&nu,&u1[0],&u2[0],&y1[0],&y2[0]);

  /*Appel routine fft*/
  ierr=calc_fft_c(&fft_pr,&y1[0],&y2[0]);
 }
 /*Terminaison */
 else if (flag==5) {
  /*Libère mémoire allouée*/
  if (block->ipar[1]==1) { /*cas fftmx*/
    scicos_free(*block->work);
  }
 }

}