func (e *baseFieldElem) reduce() {
	if e[0] == bMask && e[1] == aMask {
		e[0], e[1] = 0, 0
	}
}
reduce sets e to zero if it is equal to p. This is the only case where e will
not naturally be reduce to canonical form.
package fourq
import (
	"fmt"
)
type baseFieldElem [2]uint64
func newBaseFieldElem() *baseFieldElem {
	return &baseFieldElem{}
}
func (e *baseFieldElem) String() string { return fmt.Sprintf("%16.16x %16.16x", e[1], e[0]) }
func (e *baseFieldElem) GoString() string {
	return fmt.Sprintf("baseFieldElem{0x%16.16x, 0x%16.16x}", e[0], e[1])
}
func (e *baseFieldElem) Bytes() [16]byte {
	return [16]byte{
		byte(e[0]), byte(e[0] >> 8), byte(e[0] >> 16), byte(e[0] >> 24),
		byte(e[0] >> 32), byte(e[0] >> 40), byte(e[0] >> 48), byte(e[0] >> 56),
		byte(e[1]), byte(e[1] >> 8), byte(e[1] >> 16), byte(e[1] >> 24),
		byte(e[1] >> 32), byte(e[1] >> 40), byte(e[1] >> 48), byte(e[1] >> 56),
	}
}
func (e *baseFieldElem) Set(a *baseFieldElem) { e[0], e[1] = a[0], a[1] }
func (e *baseFieldElem) SetZero()             { e[0], e[1] = 0, 0 }
func (e *baseFieldElem) SetOne()              { e[0], e[1] = 1, 0 }
func (e *baseFieldElem) SetBytes(in []byte) {
	e[0] = uint64(in[0]) | uint64(in[1])<<8 | uint64(in[2])<<16 |
		uint64(in[3])<<24 | uint64(in[4])<<32 | uint64(in[5])<<40 |
		uint64(in[6])<<48 | uint64(in[7])<<56
	e[1] = uint64(in[8]) | uint64(in[9])<<8 | uint64(in[10])<<16 |
		uint64(in[11])<<24 | uint64(in[12])<<32 | uint64(in[13])<<40 |
		uint64(in[14])<<48 | uint64(in[15])<<56
}
func (e *baseFieldElem) IsZero() bool { return e[0] == 0 && e[1] == 0 }
func (e *baseFieldElem) Neg(a *baseFieldElem) *baseFieldElem {
	e[0] = ^a[0]
	e[1] = (^a[1]) & aMask
	return e
}
func (e *baseFieldElem) chain1251(a *baseFieldElem) *baseFieldElem {
	t1 := newBaseFieldElem()
	t2 := newBaseFieldElem()
	t3 := newBaseFieldElem()
	t4 := newBaseFieldElem()
	t5 := newBaseFieldElem()
	bfeSquare(t2, a)
	bfeMul(t2, t2, a)
	bfeSquare(t3, t2)
	bfeSquare(t3, t3)
	bfeMul(t3, t3, t2)
	bfeSquare(t4, t3)
	bfeSquare(t4, t4)
	bfeSquare(t4, t4)
	bfeSquare(t4, t4)
	bfeMul(t4, t4, t3)
	bfeSquare(t5, t4)
	for i := 0; i < 7; i++ {
		bfeSquare(t5, t5)
	}
	bfeMul(t5, t5, t4)
	bfeSquare(t2, t5)
	for i := 0; i < 15; i++ {
		bfeSquare(t2, t2)
	}
	bfeMul(t2, t2, t5)
	bfeSquare(t1, t2)
	for i := 0; i < 31; i++ {
		bfeSquare(t1, t1)
	}
	bfeMul(t1, t1, t2)
	for i := 0; i < 32; i++ {
		bfeSquare(t1, t1)
	}
	bfeMul(t1, t2, t1)
	for i := 0; i < 16; i++ {
		bfeSquare(t1, t1)
	}
	bfeMul(t1, t1, t5)
	for i := 0; i < 8; i++ {
		bfeSquare(t1, t1)
	}
	bfeMul(t1, t1, t4)
	for i := 0; i < 4; i++ {
		bfeSquare(t1, t1)
	}
	bfeMul(t1, t1, t3)
	bfeSquare(t1, t1)
	bfeMul(e, t1, a)
	return e
}
func (e *baseFieldElem) Invert(a *baseFieldElem) *baseFieldElem {
	t := newBaseFieldElem().chain1251(a)
	bfeSquare(t, t)
	bfeSquare(t, t)
	bfeMul(e, t, a)
	return e
}