diff --git a/gc_latency/latency.go b/gc_latency/latency.go
new file mode 100644
index 0000000..7eb5391
--- /dev/null
+++ b/gc_latency/latency.go
@@ -0,0 +1,188 @@
+package main
+
+import (
+	"flag"
+	"fmt"
+	"math"
+	"os"
+	"runtime/trace"
+	"sort"
+	"testing"
+	"time"
+	"unsafe"
+)
+
+const (
+	bufferLen = 200000
+	msgCount0 = 1000000 // warmup
+	msgCount1 = 5000000 // run (5 million)
+)
+
+type kbyte []byte
+type circularBuffer [bufferLen]kbyte
+
+var worst time.Duration
+
+// For making sense of the bad outcome.
+var total time.Duration
+var worstIndex int
+var allStart time.Time
+var worstElapsed time.Duration
+
+var fluff kbyte
+var doFluff bool
+
+// newSlice allocates a 1k slice of bytes and initializes them all to byte(n)
+func newSlice(n int) kbyte {
+	m := make(kbyte, 1024)
+	if doFluff && n&63 == 0 {
+		fluff = make(kbyte, 1024)
+	}
+	for i := range m {
+		m[i] = byte(n)
+	}
+	return m
+}
+
+var delays []time.Duration
+
+// storeSlice stores a newly allocated 1k slice of bytes at circularBuffer[count%bufferLen]
+// (bufferLen is the length of the array circularBuffer)
+// It also checks the time needed to do this and records the worst case.
+func storeSlice(c *circularBuffer, highID int) {
+	start := time.Now()
+	c[highID%bufferLen] = newSlice(highID)
+	elapsed := time.Since(start)
+
+	candElapsed := time.Since(allStart) // Record location of worst in trace
+	if elapsed > worst {
+		worst = elapsed
+		worstIndex = highID
+		worstElapsed = candElapsed
+	}
+	total = time.Duration(total.Nanoseconds() + elapsed.Nanoseconds())
+	delays = append(delays, elapsed)
+}
+
+//go:noinline
+func work(c *circularBuffer, count int) {
+	for i := 0; i < count; i++ {
+		storeSlice(c, i)
+	}
+}
+
+//go:noinline
+func stack() {
+	var c circularBuffer
+	work(&c, msgCount1)
+}
+
+var global_c circularBuffer
+
+//go:noinline
+func global() {
+	work(&global_c, msgCount1)
+}
+
+var sink *circularBuffer
+
+//go:noinline
+func heap() {
+	c := &circularBuffer{}
+	sink = c // force to heap
+	work(c, msgCount1)
+}
+
+var traceFile string
+var howAllocated = "stack"
+var asBench bool
+
+func flags() {
+	flag.StringVar(&traceFile, "trace", traceFile, "name of trace file to create")
+	flag.StringVar(&howAllocated, "how", howAllocated, "how the buffer is allocated = {stack,heap,global}")
+	flag.BoolVar(&doFluff, "fluff", doFluff, "insert 'fluff' into allocation runs to break up sweeps")
+	flag.BoolVar(&asBench, "bench", asBench, "output in Go benchmark format")
+	flag.Parse()
+}
+
+func bench(b *testing.B) {
+	var c *circularBuffer = &circularBuffer{}
+	delays = make([]time.Duration, 0, msgCount1)
+	// Warmup
+	work(c, msgCount0)
+	c = nil
+	delays = delays[:0]
+
+	total = time.Duration(0)
+	worstElapsed = time.Duration(0)
+	worst = time.Duration(0)
+	worstIndex = 0
+
+	if traceFile != "" {
+		f, err := os.Create(traceFile)
+		if err != nil {
+			fmt.Fprintf(os.Stderr, "Could not create trace file '%s'\n", traceFile)
+			os.Exit(1)
+		}
+		defer f.Close()
+		trace.Start(f)
+		defer trace.Stop()
+	}
+	if b != nil {
+		b.ResetTimer()
+	}
+	allStart = time.Now()
+	c = nil
+
+	switch howAllocated {
+	case "stack":
+		stack()
+	case "heap":
+		heap()
+	case "global":
+		global()
+	default:
+		fmt.Fprintf(os.Stderr, "-how needs to be one of 'heap', 'stack' or 'global, saw '%s' instead\n", howAllocated)
+		os.Exit(1)
+	}
+
+	sort.Slice(delays, func(i, j int) bool { return delays[i] < delays[j] })
+
+	average, median := time.Duration(total.Nanoseconds()/msgCount1), delays[len(delays)/2]
+	p29, p39, p49, p59, p69 := delays[int(0.99*float64(len(delays)))], delays[int(0.999*float64(len(delays)))], delays[int(0.9999*float64(len(delays)))], delays[int(0.99999*float64(len(delays)))], delays[int(0.999999*float64(len(delays)))]
+
+	if asBench {
+		width := int(math.Ceil(math.Log10(float64(1 + worst))))
+		if b != nil {
+			b.ReportMetric(float64(average.Nanoseconds()), "avg-ns")
+			b.ReportMetric(float64(median), "median-ns")
+			b.ReportMetric(float64(p29), "p29-ns")
+			b.ReportMetric(float64(p39), "p39-ns")
+			b.ReportMetric(float64(p49), "p49-ns")
+			b.ReportMetric(float64(p59), "p59-ns")
+			b.ReportMetric(float64(p69), "p69-ns")
+			b.ReportMetric(float64(worst), "worst-ns")
+		} else {
+			fmt.Printf("BenchmarkAverageLatency 1 %[1]*dns\n", width, average)
+			fmt.Printf("BenchmarkMedianLatency  1 %[1]*dns\n", width, median)
+			fmt.Printf("Benchmark99Latency      1 %[1]*dns\n", width, p29)
+			fmt.Printf("Benchmark999Latency     1 %[1]*dns\n", width, p39)
+			fmt.Printf("Benchmark9999Latency    1 %[1]*dns\n", width, p49)
+			fmt.Printf("Benchmark99999Latency   1 %[1]*dns\n", width, p59)
+			fmt.Printf("Benchmark999999Latency  1 %[1]*dns\n", width, p69)
+			fmt.Printf("BenchmarkWorstLatency   1 %[1]*dns\n", width, worst)
+		}
+	} else {
+		fmt.Println("Worst allocation latency:", worst)
+		fmt.Println("Worst allocation index:", worstIndex)
+		fmt.Println("Worst allocation occurs at run elapsed time:", worstElapsed)
+		fmt.Println("Average allocation latency:", average)
+		fmt.Println("Median allocation latency:", median)
+		fmt.Println("99% allocation latency:", p29)
+		fmt.Println("99.9% allocation latency:", p39)
+		fmt.Println("99.99% allocation latency:", p49)
+		fmt.Println("99.999% allocation latency:", p59)
+		fmt.Println("Sizeof(circularBuffer) =", unsafe.Sizeof(*c))
+		fmt.Println("Approximate live memory =", unsafe.Sizeof(*c)+bufferLen*1024)
+	}
+}
diff --git a/gc_latency/latency_test.go b/gc_latency/latency_test.go
new file mode 100644
index 0000000..aba3675
--- /dev/null
+++ b/gc_latency/latency_test.go
@@ -0,0 +1,10 @@
+package main
+
+import (
+	"testing"
+)
+
+func BenchmarkGCLatency(b *testing.B) {
+	asBench = true
+	bench(b)
+}
diff --git a/gc_latency/main.go b/gc_latency/main.go
new file mode 100644
index 0000000..d6b1bd0
--- /dev/null
+++ b/gc_latency/main.go
@@ -0,0 +1,6 @@
+package main
+
+func main() {
+	flags()
+	bench(nil)
+}