move wfa.js glue script to root,

add minify wfa.js job to build command

Makefile

@@ -2,16 +2,25 @@
 
 build: clean
 	@echo "======================== Building Binary ======================="
+	minify wfa.js > dist/wfa.js
 	GOOS=js GOARCH=wasm CGO_ENABLED=0 tinygo build -no-debug -opt=2 -target=wasm -o dist/wfa.wasm .
 
 clean:
 	@echo "======================== Cleaning Project ======================"
 	go clean
-	rm -f dist/wfa.wasm
+	rm -f dist/wfa.wasm dist/wfa.js cover.prof cpu.prof mem.prof test.test
 
 test:
 	@echo "======================== Running Tests ========================="
-	go test -v -cover -coverpkg=./pkg/ -coverprofile coverage ./test/
+	go test -v -cover -coverpkg=./pkg/ -coverprofile cover.prof -cpuprofile cpu.prof -memprofile mem.prof ./test/
 	@echo "======================= Coverage Report ========================"
-	go tool cover -func=coverage
-	@rm -f coverage
+	go tool cover -func=cover.prof
+	@rm -f cover.prof
+	@echo "==================== CPU Performance Report ===================="
+	go tool pprof -top cpu.prof
+	@rm -f cpu.prof
+	@echo "=================== Memory Performance Report =================="
+	go tool pprof -top mem.prof
+	@rm -f mem.prof
+
+	@rm -f test.test

pkg/custom_slice.go

