File Manager

X7ROOT File Manager

Current Path: /opt/golang/1.22.0/src/math/big

opt / golang / 1.22.0 / src / math / big /

📁 ..
📄 accuracy_string.go(647 B)
📄 alias_test.go(8.81 KB)
📄 arith.go(8.28 KB)
📄 arith_386.s(4.04 KB)
📄 arith_amd64.go(278 B)
📄 arith_amd64.s(9.07 KB)
📄 arith_arm.s(4 KB)
📄 arith_arm64.s(11.85 KB)
📄 arith_decl.go(686 B)
📄 arith_decl_pure.go(1.01 KB)
📄 arith_decl_s390x.go(503 B)
📄 arith_loong64.s(749 B)
📄 arith_mips64x.s(763 B)
📄 arith_mipsx.s(759 B)
📄 arith_ppc64x.s(16.4 KB)
📄 arith_riscv64.s(750 B)
📄 arith_s390x.s(20.29 KB)
📄 arith_s390x_test.go(770 B)
📄 arith_test.go(19.88 KB)
📄 arith_wasm.s(613 B)
📄 bits_test.go(5.07 KB)
📄 calibrate_test.go(4.63 KB)
📄 decimal.go(6.63 KB)
📄 decimal_test.go(3.33 KB)
📄 doc.go(3.8 KB)
📄 example_rat_test.go(1.68 KB)
📄 example_test.go(4.05 KB)
📄 float.go(44.55 KB)
📄 float_test.go(51.94 KB)
📄 floatconv.go(8.35 KB)
📄 floatconv_test.go(24.27 KB)
📄 floatexample_test.go(3.63 KB)
📄 floatmarsh.go(3.65 KB)
📄 floatmarsh_test.go(4.51 KB)
📄 ftoa.go(13.5 KB)
📄 gcd_test.go(2.16 KB)
📄 hilbert_test.go(2.88 KB)
📄 int.go(33.21 KB)
📄 int_test.go(58.47 KB)
📄 intconv.go(6.7 KB)
📄 intconv_test.go(10.01 KB)
📄 intmarsh.go(2.19 KB)
📄 intmarsh_test.go(3.07 KB)
📄 link_test.go(1.4 KB)
📄 nat.go(31.81 KB)
📄 nat_test.go(26.22 KB)
📄 natconv.go(14.56 KB)
📄 natconv_test.go(16.85 KB)
📄 natdiv.go(34.36 KB)
📄 prime.go(10.39 KB)
📄 prime_test.go(7.1 KB)
📄 rat.go(13.48 KB)
📄 rat_test.go(18.89 KB)
📄 ratconv.go(12.35 KB)
📄 ratconv_test.go(19.31 KB)
📄 ratmarsh.go(2.23 KB)
📄 ratmarsh_test.go(3.34 KB)
📄 roundingmode_string.go(819 B)
📄 sqrt.go(2.79 KB)
📄 sqrt_test.go(4.81 KB)

Editing: hilbert_test.go

// Copyright 2009 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.

// A little test program and benchmark for rational arithmetics.
// Computes a Hilbert matrix, its inverse, multiplies them
// and verifies that the product is the identity matrix.

package big

import (
	"fmt"
	"testing"
)

type matrix struct {
	n, m int
	a    []*Rat
}

func (a *matrix) at(i, j int) *Rat {
	if !(0 <= i && i < a.n && 0 <= j && j < a.m) {
		panic("index out of range")
	}
	return a.a[i*a.m+j]
}

func (a *matrix) set(i, j int, x *Rat) {
	if !(0 <= i && i < a.n && 0 <= j && j < a.m) {
		panic("index out of range")
	}
	a.a[i*a.m+j] = x
}

func newMatrix(n, m int) *matrix {
	if !(0 <= n && 0 <= m) {
		panic("illegal matrix")
	}
	a := new(matrix)
	a.n = n
	a.m = m
	a.a = make([]*Rat, n*m)
	return a
}

func newUnit(n int) *matrix {
	a := newMatrix(n, n)
	for i := 0; i < n; i++ {
		for j := 0; j < n; j++ {
			x := NewRat(0, 1)
			if i == j {
				x.SetInt64(1)
			}
			a.set(i, j, x)
		}
	}
	return a
}

func newHilbert(n int) *matrix {
	a := newMatrix(n, n)
	for i := 0; i < n; i++ {
		for j := 0; j < n; j++ {
			a.set(i, j, NewRat(1, int64(i+j+1)))
		}
	}
	return a
}

func newInverseHilbert(n int) *matrix {
	a := newMatrix(n, n)
	for i := 0; i < n; i++ {
		for j := 0; j < n; j++ {
			x1 := new(Rat).SetInt64(int64(i + j + 1))
			x2 := new(Rat).SetInt(new(Int).Binomial(int64(n+i), int64(n-j-1)))
			x3 := new(Rat).SetInt(new(Int).Binomial(int64(n+j), int64(n-i-1)))
			x4 := new(Rat).SetInt(new(Int).Binomial(int64(i+j), int64(i)))

x1.Mul(x1, x2)
			x1.Mul(x1, x3)
			x1.Mul(x1, x4)
			x1.Mul(x1, x4)

if (i+j)&1 != 0 {
				x1.Neg(x1)
			}

a.set(i, j, x1)
		}
	}
	return a
}

func (a *matrix) mul(b *matrix) *matrix {
	if a.m != b.n {
		panic("illegal matrix multiply")
	}
	c := newMatrix(a.n, b.m)
	for i := 0; i < c.n; i++ {
		for j := 0; j < c.m; j++ {
			x := NewRat(0, 1)
			for k := 0; k < a.m; k++ {
				x.Add(x, new(Rat).Mul(a.at(i, k), b.at(k, j)))
			}
			c.set(i, j, x)
		}
	}
	return c
}

func (a *matrix) eql(b *matrix) bool {
	if a.n != b.n || a.m != b.m {
		return false
	}
	for i := 0; i < a.n; i++ {
		for j := 0; j < a.m; j++ {
			if a.at(i, j).Cmp(b.at(i, j)) != 0 {
				return false
			}
		}
	}
	return true
}

func (a *matrix) String() string {
	s := ""
	for i := 0; i < a.n; i++ {
		for j := 0; j < a.m; j++ {
			s += fmt.Sprintf("\t%s", a.at(i, j))
		}
		s += "\n"
	}
	return s
}

func doHilbert(t *testing.T, n int) {
	a := newHilbert(n)
	b := newInverseHilbert(n)
	I := newUnit(n)
	ab := a.mul(b)
	if !ab.eql(I) {
		if t == nil {
			panic("Hilbert failed")
		}
		t.Errorf("a   = %s\n", a)
		t.Errorf("b   = %s\n", b)
		t.Errorf("a*b = %s\n", ab)
		t.Errorf("I   = %s\n", I)
	}
}

func TestHilbert(t *testing.T) {
	doHilbert(t, 10)
}

func BenchmarkHilbert(b *testing.B) {
	for i := 0; i < b.N; i++ {
		doHilbert(nil, 10)
	}
}

X7ROOT File Manager

Editing: hilbert_test.go

Upload File

Create Folder