tests/agm.dat - master - niceSOFT/mpc

tests/agm.dat
47 строк5 KB
Andreas Enge
Rewrite AGM tests with a file and add generic tests.
07 апр 2022, 10:23
07 апр 2022, 10:232ab55ee
# Data file for mpc_agm
#
# Copyright (C) 2022 INRIA
#
# This file is part of GNU MPC.
#
# GNU MPC is free software; you can redistribute it and/or modify it under
# the terms of the GNU Lesser General Public License as published by the
# Free Software Foundation; either version 3 of the License, or (at your
#o ption) any later version.
#
# GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
# FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
# more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with this program. If not, see http://www.gnu.org/licenses/ .
#
# The line format respects the parameter order in the function prototype as
# follows:
#
# INEX_RE  INEX_IM  PREC_ROP_RE  ROP_RE  PREC_ROP_IM  ROP_IM  PREC_OP1_RE  OP1_RE  PREC_OP1_IM  OP1_IM  PREC_OP2_RE  OP2_RE  PREC_OP2_IM  OP2_IM  RND_RE  RND_IM
#
# See add.dat for more details.

# special cases of 0
0 0 53 0 53 0  53 -1 53 0  53 1 53 0 U D
0 0 53 0 53 0  53  0 53 0  53 1 53 1 U D
0 0 53 0 53 0  53  1 53 0  53 0 53 0 U D

# a=2^100 is real and b=-2^100+1 is close to -a, so that there is
# cancellation in the first step
+ + 1000 6.13343501065248469040220876057575735521422156659277284015370489065354785004941290892259011656852869912332814488649060620003153660129464591260581455777223198535415374247749106843365161158273510545819077989996849040655334209063638548567093807595920704639335161071221606819684075137748005838679614506936700e26 1000 2.78770347601954771576648332450140186390805588069947753721317890662882377009087432360024540561493332099222362713441719824215798470248682830506492225727606531299471177716509129100229573037698818619139314654444682343585665232750491711201076271867938956241318886378833024121607108490397434698394730099523728e28  200 1267650600228229401496703205376 53 0  200 -1267650600228229401496703205375 53 0 N N

# a=2^100 is real and b=2^100+I is close to a, which causes a loss in the
# error analysis when switching from a relative error to an absolute error
# (in ulp) for the imaginary part.
- + 1000 1.26765060022822940149670320537600000000000000000000000000000004930380657631323783823303533017413935457540219431393779814242309902707311989671728406142043830190405473381225698819654858454795973826194101850808953117443181105327393024632455534713154899279762797670267110875302148050274073504013709877973442e30 1000 4.99999999999999999999999999999999999999999999999999999999999980553077256683932165174799831937111748155461972746352709059016440615656003447359234530134921721757968337357997118603089092837261590319325008176491178329196315059130298011554852301129700508936432021390368960886831534563942526561512902959629403e-1  200 1267650600228229401496703205376 53 0 200 1267650600228229401496703205376 53 1 N N
# similar case for the real part when b is close to -a
+ + 1000 1.08467876592530886754246811029279356966981636533356136322144733028435886515931858750854237361072769893252665854960527081478235500456780213233344279303704706567669269542103227306403258182872546435486314875088697205523310058229087222631501223687780345556279062417814851470906262487893004445615660926562691e-2 1000 2.78905293895036173379978221446297250405753957415552085947732729471277330502161158283898862494994763003064685588768219963518185298464700472486122950691817734185737320652838329171863564571063685685782635224759876773001459257327047419614690899573835904762076822885023105374294261695197351950136206550502212e28  200 1267650600228229401496703205376 53 0 200 -1267650600228229401496703205376 53 1 N N

# Non-trivial example with an angle pi
- + 1000 1.73445281498669753232856248343970177809463913976908704246943226531366943181635177029314525135639219796336367868027485122390166753512240063464119645217760215308428521377134601470801777085207181315839747680069598739038938895527870204815473213761883712470377131465129627672808480329626971055399178573237530e0 1000 -7.02268113524478075006321405810247223230827765511578159358705685084648706076496079901548878047371750259623565862936612153628319914882054169762337559071422073790380764162955947410387987517205983880807663678020872942112760649464128028437895895107347428874323595807150644530530843707101528138966376736841785e0  20 1 20 2 20 -10 20 -20 N N

# Non-trivial example with an angle 0
- - 1000 4.25040709493227486172816431837313486679846786419019285967014766222375531274090378452528546078761717904588171358976686523664106901878258668543430057143043997188667013456002687950950378230536772481087956970495220412257232297324589475076978359364065270281502572385189827930845694706585008531069979410829173e0 1000 8.50081418986454972345632863674626973359693572838038571934029532444751062548180756905057092157523435809176342717953373047328213803756517337086860114286087994377334026912005375901900756461073544962175913940990440824514464594649178950153956718728130540563005144770379655861691389413170017062139958821658346e0  20 1 20 2 20 10 20 20 N N