slu_scomplex.h

Go to the documentation of this file.

 
#ifndef __SUPERLU_SCOMPLEX /* allow multiple inclusions */
#define __SUPERLU_SCOMPLEX
 
 
#ifndef SCOMPLEX_INCLUDE
#define SCOMPLEX_INCLUDE
 
typedef struct { float r, i; } singlecomplex;
 
#if defined(SUPERLU_TYPEDEF_COMPLEX) || DOXYGEN
typedef singlecomplex complex;
#endif
 
/* Macro definitions */
 
#define c_add(c, a, b) { (c)->r = (a)->r + (b)->r; \
                         (c)->i = (a)->i + (b)->i; }
 
#define c_sub(c, a, b) { (c)->r = (a)->r - (b)->r; \
                         (c)->i = (a)->i - (b)->i; }
 
#define cs_mult(c, a, b) { (c)->r = (a)->r * (b); \
                           (c)->i = (a)->i * (b); }
 
#define cc_mult(c, a, b) { \
        float 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 cc_conj(a, b) { \
        (a)->r = (b)->r; \
        (a)->i = -((b)->i); \
    }
 
#define c_eq(a, b)  ( (a)->r == (b)->r && (a)->i == (b)->i )
 
 
#ifdef __cplusplus
extern "C" {
#endif
 
/* Prototypes for functions in scomplex.c */
void c_div(singlecomplex *, singlecomplex *, singlecomplex *);
double c_abs(singlecomplex *);     /* exact */
double c_abs1(singlecomplex *);    /* approximate */
void c_exp(singlecomplex *, singlecomplex *);
void r_cnjg(singlecomplex *, singlecomplex *);
double r_imag(singlecomplex *);
singlecomplex c_sgn(singlecomplex *);
singlecomplex c_sqrt(singlecomplex *);
 
 
 
#ifdef __cplusplus
  }
#endif
 
#endif
 
#endif  /* __SUPERLU_SCOMPLEX */