File: polynomials.tst

package info (click to toggle)
gap-float 0.9.1%2Bds-6
  • links: PTS, VCS
  • area: main
  • in suites: bookworm, bullseye
  • size: 612 kB
  • sloc: ansic: 2,537; cpp: 1,998; xml: 184; makefile: 103; sh: 1
file content (30 lines) | stat: -rw-r--r-- 970 bytes parent folder | download | duplicates (6)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
#############################################################################
##
#W  polynomials.tst              Float Package              Laurent Bartholdi
##
#Y  Copyright (C) 2013,  Laurent Bartholdi
##
#############################################################################
##
##  This file tests polynomials
##
#############################################################################

gap> START_TEST("polynomials");
gap> x := Indeterminate(field,"x");
x
gap> x+1;; x+1.0;; 1+x;; 1-x;; 1*x;; 1.0*x;; x*1.0;; x*x;
x^2
gap> Value(x^2+1,2)=5.0;
true
gap> roots := [0.1,0.1,0.1,0.5,0.6,0.7];;
gap> poly := Product(x-roots);;
gap> newroots := RootsFloat(poly);;
gap> Filtered([1..3],i->AbsoluteValue(roots[i]-newroots[i])>1.e-6);
[  ]
gap> Filtered([4..Length(roots)],i->AbsoluteValue(roots[i]-newroots[i])>1.e-14);
[  ]
gap> STOP_TEST( "polynomials.tst", 3*10^8 );
polynomials

#E polynomials.tst . . . . . . . . . . . . . . . . . . . . . . . . .ends here