| Quellcodebibliothek Statistik Leitseite products/Sources/formale Sprachen/GAP/pkg/float/src/ (Algebra von RWTH Aachen Version 4.15.1©) Datei vom 26.7.2025 mit Größe 1 kB |
|
Quelle c_poly.C Sprache: C |
|||||||||
|
/****************************************************************************
* * c_poly.C Laurent Bartholdi * * Copyright (c) 2011, Laurent Bartholdi * **************************************************************************** * * driver for cpoly.C, using standard C++ complex numbers * * This code is currently unused, but remains there to test the root finding * method outside of GAP. * * It would be revived if someone implemented double-precision complex numbers * in float (currently there are only multi-precision complex numbers). * ****************************************************************************/ #include <real.hpp> #include <complex.hpp> #define cpoly cpoly_C #define xMAX_EXP DBL_MAX_EXP #define xMIN_EXP DBL_MIN_EXP #define xINFIN DBL_MAX typedef double xreal; typedef _complex double xcomplex; static const xreal xabs(const xcomplex &newz) { return _double(abs(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 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; #include "cpoly.C" ¤ Dauer der Verarbeitung: 0.13 Sekunden (vorverarbeitet) ¤*© Formatika GbR, Deutschland |
|
2026-03-28