Routine bas-niveau
fr - eng


cmplxop_c - routines de calculs pour des nombres complexes

Module

Contenu du fichier


/* cmplxadd_c subroutine
 * complex addition computation
 *
 * Copyright (C) 2007-2011 Alan Layec
 *
 * This file is part of modnumlib.
 *
 * modnumlib 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.
 *
 * modnumlib 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 modnumlib; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 *
 */

/* REVISION HISTORY :
 * $Log$
 */

#include "modnum_lib.h"

/* cmplxadd_c routine de calcul d'addition-soustraction de vecteurs complexes
 *
 * n           : la taille des vecteurs
 * sig         : signe de l'opération
 * [z1_r;z1_i] : adresses de départ du vecteur complexe 1
 * [z2_r;z2_i] : adresses de départ du vecteur complexe 2
 * [y_r;y_i]   : adresses de départ du vecteur complexe résultat
 *
 * rmq : doit exister en version BLAS(!?)
 */

void cmplxadd_c(int *n,int *sig,double *z1_r,double *z1_i,double *z2_r,double *z2_i,double *y_r,double *y_i)
{
 /*déclaration des variables*/
 int i;

 /*réalise multiplication complex*/
 for(i=0;i<(*n);i++) {
  y_r[i]=z1_r[i]+(*sig)*z2_r[i];
  y_i[i]=z1_i[i]+(*sig)*z2_i[i];
 }
 return;
}

/* cmplxcopy_c routine qui copie deux tableaux dans deux autres
 * n     : longueur du tableau
 * u1    : tableau en entrée (double)
 * u2    : tableau en entrée (double)
 * y1    : tableau en sortie (double)
 * y2    : tableau en sortie (double)
 */

void cmplxcopy_c(int *n,double *u1,double *u2,double *y1,double *y2)
{
 /*Déclaration des variables*/
 int i;

 for(i=0;i<(*n);i++) {
   y1[i] = u1[i];
 }
 for(i=0;i<(*n);i++) {
   y2[i] = u2[i];
 }
 return;
}

/*
 * cmplxmult_c routine de calcul de multiplication de vecteurs complexes
 *
 * n           : la taille des vecteurs
 * [z1_r;z1_i] : adresses de départ du vecteur complexe 1
 * [z2_r;z2_i] : adresses de départ du vecteur complexe 2
 * [y_r;y_i]   : adresses de départ du vecteur complexe résultat
 *
 * rmq : doit exister en version BLAS(!?)
 */

void cmplxmult_c(int *n,double *z1_r,double *z1_i,double *z2_r,double *z2_i,double *y_r,double *y_i)
{
 /*déclaration des variables*/
 int i;

 /*réalise multiplication complex*/ 
 /*F2C(wmmul)(&z1_r[0],&z1_i[0],(n1=1,&n1),&z2_r[0],&z2_i[0],(m=1,&m),&y_r[0],&y_i[0],(m=1,&m),n,(m=1,&m),n);*/
 for(i=0;i<(*n);i++) {
  y_r[i]=z1_r[i]*z2_r[i]-z1_i[i]*z2_i[i];
  y_i[i]=z1_r[i]*z2_i[i]+z1_i[i]*z2_r[i];
 }
 return;
}

Auteurs

A. Layec A. Layec