script/print-mpb-powers-of-10.go

// Copyright 2020 The Wuffs Authors.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//    https://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

// +build ignore

package main

// print-mpb-powers-of-10.go prints the
// wuffs_base__private_implementation__medium_prec_bin__powers_of_10 tables.
//
// Usage: go run print-mpb-powers-of-10.go -comments

import (
	"flag"
	"fmt"
	"math/big"
	"os"
)

var (
	comments = flag.Bool("comments", false, "whether to print comments")
)

func main() {
	if err := main1(); err != nil {
		os.Stderr.WriteString(err.Error() + "\n")
		os.Exit(1)
	}
}

func main1() error {
	flag.Parse()

	const bigCount = 1 + ((+340 - -348) / 8)
	fmt.Printf("static const uint32_t "+
		"wuffs_base__private_implementation__big_powers_of_10[%d] = {\n", 3*bigCount)
	for e := -348; e <= +340; e += 8 {
		if err := do(e); err != nil {
			return err
		}
	}
	fmt.Printf("};\n\n")

	fmt.Printf("static const uint32_t " +
		"wuffs_base__private_implementation__small_powers_of_10[24] = {\n")
	for e := 0; e <= 7; e += 1 {
		if err := do(e); err != nil {
			return err
		}
	}
	fmt.Printf("};\n")

	return nil
}

var (
	one   = big.NewInt(1)
	ten   = big.NewInt(10)
	two64 = big.NewInt(0).Lsh(one, 64)
)

// N is large enough so that (1<<N) is bigger than 1e348.
const N = 2048

func do(e int) error {
	z := big.NewInt(0).Lsh(one, N)
	if e >= 0 {
		exp := big.NewInt(0).Exp(ten, big.NewInt(int64(+e)), nil)
		z.Mul(z, exp)
	} else {
		exp := big.NewInt(0).Exp(ten, big.NewInt(int64(-e)), nil)
		z.Div(z, exp)
	}

	roundUp := false
	n := int32(-N)
	for z.Cmp(two64) >= 0 {
		roundUp = z.Bit(0) > 0
		z.Rsh(z, 1)
		n++
	}
	if roundUp {
		z.Add(z, one)
	}
	hex := fmt.Sprintf("%X", z)
	if len(hex) != 16 {
		return fmt.Errorf("invalid hexadecimal representation %q", hex)
	}

	fmt.Printf("    0x%s, 0x%s, 0x%08X,", hex[8:], hex[:8], uint32(n))
	if *comments {
		fmt.Printf("  // 1e%-04d ≈ (0x%s ", e, hex)
		if n >= 0 {
			fmt.Printf("<< %4d)", +n)
		} else {
			fmt.Printf(">> %4d)", -n)
		}
	}

	fmt.Println()
	return nil
}