dcomplex.h

Go to the documentation of this file.

 
/* 
 * This header file is to be included in source files z*.c
 */
#ifndef __SUPERLU_DCOMPLEX /* allow multiple inclusions */
#define __SUPERLU_DCOMPLEX
 
#include <mpi.h>
#include "superlu_defs.h"
 
typedef struct { double r, i; } doublecomplex;
 
/*
 * These variables will be defined to be MPI datatypes for complex
 * and double complex. I'm too lazy to declare them external in every
 * file that needs them.
 * 
 * Use WINOWS_EXPORT_ALL_SYMBOLS on windows to export all symbols when
 * building a shared library.
 * Introduce macro SUPERLU_DIST_EXPORT to correctly export the only
 * remaining data symbol SuperLU_MPI_DOUBLE_COMPLEX.
 */
//extern SUPERLU_DIST_EXPORT MPI_Datatype SuperLU_MPI_DOUBLE_COMPLEX;
 
/* Macro definitions */
 
#define z_copy(c, a) { (c)->r = (a)->r ; \
             (c)->i = (a)->i ; }
 
#define z_add(c, a, b) { (c)->r = (a)->r + (b)->r; \
             (c)->i = (a)->i + (b)->i; }
 
#define z_sub(c, a, b) { (c)->r = (a)->r - (b)->r; \
             (c)->i = (a)->i - (b)->i; }
 
#define zd_mult(c, a, b) { (c)->r = (a)->r * (b); \
                           (c)->i = (a)->i * (b); }
 
#define zz_mult(c, a, b) { \
    double cr, ci; \
        cr = (a)->r * (b)->r - (a)->i * (b)->i; \
        ci = (a)->i * (b)->r + (a)->r * (b)->i; \
        (c)->r = cr; \
        (c)->i = ci; \
    }
 
#define z_eq(a, b)  ( (a)->r == (b)->r && (a)->i == (b)->i )
 
 
#ifdef __cplusplus
extern "C" {
#endif
 
/* Prototypes for functions in dcomplex.c */
void   slud_z_div(doublecomplex *, doublecomplex *, doublecomplex *);
double slud_z_abs(doublecomplex *);     /* exact */
double slud_z_abs1(doublecomplex *);    /* approximate */
 
 
#ifdef __cplusplus
  }
#endif
 
 
#endif  /* __SUPERLU_DCOMPLEX */