2015-03-20 14:35:02 -06:00
|
|
|
// Copyright 2015 The Go Authors. All rights reserved.
|
|
|
|
// Use of this source code is governed by a BSD-style
|
|
|
|
// license that can be found in the LICENSE file.
|
|
|
|
|
|
|
|
package big_test
|
|
|
|
|
|
|
|
import (
|
2015-05-22 16:21:56 -06:00
|
|
|
"cmd/compile/internal/big"
|
2015-03-20 14:35:02 -06:00
|
|
|
"fmt"
|
|
|
|
"math"
|
|
|
|
)
|
|
|
|
|
|
|
|
func ExampleFloat_Add() {
|
|
|
|
// Operating on numbers of different precision.
|
|
|
|
var x, y, z big.Float
|
|
|
|
x.SetInt64(1000) // x is automatically set to 64bit precision
|
|
|
|
y.SetFloat64(2.718281828) // y is automatically set to 53bit precision
|
|
|
|
z.SetPrec(32)
|
|
|
|
z.Add(&x, &y)
|
2015-06-03 16:11:05 -06:00
|
|
|
fmt.Printf("x = %.10g (%s, prec = %d, acc = %s)\n", &x, x.Text('p', 0), x.Prec(), x.Acc())
|
|
|
|
fmt.Printf("y = %.10g (%s, prec = %d, acc = %s)\n", &y, y.Text('p', 0), y.Prec(), y.Acc())
|
|
|
|
fmt.Printf("z = %.10g (%s, prec = %d, acc = %s)\n", &z, z.Text('p', 0), z.Prec(), z.Acc())
|
2015-03-20 14:35:02 -06:00
|
|
|
// Output:
|
2015-06-03 16:11:05 -06:00
|
|
|
// x = 1000 (0x.fap+10, prec = 64, acc = Exact)
|
|
|
|
// y = 2.718281828 (0x.adf85458248cd8p+2, prec = 53, acc = Exact)
|
|
|
|
// z = 1002.718282 (0x.faadf854p+10, prec = 32, acc = Below)
|
2015-03-20 14:35:02 -06:00
|
|
|
}
|
|
|
|
|
|
|
|
func Example_Shift() {
|
|
|
|
// Implementing Float "shift" by modifying the (binary) exponents directly.
|
|
|
|
for s := -5; s <= 5; s++ {
|
|
|
|
x := big.NewFloat(0.5)
|
|
|
|
x.SetMantExp(x, x.MantExp(nil)+s) // shift x by s
|
|
|
|
fmt.Println(x)
|
|
|
|
}
|
|
|
|
// Output:
|
|
|
|
// 0.015625
|
|
|
|
// 0.03125
|
|
|
|
// 0.0625
|
|
|
|
// 0.125
|
|
|
|
// 0.25
|
|
|
|
// 0.5
|
|
|
|
// 1
|
|
|
|
// 2
|
|
|
|
// 4
|
|
|
|
// 8
|
|
|
|
// 16
|
|
|
|
}
|
|
|
|
|
|
|
|
func ExampleFloat_Cmp() {
|
|
|
|
inf := math.Inf(1)
|
|
|
|
zero := 0.0
|
|
|
|
|
2015-04-02 15:47:03 -06:00
|
|
|
operands := []float64{-inf, -1.2, -zero, 0, +1.2, +inf}
|
2015-03-20 14:35:02 -06:00
|
|
|
|
2015-04-02 15:47:03 -06:00
|
|
|
fmt.Println(" x y cmp")
|
|
|
|
fmt.Println("---------------")
|
2015-03-20 14:35:02 -06:00
|
|
|
for _, x64 := range operands {
|
|
|
|
x := big.NewFloat(x64)
|
|
|
|
for _, y64 := range operands {
|
|
|
|
y := big.NewFloat(y64)
|
2015-06-03 16:11:05 -06:00
|
|
|
fmt.Printf("%4g %4g %3d\n", x, y, x.Cmp(y))
|
2015-03-20 14:35:02 -06:00
|
|
|
}
|
|
|
|
fmt.Println()
|
|
|
|
}
|
|
|
|
|
|
|
|
// Output:
|
2015-04-02 15:47:03 -06:00
|
|
|
// x y cmp
|
|
|
|
// ---------------
|
|
|
|
// -Inf -Inf 0
|
|
|
|
// -Inf -1.2 -1
|
|
|
|
// -Inf -0 -1
|
|
|
|
// -Inf 0 -1
|
|
|
|
// -Inf 1.2 -1
|
|
|
|
// -Inf +Inf -1
|
2015-03-20 14:35:02 -06:00
|
|
|
//
|
2015-04-02 15:47:03 -06:00
|
|
|
// -1.2 -Inf 1
|
|
|
|
// -1.2 -1.2 0
|
|
|
|
// -1.2 -0 -1
|
|
|
|
// -1.2 0 -1
|
|
|
|
// -1.2 1.2 -1
|
|
|
|
// -1.2 +Inf -1
|
2015-03-20 14:35:02 -06:00
|
|
|
//
|
2015-04-02 15:47:03 -06:00
|
|
|
// -0 -Inf 1
|
|
|
|
// -0 -1.2 1
|
|
|
|
// -0 -0 0
|
|
|
|
// -0 0 0
|
|
|
|
// -0 1.2 -1
|
|
|
|
// -0 +Inf -1
|
2015-03-20 14:35:02 -06:00
|
|
|
//
|
2015-04-02 15:47:03 -06:00
|
|
|
// 0 -Inf 1
|
|
|
|
// 0 -1.2 1
|
|
|
|
// 0 -0 0
|
|
|
|
// 0 0 0
|
|
|
|
// 0 1.2 -1
|
|
|
|
// 0 +Inf -1
|
2015-03-20 14:35:02 -06:00
|
|
|
//
|
2015-04-02 15:47:03 -06:00
|
|
|
// 1.2 -Inf 1
|
|
|
|
// 1.2 -1.2 1
|
|
|
|
// 1.2 -0 1
|
|
|
|
// 1.2 0 1
|
|
|
|
// 1.2 1.2 0
|
|
|
|
// 1.2 +Inf -1
|
2015-03-20 14:35:02 -06:00
|
|
|
//
|
2015-04-02 15:47:03 -06:00
|
|
|
// +Inf -Inf 1
|
|
|
|
// +Inf -1.2 1
|
|
|
|
// +Inf -0 1
|
|
|
|
// +Inf 0 1
|
|
|
|
// +Inf 1.2 1
|
|
|
|
// +Inf +Inf 0
|
2015-03-20 14:35:02 -06:00
|
|
|
}
|