Quelle cxsc_poly.C Sprache: C

/****************************************************************************
*
* cxsc_poly.C                                              Laurent Bartholdi
*
* Copyright (c) 2011, Laurent Bartholdi
*
****************************************************************************
*
* driver for cpoly.C, using cxsc complex numbers
*
****************************************************************************/

#include <real.hpp>
#include <complex.hpp>

#include "cxsc_poly.h"

#define xMAX_EXP DBL_MAX_EXP
#define xMIN_EXP DBL_MIN_EXP
#define xINFIN DBL_MAX

typedef double xreal;
#define xcomplex cxsc::complex

static const xreal xabs(const xcomplex &newz)
{
  return _double(abs(newz));
}

static const xreal xabs(const xreal &newz){
  return fabs(newz);
}

static const xreal xnorm(const xcomplex &newz)
{
  return _double(abs2(newz));
}

static const xreal xroot(xreal x, int n)
{
  return pow(x,1.0/n);
}

static const xreal xsqrt(xreal x)
{
  return sqrt(x);
}

static const int xbits(const xcomplex &z)
{
  return 53;
}

static const xreal xeta(const xcomplex &z)
{
  return _double(cxsc::Epsilon);
}

static const long int xlogb(const xcomplex &newz)
{
  long e0 = expo(Re(newz)), e1 = expo(Im(newz));
  return e0 > e1 ? e0 : e1;
}

static void xscalbln(xcomplex *z, long int a)
{
  cxsc::real re = Re(*z), im = Im(*z);

  times2pown(re,a);
  times2pown(im,a);
  *z = cxsc::complex(re,im);
}

int default_prec;

#define cpoly cpoly_CXSC
#include "cpoly.C"

quality92%

¤ Dauer der Verarbeitung: 0.12 Sekunden (vorverarbeitet) ¤

Wurzel

Suchen

Beweissystem der NASA

Beweissystem Isabelle

NIST Cobol Testsuite

Cephes Mathematical Library

Wiener Entwicklungsmethode

Haftungshinweis

Die Informationen auf dieser Webseite wurden nach bestem Wissen sorgfältig zusammengestellt. Es wird jedoch weder Vollständigkeit, noch Richtigkeit, noch Qualität der bereit gestellten Informationen zugesichert.

Bemerkung:

Die farbliche Syntaxdarstellung ist noch experimentell.