@@ -16,16 +16,12 @@ func (a *IntegerSlice[T]) TranslateIndex(idx int) int {
 
 func (a *IntegerSlice[T]) Valid(idx int) bool {
 	actualIdx := a.TranslateIndex(idx)
-	if actualIdx < len(a.valid) { // idx is in the slice
-		return a.valid[actualIdx]
-	} else { // idx is out of the slice
-		return false
-	}
+	return 0 <= actualIdx && actualIdx < len(a.valid) && a.valid[actualIdx]
 }
 
 func (a *IntegerSlice[T]) Get(idx int) T {
 	actualIdx := a.TranslateIndex(idx)
-	if actualIdx < len(a.valid) { // idx is in the slice
+	if 0 <= actualIdx && actualIdx < len(a.valid) && a.valid[actualIdx] { // idx is in the slice
 		return a.data[actualIdx]
 	} else { // idx is out of the slice
 		return a.defaultValue
@@ -36,12 +32,12 @@ func (a *IntegerSlice[T]) Set(idx int, value T) {
 	actualIdx := a.TranslateIndex(idx)
 	if actualIdx >= len(a.valid) { // idx is outside the slice
 		// expand data array to actualIdx
-		newData := make([]T, actualIdx+1)
+		newData := make([]T, 2*actualIdx+1)
 		copy(newData, a.data)
 		a.data = newData
 
 		// expand valid array to actualIdx
-		newValid := make([]bool, actualIdx+1)
+		newValid := make([]bool, 2*actualIdx+1)
 		copy(newValid, a.valid)
 		a.valid = newValid
 	}
@@ -50,6 +46,20 @@ func (a *IntegerSlice[T]) Set(idx int, value T) {
 	a.valid[actualIdx] = true
 }
 
+func (a *IntegerSlice[T]) Preallocate(lo int, hi int) {
+	actualLo := a.TranslateIndex(lo)
+	actualHi := a.TranslateIndex(hi)
+	size := max(actualHi, actualLo)
+
+	// expand data array to actualIdx
+	newData := make([]T, size+1)
+	a.data = newData
+
+	// expand valid array to actualIdx
+	newValid := make([]bool, size+1)
+	a.valid = newValid
+}
+
 type PositiveSlice[T any] struct {
 	data         []T
 	valid        []bool
@@ -62,16 +72,12 @@ func (a *PositiveSlice[T]) TranslateIndex(idx int) int {
 
 func (a *PositiveSlice[T]) Valid(idx int) bool {
 	actualIdx := a.TranslateIndex(idx)
-	if actualIdx >= 0 && actualIdx < len(a.valid) { // idx is in the slice
-		return a.valid[actualIdx]
-	} else { // idx is out of the slice
-		return false
-	}
+	return 0 <= actualIdx && actualIdx < len(a.valid) && a.valid[actualIdx]
 }
 
 func (a *PositiveSlice[T]) Get(idx int) T {
 	actualIdx := a.TranslateIndex(idx)
-	if actualIdx >= 0 && actualIdx < len(a.valid) { // idx is in the slice
+	if 0 <= actualIdx && actualIdx < len(a.valid) && a.valid[actualIdx] { // idx is in the slice
 		return a.data[actualIdx]
 	} else { // idx is out of the slice
 		return a.defaultValue
@@ -82,12 +88,12 @@ func (a *PositiveSlice[T]) Set(idx int, value T) {
 	actualIdx := a.TranslateIndex(idx)
 	if actualIdx < 0 || actualIdx >= len(a.valid) { // idx is outside the slice
 		// expand data array to actualIdx
-		newData := make([]T, actualIdx+1)
+		newData := make([]T, 2*actualIdx+1)
 		copy(newData, a.data)
 		a.data = newData
 
 		// expand valid array to actualIdx
-		newValid := make([]bool, actualIdx+1)
+		newValid := make([]bool, 2*actualIdx+1)
 		copy(newValid, a.valid)
 		a.valid = newValid
 	}
@@ -95,3 +101,15 @@ func (a *PositiveSlice[T]) Set(idx int, value T) {
 	a.data[actualIdx] = value
 	a.valid[actualIdx] = true
 }
+
+func (a *PositiveSlice[T]) Preallocate(hi int) {
+	size := hi
+
+	// expand data array to actualIdx
+	newData := make([]T, size+1)
+	a.data = newData
+
+	// expand valid array to actualIdx
+	newValid := make([]bool, size+1)
+	a.valid = newValid
+}

pkg/types.go

@@ -5,6 +5,11 @@ import (
 	"math"
 )
 
+type Result struct {
+	Score int
+	CIGAR string
+}
+
 type Penalty struct {
 	M int
 	X int
@@ -32,14 +37,14 @@ type WavefrontComponent struct {
 	A  *PositiveSlice[*IntegerSlice[traceback]] // compact CIGAR for backtrace for each wavefront
 }
 
-func NewWavefrontComponent() WavefrontComponent {
+func NewWavefrontComponent(preallocateSize int) WavefrontComponent {
 	// new wavefront component = {
 	// lo = [0]
 	// hi = [0]
 	// W = []
 	// A = []
 	// }
-	return WavefrontComponent{
+	w := WavefrontComponent{
 		lo: &PositiveSlice[int]{
 			data:  []int{0},
 			valid: []bool{true},
@@ -48,98 +53,73 @@ func NewWavefrontComponent() WavefrontComponent {
 			data:  []int{0},
 			valid: []bool{true},
 		},
-		W: &PositiveSlice[*IntegerSlice[int]]{},
-		A: &PositiveSlice[*IntegerSlice[traceback]]{},
+		W: &PositiveSlice[*IntegerSlice[int]]{
+			defaultValue: &IntegerSlice[int]{
+				data:  []int{},
+				valid: []bool{},
+			},
+		},
+		A: &PositiveSlice[*IntegerSlice[traceback]]{
+			defaultValue: &IntegerSlice[traceback]{
+				data:  []traceback{},
+				valid: []bool{},
+			},
+		},
 	}
+
+	w.lo.Preallocate(preallocateSize)
+	w.hi.Preallocate(preallocateSize)
+	w.W.Preallocate(preallocateSize)
+	w.A.Preallocate(preallocateSize)
+
+	return w
 }
 
 // get value for wavefront=score, diag=k => returns ok, value
 func (w *WavefrontComponent) GetVal(score int, k int) (bool, int) {
-	// if W[score][k] is valid
-	if w.W.Valid(score) && w.W.Get(score).Valid(k) {
-		// return W[score][k]
-		return true, w.W.Get(score).Get(k)
-	} else {
-		return false, 0
-	}
+	return w.W.Valid(score) && w.W.Get(score).Valid(k), w.W.Get(score).Get(k)
 }
 
 // set value for wavefront=score, diag=k
 func (w *WavefrontComponent) SetVal(score int, k int, val int) {
-	// if W[score] is valid
-	if w.W.Valid(score) {
-		// W[score][k] = val
-		w.W.Get(score).Set(k, val)
-	} else {
-		// W[score] = []
-		w.W.Set(score, &IntegerSlice[int]{})
-		// W[score][k] = val
-		w.W.Get(score).Set(k, val)
-	}
+	w.W.Get(score).Set(k, val)
 }
 
 // get alignment traceback for wavefront=score, diag=k => returns ok, value
 func (w *WavefrontComponent) GetTraceback(score int, k int) (bool, traceback) {
-	// if W[score][k] is valid
-	if w.A.Valid(score) && w.A.Get(score).Valid(k) {
-		// return W[score][k]
-		return true, w.A.Get(score).Get(k)
-	} else {
-		return false, 0
-	}
+	return w.A.Valid(score) && w.A.Get(score).Valid(k), w.A.Get(score).Get(k)
 }
 
 // set alignment traceback for wavefront=score, diag=k
 func (w *WavefrontComponent) SetTraceback(score int, k int, val traceback) {
-	// if A[score] is valid
-	if w.A.Valid(score) {
-		// A[score][k] = val
-		w.A.Get(score).Set(k, val)
-	} else {
-		// W[score] = []
-		w.A.Set(score, &IntegerSlice[traceback]{})
-		// W[score][k] = val
-		w.A.Get(score).Set(k, val)
-	}
+	w.A.Get(score).Set(k, val)
 }
 
 // get hi for wavefront=score
-func (w *WavefrontComponent) GetHi(score int) (bool, int) {
-	// if hi[score] is valid
-	if w.hi.Valid(score) {
-		// return hi[score]
-		return true, w.hi.Get(score)
+func (w *WavefrontComponent) GetLoHi(score int) (bool, int, int) {
+	// if lo[score] and hi[score] are valid
+	if w.lo.Valid(score) && w.hi.Valid(score) {
+		// return lo[score] hi[score]
+		return true, w.lo.Get(score), w.hi.Get(score)
 	} else {
-		return false, 0
+		return false, 0, 0
 	}
 }
 
 // set hi for wavefront=score
-func (w *WavefrontComponent) SetHi(score int, hi int) {
-	// hi[score] = hi
-	w.hi.Set(score, hi)
-}
-
-// get lo for wavefront=score
-func (w *WavefrontComponent) GetLo(score int) (bool, int) {
-	// if lo[score] is valid
-	if w.lo.Valid(score) {
-		// return lo[score]
-		return true, w.lo.Get(score)
-	} else {
-		return false, 0
-	}
-}
-
-// set hi for wavefront=score
-func (w *WavefrontComponent) SetLo(score int, lo int) {
+func (w *WavefrontComponent) SetLoHi(score int, lo int, hi int) {
 	// lo[score] = lo
 	w.lo.Set(score, lo)
-}
+	// hi[score] = hi
+	w.hi.Set(score, hi)
 
-type Result struct {
-	Score int
-	CIGAR string
+	// preemptively setup w.A
+	w.A.Set(score, &IntegerSlice[traceback]{})
+	w.A.Get(score).Preallocate(lo, hi)
+
+	// preemptively setup w.W
+	w.W.Set(score, &IntegerSlice[int]{})
+	w.W.Get(score).Preallocate(lo, hi)
 }
 
 func (w *WavefrontComponent) String(score int) string {

pkg/utils.go

@@ -5,14 +5,12 @@ import (
 	"unicode/utf8"
 )
 
-func SafeMin(valids []bool, values []int) (bool, int) {
-	ok, idx := SafeArgMin(valids, values)
-	return ok, values[idx]
+func SafeMin(values []int, idx int) int {
+	return values[idx]
 }
 
-func SafeMax(valids []bool, values []int) (bool, int) {
-	ok, idx := SafeArgMax(valids, values)
-	return ok, values[idx]
+func SafeMax(values []int, idx int) int {
+	return values[idx]
 }
 
 func SafeArgMax(valids []bool, values []int) (bool, int) {
@@ -66,50 +64,34 @@ func Splice(s string, c rune, idx int) string {
 	return s[:idx] + string(c) + s[idx:]
 }
 
-func NextLo(M WavefrontComponent, I WavefrontComponent, D WavefrontComponent, score int, penalties Penalty) int {
+func NextLoHi(M WavefrontComponent, I WavefrontComponent, D WavefrontComponent, score int, penalties Penalty) (int, int) {
 	x := penalties.X
 	o := penalties.O
 	e := penalties.E
 
-	a_ok, a := M.GetLo(score - x)
-	b_ok, b := M.GetLo(score - o - e)
-	c_ok, c := I.GetLo(score - e)
-	d_ok, d := D.GetLo(score - e)
+	a_ok, a_lo, a_hi := M.GetLoHi(score - x)
+	b_ok, b_lo, b_hi := M.GetLoHi(score - o - e)
+	c_ok, c_lo, c_hi := I.GetLoHi(score - e)
+	d_ok, d_lo, d_hi := D.GetLoHi(score - e)
 
-	ok, lo := SafeMin(
+	ok_lo, idx := SafeArgMin(
 		[]bool{a_ok, b_ok, c_ok, d_ok},
-		[]int{a, b, c, d},
+		[]int{a_lo, b_lo, c_lo, d_lo},
 	)
-	lo--
-	if ok {
-		M.SetLo(score, lo)
-		I.SetLo(score, lo)
-		D.SetLo(score, lo)
-	}
-	return lo
-}
+	lo := SafeMin([]int{a_lo, b_lo, c_lo, d_lo}, idx) - 1
 
-func NextHi(M WavefrontComponent, I WavefrontComponent, D WavefrontComponent, score int, penalties Penalty) int {
-	x := penalties.X
-	o := penalties.O
-	e := penalties.E
-
-	a_ok, a := M.GetHi(score - x)
-	b_ok, b := M.GetHi(score - o - e)
-	c_ok, c := I.GetHi(score - e)
-	d_ok, d := D.GetHi(score - e)
-
-	ok, hi := SafeMax(
+	ok_hi, idx := SafeArgMax(
 		[]bool{a_ok, b_ok, c_ok, d_ok},
-		[]int{a, b, c, d},
+		[]int{a_hi, b_hi, c_hi, d_hi},
 	)
-	hi++
-	if ok {
-		M.SetHi(score, hi)
-		I.SetHi(score, hi)
-		D.SetHi(score, hi)
+	hi := SafeMax([]int{a_hi, b_hi, c_hi, d_hi}, idx) + 1
+
+	if ok_lo && ok_hi {
+		M.SetLoHi(score, lo, hi)
+		I.SetLoHi(score, lo, hi)
+		D.SetLoHi(score, lo, hi)
 	}
-	return hi
+	return lo, hi
 }
 
 func NextI(M WavefrontComponent, I WavefrontComponent, score int, k int, penalties Penalty) {
@@ -119,13 +101,10 @@ func NextI(M WavefrontComponent, I WavefrontComponent, score int, k int, penalti
 	a_ok, a := M.GetVal(score-o-e, k-1)
 	b_ok, b := I.GetVal(score-e, k-1)
 
-	ok, nextIVal := SafeMax([]bool{a_ok, b_ok}, []int{a, b})
-	if ok {
-		I.SetVal(score, k, nextIVal+1) // important that the +1 is here
-	}
-
 	ok, nextITraceback := SafeArgMax([]bool{a_ok, b_ok}, []int{a, b})
+	nextIVal := SafeMax([]int{a, b}, nextITraceback) + 1 // important that the +1 is here
 	if ok {
+		I.SetVal(score, k, nextIVal)
 		I.SetTraceback(score, k, []traceback{OpenIns, ExtdIns}[nextITraceback])
 	}
 }
@@ -137,13 +116,13 @@ func NextD(M WavefrontComponent, D WavefrontComponent, score int, k int, penalti
 	a_ok, a := M.GetVal(score-o-e, k+1)
 	b_ok, b := D.GetVal(score-e, k+1)
 
-	ok, nextDVal := SafeMax([]bool{a_ok, b_ok}, []int{a, b})
-	if ok {
-		D.SetVal(score, k, nextDVal) // nothing special
-	}
-
-	ok, nextDTraceback := SafeArgMax([]bool{a_ok, b_ok}, []int{a, b})
+	ok, nextDTraceback := SafeArgMax(
+		[]bool{a_ok, b_ok},
+		[]int{a, b},
+	)
+	nextDVal := SafeMax([]int{a, b}, nextDTraceback) // nothing special
 	if ok {
+		D.SetVal(score, k, nextDVal)
 		D.SetTraceback(score, k, []traceback{OpenDel, ExtdDel}[nextDTraceback])
 	}
 }
@@ -156,13 +135,11 @@ func NextM(M WavefrontComponent, I WavefrontComponent, D WavefrontComponent, sco
 	b_ok, b := I.GetVal(score, k)
 	c_ok, c := D.GetVal(score, k)
 
-	ok, nextMVal := SafeMax([]bool{a_ok, b_ok, c_ok}, []int{a, b, c})
+	ok, nextMTraceback := SafeArgMax([]bool{a_ok, b_ok, c_ok}, []int{a, b, c})
+	nextMVal := SafeMax([]int{a, b, c}, nextMTraceback)
+
 	if ok {
 		M.SetVal(score, k, nextMVal)
-	}
-
-	ok, nextMTraceback := SafeArgMax([]bool{a_ok, b_ok, c_ok}, []int{a, b, c})
-	if ok {
 		M.SetTraceback(score, k, []traceback{Sub, Ins, Del}[nextMTraceback])
 	}
 }

pkg/wfa.go

@@ -6,13 +6,13 @@ func WFAlign(s1 string, s2 string, penalties Penalty, doCIGAR bool) Result {
 	A_k := m - n
 	A_offset := m
 	score := 0
-	M := NewWavefrontComponent()
+	estimatedScore := (max(n, m) * max(penalties.M, penalties.X, penalties.O, penalties.E)) / 4
+	M := NewWavefrontComponent(estimatedScore)
+	M.SetLoHi(0, 0, 0)
 	M.SetVal(0, 0, 0)
-	M.SetHi(0, 0)
-	M.SetLo(0, 0)
 	M.SetTraceback(0, 0, End)
-	I := NewWavefrontComponent()
-	D := NewWavefrontComponent()
+	I := NewWavefrontComponent(estimatedScore)
+	D := NewWavefrontComponent(estimatedScore)
 
 	for {
 		WFExtend(M, s1, n, s2, m, score)
@@ -36,8 +36,7 @@ func WFAlign(s1 string, s2 string, penalties Penalty, doCIGAR bool) Result {
 }
 
 func WFExtend(M WavefrontComponent, s1 string, n int, s2 string, m int, score int) {
-	_, lo := M.GetLo(score)
-	_, hi := M.GetHi(score)
+	_, lo, hi := M.GetLoHi(score)
 	for k := lo; k <= hi; k++ {
 		// v = M[score][k] - k
 		// h = M[score][k]
@@ -58,11 +57,8 @@ func WFExtend(M WavefrontComponent, s1 string, n int, s2 string, m int, score in
 }
 
 func WFNext(M WavefrontComponent, I WavefrontComponent, D WavefrontComponent, score int, penalties Penalty) {
-	// get this score's lo
-	lo := NextLo(M, I, D, score, penalties)
-
-	// get this score's hi
-	hi := NextHi(M, I, D, score, penalties)
+	// get this score's lo, hi
+	lo, hi := NextLoHi(M, I, D, score, penalties)
 
 	for k := lo; k <= hi; k++ {
 		NextI(M, I, score, k, penalties)

wfa.js

@@ -0,0 +1,556 @@
+// wasm_exec.js from tinygo
+
+// Copyright 2018 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.
+//
+// This file has been modified for use by the TinyGo compiler.
+
+(() => {
+	// Map multiple JavaScript environments to a single common API,
+	// preferring web standards over Node.js API.
+	//
+	// Environments considered:
+	// - Browsers
+	// - Node.js
+	// - Electron
+	// - Parcel
+
+	if (typeof global !== "undefined") {
+		// global already exists
+	} else if (typeof window !== "undefined") {
+		window.global = window;
+	} else if (typeof self !== "undefined") {
+		self.global = self;
+	} else {
+		throw new Error("cannot export Go (neither global, window nor self is defined)");
+	}
+
+	if (!global.require && typeof require !== "undefined") {
+		global.require = require;
+	}
+
+	if (!global.fs && global.require) {
+		global.fs = require("node:fs");
+	}
+
+	const enosys = () => {
+		const err = new Error("not implemented");
+		err.code = "ENOSYS";
+		return err;
+	};
+
+	if (!global.fs) {
+		let outputBuf = "";
+		global.fs = {
+			constants: { O_WRONLY: -1, O_RDWR: -1, O_CREAT: -1, O_TRUNC: -1, O_APPEND: -1, O_EXCL: -1 }, // unused
+			writeSync(fd, buf) {
+				outputBuf += decoder.decode(buf);
+				const nl = outputBuf.lastIndexOf("\n");
+				if (nl != -1) {
+					console.log(outputBuf.substr(0, nl));
+					outputBuf = outputBuf.substr(nl + 1);
+				}
+				return buf.length;
+			},
+			write(fd, buf, offset, length, position, callback) {
+				if (offset !== 0 || length !== buf.length || position !== null) {
+					callback(enosys());
+					return;
+				}
+				const n = this.writeSync(fd, buf);
+				callback(null, n);
+			},
+			chmod(path, mode, callback) { callback(enosys()); },
+			chown(path, uid, gid, callback) { callback(enosys()); },
+			close(fd, callback) { callback(enosys()); },
+			fchmod(fd, mode, callback) { callback(enosys()); },
+			fchown(fd, uid, gid, callback) { callback(enosys()); },
+			fstat(fd, callback) { callback(enosys()); },
+			fsync(fd, callback) { callback(null); },
+			ftruncate(fd, length, callback) { callback(enosys()); },
+			lchown(path, uid, gid, callback) { callback(enosys()); },
+			link(path, link, callback) { callback(enosys()); },
+			lstat(path, callback) { callback(enosys()); },
+			mkdir(path, perm, callback) { callback(enosys()); },
+			open(path, flags, mode, callback) { callback(enosys()); },
+			read(fd, buffer, offset, length, position, callback) { callback(enosys()); },
+			readdir(path, callback) { callback(enosys()); },
+			readlink(path, callback) { callback(enosys()); },
+			rename(from, to, callback) { callback(enosys()); },
+			rmdir(path, callback) { callback(enosys()); },
+			stat(path, callback) { callback(enosys()); },
+			symlink(path, link, callback) { callback(enosys()); },
+			truncate(path, length, callback) { callback(enosys()); },
+			unlink(path, callback) { callback(enosys()); },
+			utimes(path, atime, mtime, callback) { callback(enosys()); },
+		};
+	}
+
+	if (!global.process) {
+		global.process = {
+			getuid() { return -1; },
+			getgid() { return -1; },
+			geteuid() { return -1; },
+			getegid() { return -1; },
+			getgroups() { throw enosys(); },
+			pid: -1,
+			ppid: -1,
+			umask() { throw enosys(); },
+			cwd() { throw enosys(); },
+			chdir() { throw enosys(); },
+		}
+	}
+
+	if (!global.crypto) {
+		const nodeCrypto = require("node:crypto");
+		global.crypto = {
+			getRandomValues(b) {
+				nodeCrypto.randomFillSync(b);
+			},
+		};
+	}
+
+	if (!global.performance) {
+		global.performance = {
+			now() {
+				const [sec, nsec] = process.hrtime();
+				return sec * 1000 + nsec / 1000000;
+			},
+		};
+	}
+
+	if (!global.TextEncoder) {
+		global.TextEncoder = require("node:util").TextEncoder;
+	}
+
+	if (!global.TextDecoder) {
+		global.TextDecoder = require("node:util").TextDecoder;
+	}
+
+	// End of polyfills for common API.
+
+	const encoder = new TextEncoder("utf-8");
+	const decoder = new TextDecoder("utf-8");
+	let reinterpretBuf = new DataView(new ArrayBuffer(8));
+	var logLine = [];
+
+	global.Go = class {
+		constructor() {
+			this._callbackTimeouts = new Map();
+			this._nextCallbackTimeoutID = 1;
+
+			const mem = () => {
+				// The buffer may change when requesting more memory.
+				return new DataView(this._inst.exports.memory.buffer);
+			}
+
+			const unboxValue = (v_ref) => {
+				reinterpretBuf.setBigInt64(0, v_ref, true);
+				const f = reinterpretBuf.getFloat64(0, true);
+				if (f === 0) {
+					return undefined;
+				}
+				if (!isNaN(f)) {
+					return f;
+				}
+
+				const id = v_ref & 0xffffffffn;
+				return this._values[id];
+			}
+
+
+			const loadValue = (addr) => {
+				let v_ref = mem().getBigUint64(addr, true);
+				return unboxValue(v_ref);
+			}
+
+			const boxValue = (v) => {
+				const nanHead = 0x7FF80000n;
+
+				if (typeof v === "number") {
+					if (isNaN(v)) {
+						return nanHead << 32n;
+					}
+					if (v === 0) {
+						return (nanHead << 32n) | 1n;
+					}
+					reinterpretBuf.setFloat64(0, v, true);
+					return reinterpretBuf.getBigInt64(0, true);
+				}
+
+				switch (v) {
+					case undefined:
+						return 0n;
+					case null:
+						return (nanHead << 32n) | 2n;
+					case true:
+						return (nanHead << 32n) | 3n;
+					case false:
+						return (nanHead << 32n) | 4n;
+				}
+
+				let id = this._ids.get(v);
+				if (id === undefined) {
+					id = this._idPool.pop();
+					if (id === undefined) {
+						id = BigInt(this._values.length);
+					}
+					this._values[id] = v;
+					this._goRefCounts[id] = 0;
+					this._ids.set(v, id);
+				}
+				this._goRefCounts[id]++;
+				let typeFlag = 1n;
+				switch (typeof v) {
+					case "string":
+						typeFlag = 2n;
+						break;
+					case "symbol":
+						typeFlag = 3n;
+						break;
+					case "function":
+						typeFlag = 4n;
+						break;
+				}
+				return id | ((nanHead | typeFlag) << 32n);
+			}
+
+			const storeValue = (addr, v) => {
+				let v_ref = boxValue(v);
+				mem().setBigUint64(addr, v_ref, true);
+			}
+
+			const loadSlice = (array, len, cap) => {
+				return new Uint8Array(this._inst.exports.memory.buffer, array, len);
+			}
+
+			const loadSliceOfValues = (array, len, cap) => {
+				const a = new Array(len);
+				for (let i = 0; i < len; i++) {
+					a[i] = loadValue(array + i * 8);
+				}
+				return a;
+			}
+
+			const loadString = (ptr, len) => {
+				return decoder.decode(new DataView(this._inst.exports.memory.buffer, ptr, len));
+			}
+
+			const timeOrigin = Date.now() - performance.now();
+			this.importObject = {
+				wasi_snapshot_preview1: {
+					// https://github.com/WebAssembly/WASI/blob/main/phases/snapshot/docs.md#fd_write
+					fd_write: function(fd, iovs_ptr, iovs_len, nwritten_ptr) {
+						let nwritten = 0;
+						if (fd == 1) {
+							for (let iovs_i=0; iovs_i<iovs_len;iovs_i++) {
+								let iov_ptr = iovs_ptr+iovs_i*8; // assuming wasm32
+								let ptr = mem().getUint32(iov_ptr + 0, true);
+								let len = mem().getUint32(iov_ptr + 4, true);
+								nwritten += len;
+								for (let i=0; i<len; i++) {
+									let c = mem().getUint8(ptr+i);
+									if (c == 13) { // CR
+										// ignore
+									} else if (c == 10) { // LF
+										// write line
+										let line = decoder.decode(new Uint8Array(logLine));
+										logLine = [];
+										console.log(line);
+									} else {
+										logLine.push(c);
+									}
+								}
+							}
+						} else {
+							console.error('invalid file descriptor:', fd);
+						}
+						mem().setUint32(nwritten_ptr, nwritten, true);
+						return 0;
+					},
+					fd_close: () => 0,      // dummy
+					fd_fdstat_get: () => 0, // dummy
+					fd_seek: () => 0,       // dummy
+					"proc_exit": (code) => {
+						if (global.process) {
+							// Node.js
+							process.exit(code);
+						} else {
+							// Can't exit in a browser.
+							throw 'trying to exit with code ' + code;
+						}
+					},
+					random_get: (bufPtr, bufLen) => {
+						crypto.getRandomValues(loadSlice(bufPtr, bufLen));
+						return 0;
+					},
+				},
+				gojs: {
+					// func ticks() float64
+					"runtime.ticks": () => {
+						return timeOrigin + performance.now();
+					},
+
+					// func sleepTicks(timeout float64)
+					"runtime.sleepTicks": (timeout) => {
+						// Do not sleep, only reactivate scheduler after the given timeout.
+						setTimeout(this._inst.exports.go_scheduler, timeout);
+					},
+
+					// func finalizeRef(v ref)
+					"syscall/js.finalizeRef": (v_ref) => {
+						// Note: TinyGo does not support finalizers so this should never be
+						// called.
+						//console.error('syscall/js.finalizeRef not implemented');
+						// for whatever reason this is called by wfajs but doesnt impact the results at all??
+					},
+
+					// func stringVal(value string) ref
+					"syscall/js.stringVal": (value_ptr, value_len) => {
+						const s = loadString(value_ptr, value_len);
+						return boxValue(s);
+					},
+
+					// func valueGet(v ref, p string) ref
+					"syscall/js.valueGet": (v_ref, p_ptr, p_len) => {
+						let prop = loadString(p_ptr, p_len);
+						let v = unboxValue(v_ref);
+						let result = Reflect.get(v, prop);
+						return boxValue(result);
+					},
+
+					// func valueSet(v ref, p string, x ref)
+					"syscall/js.valueSet": (v_ref, p_ptr, p_len, x_ref) => {
+						const v = unboxValue(v_ref);
+						const p = loadString(p_ptr, p_len);
+						const x = unboxValue(x_ref);
+						Reflect.set(v, p, x);
+					},
+
+					// func valueDelete(v ref, p string)
+					"syscall/js.valueDelete": (v_ref, p_ptr, p_len) => {
+						const v = unboxValue(v_ref);
+						const p = loadString(p_ptr, p_len);
+						Reflect.deleteProperty(v, p);
+					},
+
+					// func valueIndex(v ref, i int) ref
+					"syscall/js.valueIndex": (v_ref, i) => {
+						return boxValue(Reflect.get(unboxValue(v_ref), i));
+					},
+
+					// valueSetIndex(v ref, i int, x ref)
+					"syscall/js.valueSetIndex": (v_ref, i, x_ref) => {
+						Reflect.set(unboxValue(v_ref), i, unboxValue(x_ref));
+					},
+
+					// func valueCall(v ref, m string, args []ref) (ref, bool)
+					"syscall/js.valueCall": (ret_addr, v_ref, m_ptr, m_len, args_ptr, args_len, args_cap) => {
+						const v = unboxValue(v_ref);
+						const name = loadString(m_ptr, m_len);
+						const args = loadSliceOfValues(args_ptr, args_len, args_cap);
+						try {
+							const m = Reflect.get(v, name);
+							storeValue(ret_addr, Reflect.apply(m, v, args));
+							mem().setUint8(ret_addr + 8, 1);
+						} catch (err) {
+							storeValue(ret_addr, err);
+							mem().setUint8(ret_addr + 8, 0);
+						}
+					},
+
+					// func valueInvoke(v ref, args []ref) (ref, bool)
+					"syscall/js.valueInvoke": (ret_addr, v_ref, args_ptr, args_len, args_cap) => {
+						try {
+							const v = unboxValue(v_ref);
+							const args = loadSliceOfValues(args_ptr, args_len, args_cap);
+							storeValue(ret_addr, Reflect.apply(v, undefined, args));
+							mem().setUint8(ret_addr + 8, 1);
+						} catch (err) {
+							storeValue(ret_addr, err);
+							mem().setUint8(ret_addr + 8, 0);
+						}
+					},
+
+					// func valueNew(v ref, args []ref) (ref, bool)
+					"syscall/js.valueNew": (ret_addr, v_ref, args_ptr, args_len, args_cap) => {
+						const v = unboxValue(v_ref);
+						const args = loadSliceOfValues(args_ptr, args_len, args_cap);
+						try {
+							storeValue(ret_addr, Reflect.construct(v, args));
+							mem().setUint8(ret_addr + 8, 1);
+						} catch (err) {
+							storeValue(ret_addr, err);
+							mem().setUint8(ret_addr+ 8, 0);
+						}
+					},
+
+					// func valueLength(v ref) int
+					"syscall/js.valueLength": (v_ref) => {
+						return unboxValue(v_ref).length;
+					},
+
+					// valuePrepareString(v ref) (ref, int)
+					"syscall/js.valuePrepareString": (ret_addr, v_ref) => {
+						const s = String(unboxValue(v_ref));
+						const str = encoder.encode(s);
+						storeValue(ret_addr, str);
+						mem().setInt32(ret_addr + 8, str.length, true);
+					},
+
+					// valueLoadString(v ref, b []byte)
+					"syscall/js.valueLoadString": (v_ref, slice_ptr, slice_len, slice_cap) => {
+						const str = unboxValue(v_ref);
+						loadSlice(slice_ptr, slice_len, slice_cap).set(str);
+					},
+
+					// func valueInstanceOf(v ref, t ref) bool
+					"syscall/js.valueInstanceOf": (v_ref, t_ref) => {
+ 						return unboxValue(v_ref) instanceof unboxValue(t_ref);
+					},
+
+					// func copyBytesToGo(dst []byte, src ref) (int, bool)
+					"syscall/js.copyBytesToGo": (ret_addr, dest_addr, dest_len, dest_cap, src_ref) => {
+						let num_bytes_copied_addr = ret_addr;
+						let returned_status_addr = ret_addr + 4; // Address of returned boolean status variable
+
+						const dst = loadSlice(dest_addr, dest_len);
+						const src = unboxValue(src_ref);
+						if (!(src instanceof Uint8Array || src instanceof Uint8ClampedArray)) {
+							mem().setUint8(returned_status_addr, 0); // Return "not ok" status
+							return;
+						}
+						const toCopy = src.subarray(0, dst.length);
+						dst.set(toCopy);
+						mem().setUint32(num_bytes_copied_addr, toCopy.length, true);
+						mem().setUint8(returned_status_addr, 1); // Return "ok" status
+					},
+
+					// copyBytesToJS(dst ref, src []byte) (int, bool)
+					// Originally copied from upstream Go project, then modified:
+					//   https://github.com/golang/go/blob/3f995c3f3b43033013013e6c7ccc93a9b1411ca9/misc/wasm/wasm_exec.js#L404-L416
+					"syscall/js.copyBytesToJS": (ret_addr, dst_ref, src_addr, src_len, src_cap) => {
+						let num_bytes_copied_addr = ret_addr;
+						let returned_status_addr = ret_addr + 4; // Address of returned boolean status variable
+
+						const dst = unboxValue(dst_ref);
+						const src = loadSlice(src_addr, src_len);
+						if (!(dst instanceof Uint8Array || dst instanceof Uint8ClampedArray)) {
+							mem().setUint8(returned_status_addr, 0); // Return "not ok" status
+							return;
+						}
+						const toCopy = src.subarray(0, dst.length);
+						dst.set(toCopy);
+						mem().setUint32(num_bytes_copied_addr, toCopy.length, true);
+						mem().setUint8(returned_status_addr, 1); // Return "ok" status
+					},
+				}
+			};
+
+			// Go 1.20 uses 'env'. Go 1.21 uses 'gojs'.
+			// For compatibility, we use both as long as Go 1.20 is supported.
+			this.importObject.env = this.importObject.gojs;
+		}
+
+		async run(instance) {
+			this._inst = instance;
+			this._values = [ // JS values that Go currently has references to, indexed by reference id
+				NaN,
+				0,
+				null,
+				true,
+				false,
+				global,
+				this,
+			];
+			this._goRefCounts = []; // number of references that Go has to a JS value, indexed by reference id
+			this._ids = new Map();  // mapping from JS values to reference ids
+			this._idPool = [];      // unused ids that have been garbage collected
+			this.exited = false;    // whether the Go program has exited
+
+			while (true) {
+				const callbackPromise = new Promise((resolve) => {
+					this._resolveCallbackPromise = () => {
+						if (this.exited) {
+							throw new Error("bad callback: Go program has already exited");
+						}
+						setTimeout(resolve, 0); // make sure it is asynchronous
+					};
+				});
+				this._inst.exports._start();
+				if (this.exited) {
+					break;
+				}
+				await callbackPromise;
+			}
+		}
+
+		_resume() {
+			if (this.exited) {
+				throw new Error("Go program has already exited");
+			}
+			this._inst.exports.resume();
+			if (this.exited) {
+				this._resolveExitPromise();
+			}
+		}
+
+		_makeFuncWrapper(id) {
+			const go = this;
+			return function () {
+				const event = { id: id, this: this, args: arguments };
+				go._pendingEvent = event;
+				go._resume();
+				return event.result;
+			};
+		}
+	}
+
+	if (
+		global.require &&
+		global.require.main === module &&
+		global.process &&
+		global.process.versions &&
+		!global.process.versions.electron
+	) {
+		if (process.argv.length != 3) {
+			console.error("usage: go_js_wasm_exec [wasm binary] [arguments]");
+			process.exit(1);
+		}
+
+		const go = new Go();
+		WebAssembly.instantiate(fs.readFileSync(process.argv[2]), go.importObject).then((result) => {
+			return go.run(result.instance);
+		}).catch((err) => {
+			console.error(err);
+			process.exit(1);
+		});
+	}
+})();
+
+// wasm setup code
+export default function init (path) {
+	return new Promise ((res) => {
+		const go = new Go();
+		var wasm;
+		
+		if ('instantiateStreaming' in WebAssembly) {
+			WebAssembly.instantiateStreaming(fetch(path), go.importObject).then(function (obj) {
+				wasm = obj.instance;
+				go.run(wasm);
+				res()
+			})
+		} else {
+			fetch(path).then(resp =>
+				resp.arrayBuffer()
+			).then(bytes =>
+				WebAssembly.instantiate(bytes, go.importObject).then(function (obj) {
+					wasm = obj.instance;
+					go.run(wasm);
+					res()
+				})
+			)
+		}
+	}) 
+}