aboutsummaryrefslogtreecommitdiff
path: root/sketch.go
diff options
context:
space:
mode:
authorn-peugnet <n.peugnet@free.fr>2021-09-09 12:09:18 +0200
committern-peugnet <n.peugnet@free.fr>2021-09-09 12:10:45 +0200
commit8a03c46bf24b5a1fa1d2080ac4f763532db01bbe (patch)
tree069554f3e1e3e235a22d13dbb4a4a555b2d6e0d6 /sketch.go
parentf061a7031181ef53d034c46b696156c143451cce (diff)
downloaddna-backup-8a03c46bf24b5a1fa1d2080ac4f763532db01bbe.tar.gz
dna-backup-8a03c46bf24b5a1fa1d2080ac4f763532db01bbe.zip
export sketch in its own package
so that tests can be cached and to make sure it is independant of the rest of the code also move tests in testdata as this folder is ignored by go test by default
Diffstat (limited to 'sketch.go')
-rw-r--r--sketch.go80
1 files changed, 0 insertions, 80 deletions
diff --git a/sketch.go b/sketch.go
deleted file mode 100644
index dca813d..0000000
--- a/sketch.go
+++ /dev/null
@@ -1,80 +0,0 @@
-package main
-
-import (
- "bytes"
- "encoding/binary"
- "io"
- "log"
- "sync"
-
- "github.com/chmduquesne/rollinghash/rabinkarp64"
-)
-
-type Sketch []uint64
-
-type ReadByteReader interface {
- io.Reader
- io.ByteReader
-}
-
-const fBytes = 8
-
-// SketchChunk produces a sketch for a chunk based on wSize: the window size,
-// sfCount: the number of super-features, and fCount: the number of feature
-// per super-feature
-func SketchChunk(chunk Chunk, pol rabinkarp64.Pol, chunkSize int, wSize int, sfCount int, fCount int) (Sketch, error) {
- var wg sync.WaitGroup
- var fSize = FeatureSize(chunkSize, sfCount, fCount)
- superfeatures := make([]uint64, 0, sfCount)
- features := make([]uint64, 0, fCount*sfCount)
- sfBuff := make([]byte, fBytes*fCount)
- r := chunk.Reader()
- for f := 0; f < chunk.Len()/fSize; f++ {
- var fBuff bytes.Buffer
- n, err := io.CopyN(&fBuff, r, int64(fSize))
- if err != nil {
- log.Println(n, err)
- }
- features = append(features, 0)
- wg.Add(1)
- go calcFeature(&wg, pol, &fBuff, wSize, fSize, &features[f])
- }
- hasher := rabinkarp64.NewFromPol(pol)
- wg.Wait()
- for sf := 0; sf < len(features)/fCount; sf++ {
- for i := 0; i < fCount; i++ {
- binary.LittleEndian.PutUint64(sfBuff[i*fBytes:(i+1)*fBytes], features[i+sf*fCount])
- }
- hasher.Reset()
- hasher.Write(sfBuff)
- superfeatures = append(superfeatures, hasher.Sum64())
- }
- return superfeatures, nil
-}
-
-func calcFeature(wg *sync.WaitGroup, p rabinkarp64.Pol, r ReadByteReader, wSize int, fSize int, result *uint64) {
- defer wg.Done()
- hasher := rabinkarp64.NewFromPol(p)
- n, err := io.CopyN(hasher, r, int64(wSize))
- if err != nil {
- log.Println(n, err)
- }
- max := hasher.Sum64()
- for w := 0; w < fSize-wSize; w++ {
- b, _ := r.ReadByte()
- hasher.Roll(b)
- h := hasher.Sum64()
- if h > max {
- max = h
- }
- }
- *result = max
-}
-
-func SuperFeatureSize(chunkSize int, sfCount int, fCount int) int {
- return FeatureSize(chunkSize, sfCount, fCount) * sfCount
-}
-
-func FeatureSize(chunkSize int, sfCount int, fCount int) int {
- return chunkSize / (sfCount * fCount)
-}