complex.cc File Reference

Detailed Description

Implementation of complex number functions.

Copyright

• (c) 2002 The University of Chicago, as Operator of Argonne National Laboratory.
• (c) 2002 The Regents of the University of California, as Operator of Los Alamos National Laboratory.

License

This file is distributed under the terms of the Software License Agreement found in the file LICENSE included with this distribution.

Author

M. Borland, R. Soliday

Definition in file complex.cc.

#include <complex>
#include "mdb.h"

Go to the source code of this file.

Functions
std::complex< double >	complexErf (std::complex< double > z, long *flag)
	Computes the complex error function of a given complex number.
std::complex< double >	cexpi (double p)
	Computes the complex exponential of an imaginary number.
std::complex< double >	cipowr (std::complex< double > a, int n)
	Raises a complex number to an integer power.
void	complex_multiply (double *r0, double *i0, double r1, double i1, double r2, double i2)
	Multiplies two complex numbers.
void	complex_divide (double *r0, double *i0, double r1, double i1, double r2, double i2, double threshold)
	Divides two complex numbers.

Function Documentation

cexpi()

std::complex< double > cexpi

(

double

p

)

Computes the complex exponential of an imaginary number.

Parameters

p	The real number representing the imaginary part.

Returns

The complex exponential of the input.

Definition at line 43 of file complex.cc.

                                   {
  std::complex<double> a;
  a = cos(p) + std::complex<double>(0, 1) * sin(p);
  return (a);
}

cipowr()

std::complex< double > cipowr	(	std::complex< double >	a,
		int	n )

Raises a complex number to an integer power.

Parameters

a	The base complex number.
n	The exponent (integer).

Returns

The complex number a raised to the power n.

Definition at line 56 of file complex.cc.

                                                       {
  int i;
  std::complex<double> p(1, 0);
 
  if (n >= 0) {
    for (i = 0; i < n; i++)
      p = p * a;
    return (p);
  }
  a = p / a;
  return (cipowr(a, -n));
}

complex_divide()

void complex_divide	(	double *	r0,
		double *	i0,
		double	r1,
		double	i1,
		double	r2,
		double	i2,
		double	threshold )

Divides two complex numbers.

Deprecated

These routines are obsolete, really, but some code uses them.

Parameters

r0	Pointer to store the real part of the result.
i0	Pointer to store the imaginary part of the result.
r1	Real part of the numerator complex number.
i1	Imaginary part of the numerator complex number.
r2	Real part of the denominator complex number.
i2	Imaginary part of the denominator complex number.
threshold	The threshold to prevent division by very small numbers.

Definition at line 105 of file complex.cc.

                    {
  double tempr, denom;
 
  if ((denom = sqr(r2) + sqr(i2)) < threshold)
    denom = threshold;
  i2 = -i2;
  tempr = (r1 * r2 - i1 * i2) / denom;
  *i0 = (r1 * i2 + i1 * r2) / denom;
  *r0 = tempr;
}

complex_multiply()

void complex_multiply	(	double *	r0,
		double *	i0,
		double	r1,
		double	i1,
		double	r2,
		double	i2 )

Multiplies two complex numbers.

Deprecated

These routines are obsolete, really, but some code uses them.

Parameters

r0	Pointer to store the real part of the result.
i0	Pointer to store the imaginary part of the result.
r1	Real part of the first complex number.
i1	Imaginary part of the first complex number.
r2	Real part of the second complex number.
i2	Imaginary part of the second complex number.

Definition at line 81 of file complex.cc.

                        {
  double tempr;
 
  tempr = r1 * r2 - i1 * i2;
  *i0 = r1 * i2 + i1 * r2;
  *r0 = tempr;
}

complexErf()

std::complex< double > complexErf	(	std::complex< double >	z,
		long *	flag )

Computes the complex error function of a given complex number.

Parameters

z	The complex number input.
flag	Pointer to a flag variable to store computation status.

Returns

The complex error function value of z.

Definition at line 26 of file complex.cc.

                                                                {
  double xi, yi;
  double u = 0, v = 0;
  long lflag = 0;
  xi = z.real();
  yi = z.imag();
  wofz(&xi, &yi, &u, &v, &lflag);
  *flag = lflag;
  return std::complex<double>(u, v);
}