src/slices/slices.go

   1 // Copyright 2021 The Go Authors. All rights reserved.
   2 // Use of this source code is governed by a BSD-style
   3 // license that can be found in the LICENSE file.
   4
   5 // Package slices defines various functions useful with slices of any type.
   6 package slices
   7
   8 import (
   9         "cmp"
  10         "unsafe"
  11 )
  12
  13 // Equal reports whether two slices are equal: the same length and all
  14 // elements equal. If the lengths are different, Equal returns false.
  15 // Otherwise, the elements are compared in increasing index order, and the
  16 // comparison stops at the first unequal pair.
  17 // Floating point NaNs are not considered equal.
  18 func Equal[S ~[]E, E comparable](s1, s2 S) bool {
  19         if len(s1) != len(s2) {
  20                 return false
  21         }
  22         for i := range s1 {
  23                 if s1[i] != s2[i] {
  24                         return false
  25                 }
  26         }
  27         return true
  28 }
  29
  30 // EqualFunc reports whether two slices are equal using an equality
  31 // function on each pair of elements. If the lengths are different,
  32 // EqualFunc returns false. Otherwise, the elements are compared in
  33 // increasing index order, and the comparison stops at the first index
  34 // for which eq returns false.
  35 func EqualFunc[S1 ~[]E1, S2 ~[]E2, E1, E2 any](s1 S1, s2 S2, eq func(E1, E2) bool) bool {
  36         if len(s1) != len(s2) {
  37                 return false
  38         }
  39         for i, v1 := range s1 {
  40                 v2 := s2[i]
  41                 if !eq(v1, v2) {
  42                         return false
  43                 }
  44         }
  45         return true
  46 }
  47
  48 // Compare compares the elements of s1 and s2, using [cmp.Compare] on each pair
  49 // of elements. The elements are compared sequentially, starting at index 0,
  50 // until one element is not equal to the other.
  51 // The result of comparing the first non-matching elements is returned.
  52 // If both slices are equal until one of them ends, the shorter slice is
  53 // considered less than the longer one.
  54 // The result is 0 if s1 == s2, -1 if s1 < s2, and +1 if s1 > s2.
  55 func Compare[S ~[]E, E cmp.Ordered](s1, s2 S) int {
  56         for i, v1 := range s1 {
  57                 if i >= len(s2) {
  58                         return +1
  59                 }
  60                 v2 := s2[i]
  61                 if c := cmp.Compare(v1, v2); c != 0 {
  62                         return c
  63                 }
  64         }
  65         if len(s1) < len(s2) {
  66                 return -1
  67         }
  68         return 0
  69 }
  70
  71 // CompareFunc is like [Compare] but uses a custom comparison function on each
  72 // pair of elements.
  73 // The result is the first non-zero result of cmp; if cmp always
  74 // returns 0 the result is 0 if len(s1) == len(s2), -1 if len(s1) < len(s2),
  75 // and +1 if len(s1) > len(s2).
  76 func CompareFunc[S1 ~[]E1, S2 ~[]E2, E1, E2 any](s1 S1, s2 S2, cmp func(E1, E2) int) int {
  77         for i, v1 := range s1 {
  78                 if i >= len(s2) {
  79                         return +1
  80                 }
  81                 v2 := s2[i]
  82                 if c := cmp(v1, v2); c != 0 {
  83                         return c
  84                 }
  85         }
  86         if len(s1) < len(s2) {
  87                 return -1
  88         }
  89         return 0
  90 }
  91
  92 // Index returns the index of the first occurrence of v in s,
  93 // or -1 if not present.
  94 func Index[S ~[]E, E comparable](s S, v E) int {
  95         for i := range s {
  96                 if v == s[i] {
  97                         return i
  98                 }
  99         }
 100         return -1
 101 }
 102
 103 // IndexFunc returns the first index i satisfying f(s[i]),
 104 // or -1 if none do.
 105 func IndexFunc[S ~[]E, E any](s S, f func(E) bool) int {
 106         for i := range s {
 107                 if f(s[i]) {
 108                         return i
 109                 }
 110         }
 111         return -1
 112 }
 113
 114 // Contains reports whether v is present in s.
 115 func Contains[S ~[]E, E comparable](s S, v E) bool {
 116         return Index(s, v) >= 0
 117 }
 118
 119 // ContainsFunc reports whether at least one
 120 // element e of s satisfies f(e).
 121 func ContainsFunc[S ~[]E, E any](s S, f func(E) bool) bool {
 122         return IndexFunc(s, f) >= 0
 123 }
 124
 125 // Insert inserts the values v... into s at index i,
 126 // returning the modified slice.
 127 // The elements at s[i:] are shifted up to make room.
 128 // In the returned slice r, r[i] == v[0],
 129 // and r[i+len(v)] == value originally at r[i].
 130 // Insert panics if i is out of range.
 131 // This function is O(len(s) + len(v)).
 132 func Insert[S ~[]E, E any](s S, i int, v ...E) S {
 133         n := len(s)
 134         m := len(v)
 135         if m == 0 {
 136                 // Panic if i is not in the range [0:n] inclusive.
 137                 // See issue 63913.
 138                 _ = s[:n:n][i:]
 139                 return s
 140         }
 141         if i == n {
 142                 return append(s, v...)
 143         }
 144         if n+m > cap(s) {
 145                 // Use append rather than make so that we bump the size of
 146                 // the slice up to the next storage class.
 147                 // This is what Grow does but we don't call Grow because
 148                 // that might copy the values twice.
 149                 s2 := append(s[:i], make(S, n+m-i)...)
 150                 copy(s2[i:], v)
 151                 copy(s2[i+m:], s[i:])
 152                 return s2
 153         }
 154         s = s[:n+m]
 155
 156         // before:
 157         // s: aaaaaaaabbbbccccccccdddd
 158         //            ^   ^       ^   ^
 159         //            i  i+m      n  n+m
 160         // after:
 161         // s: aaaaaaaavvvvbbbbcccccccc
 162         //            ^   ^       ^   ^
 163         //            i  i+m      n  n+m
 164         //
 165         // a are the values that don't move in s.
 166         // v are the values copied in from v.
 167         // b and c are the values from s that are shifted up in index.
 168         // d are the values that get overwritten, never to be seen again.
 169
 170         if !overlaps(v, s[i+m:]) {
 171                 // Easy case - v does not overlap either the c or d regions.
 172                 // (It might be in some of a or b, or elsewhere entirely.)
 173                 // The data we copy up doesn't write to v at all, so just do it.
 174
 175                 copy(s[i+m:], s[i:])
 176
 177                 // Now we have
 178                 // s: aaaaaaaabbbbbbbbcccccccc
 179                 //            ^   ^       ^   ^
 180                 //            i  i+m      n  n+m
 181                 // Note the b values are duplicated.
 182
 183                 copy(s[i:], v)
 184
 185                 // Now we have
 186                 // s: aaaaaaaavvvvbbbbcccccccc
 187                 //            ^   ^       ^   ^
 188                 //            i  i+m      n  n+m
 189                 // That's the result we want.
 190                 return s
 191         }
 192
 193         // The hard case - v overlaps c or d. We can't just shift up
 194         // the data because we'd move or clobber the values we're trying
 195         // to insert.
 196         // So instead, write v on top of d, then rotate.
 197         copy(s[n:], v)
 198
 199         // Now we have
 200         // s: aaaaaaaabbbbccccccccvvvv
 201         //            ^   ^       ^   ^
 202         //            i  i+m      n  n+m
 203
 204         rotateRight(s[i:], m)
 205
 206         // Now we have
 207         // s: aaaaaaaavvvvbbbbcccccccc
 208         //            ^   ^       ^   ^
 209         //            i  i+m      n  n+m
 210         // That's the result we want.
 211         return s
 212 }
 213
 214 // Delete removes the elements s[i:j] from s, returning the modified slice.
 215 // Delete panics if s[i:j] is not a valid slice of s.
 216 // Delete is O(len(s)-j), so if many items must be deleted, it is better to
 217 // make a single call deleting them all together than to delete one at a time.
 218 // Delete might not modify the elements s[len(s)-(j-i):len(s)]. If those
 219 // elements contain pointers you might consider zeroing those elements so that
 220 // objects they reference can be garbage collected.
 221 func Delete[S ~[]E, E any](s S, i, j int) S {
 222         _ = s[i:j] // bounds check
 223
 224         return append(s[:i], s[j:]...)
 225 }
 226
 227 // DeleteFunc removes any elements from s for which del returns true,
 228 // returning the modified slice.
 229 // When DeleteFunc removes m elements, it might not modify the elements
 230 // s[len(s)-m:len(s)]. If those elements contain pointers you might consider
 231 // zeroing those elements so that objects they reference can be garbage
 232 // collected.
 233 func DeleteFunc[S ~[]E, E any](s S, del func(E) bool) S {
 234         i := IndexFunc(s, del)
 235         if i == -1 {
 236                 return s
 237         }
 238         // Don't start copying elements until we find one to delete.
 239         for j := i + 1; j < len(s); j++ {
 240                 if v := s[j]; !del(v) {
 241                         s[i] = v
 242                         i++
 243                 }
 244         }
 245         return s[:i]
 246 }
 247
 248 // Replace replaces the elements s[i:j] by the given v, and returns the
 249 // modified slice. Replace panics if s[i:j] is not a valid slice of s.
 250 func Replace[S ~[]E, E any](s S, i, j int, v ...E) S {
 251         _ = s[i:j] // verify that i:j is a valid subslice
 252
 253         if i == j {
 254                 return Insert(s, i, v...)
 255         }
 256         if j == len(s) {
 257                 return append(s[:i], v...)
 258         }
 259
 260         tot := len(s[:i]) + len(v) + len(s[j:])
 261         if tot > cap(s) {
 262                 // Too big to fit, allocate and copy over.
 263                 s2 := append(s[:i], make(S, tot-i)...) // See Insert
 264                 copy(s2[i:], v)
 265                 copy(s2[i+len(v):], s[j:])
 266                 return s2
 267         }
 268
 269         r := s[:tot]
 270
 271         if i+len(v) <= j {
 272                 // Easy, as v fits in the deleted portion.
 273                 copy(r[i:], v)
 274                 copy(r[i+len(v):], s[j:])
 275                 return r
 276         }
 277
 278         // We are expanding (v is bigger than j-i).
 279         // The situation is something like this:
 280         // (example has i=4,j=8,len(s)=16,len(v)=6)
 281         // s: aaaaxxxxbbbbbbbbyy
 282         //        ^   ^       ^ ^
 283         //        i   j  len(s) tot
 284         // a: prefix of s
 285         // x: deleted range
 286         // b: more of s
 287         // y: area to expand into
 288
 289         if !overlaps(r[i+len(v):], v) {
 290                 // Easy, as v is not clobbered by the first copy.
 291                 copy(r[i+len(v):], s[j:])
 292                 copy(r[i:], v)
 293                 return r
 294         }
 295
 296         // This is a situation where we don't have a single place to which
 297         // we can copy v. Parts of it need to go to two different places.
 298         // We want to copy the prefix of v into y and the suffix into x, then
 299         // rotate |y| spots to the right.
 300         //
 301         //        v[2:]      v[:2]
 302         //         |           |
 303         // s: aaaavvvvbbbbbbbbvv
 304         //        ^   ^       ^ ^
 305         //        i   j  len(s) tot
 306         //
 307         // If either of those two destinations don't alias v, then we're good.
 308         y := len(v) - (j - i) // length of y portion
 309
 310         if !overlaps(r[i:j], v) {
 311                 copy(r[i:j], v[y:])
 312                 copy(r[len(s):], v[:y])
 313                 rotateRight(r[i:], y)
 314                 return r
 315         }
 316         if !overlaps(r[len(s):], v) {
 317                 copy(r[len(s):], v[:y])
 318                 copy(r[i:j], v[y:])
 319                 rotateRight(r[i:], y)
 320                 return r
 321         }
 322
 323         // Now we know that v overlaps both x and y.
 324         // That means that the entirety of b is *inside* v.
 325         // So we don't need to preserve b at all; instead we
 326         // can copy v first, then copy the b part of v out of
 327         // v to the right destination.
 328         k := startIdx(v, s[j:])
 329         copy(r[i:], v)
 330         copy(r[i+len(v):], r[i+k:])
 331         return r
 332 }
 333
 334 // Clone returns a copy of the slice.
 335 // The elements are copied using assignment, so this is a shallow clone.
 336 func Clone[S ~[]E, E any](s S) S {
 337         // The s[:0:0] preserves nil in case it matters.
 338         return append(s[:0:0], s...)
 339 }
 340
 341 // Compact replaces consecutive runs of equal elements with a single copy.
 342 // This is like the uniq command found on Unix.
 343 // Compact modifies the contents of the slice s and returns the modified slice,
 344 // which may have a smaller length.
 345 // When Compact discards m elements in total, it might not modify the elements
 346 // s[len(s)-m:len(s)]. If those elements contain pointers you might consider
 347 // zeroing those elements so that objects they reference can be garbage collected.
 348 func Compact[S ~[]E, E comparable](s S) S {
 349         if len(s) < 2 {
 350                 return s
 351         }
 352         i := 1
 353         for k := 1; k < len(s); k++ {
 354                 if s[k] != s[k-1] {
 355                         if i != k {
 356                                 s[i] = s[k]
 357                         }
 358                         i++
 359                 }
 360         }
 361         return s[:i]
 362 }
 363
 364 // CompactFunc is like [Compact] but uses an equality function to compare elements.
 365 // For runs of elements that compare equal, CompactFunc keeps the first one.
 366 func CompactFunc[S ~[]E, E any](s S, eq func(E, E) bool) S {
 367         if len(s) < 2 {
 368                 return s
 369         }
 370         i := 1
 371         for k := 1; k < len(s); k++ {
 372                 if !eq(s[k], s[k-1]) {
 373                         if i != k {
 374                                 s[i] = s[k]
 375                         }
 376                         i++
 377                 }
 378         }
 379         return s[:i]
 380 }
 381
 382 // Grow increases the slice's capacity, if necessary, to guarantee space for
 383 // another n elements. After Grow(n), at least n elements can be appended
 384 // to the slice without another allocation. If n is negative or too large to
 385 // allocate the memory, Grow panics.
 386 func Grow[S ~[]E, E any](s S, n int) S {
 387         if n < 0 {
 388                 panic("cannot be negative")
 389         }
 390         if n -= cap(s) - len(s); n > 0 {
 391                 s = append(s[:cap(s)], make([]E, n)...)[:len(s)]
 392         }
 393         return s
 394 }
 395
 396 // Clip removes unused capacity from the slice, returning s[:len(s):len(s)].
 397 func Clip[S ~[]E, E any](s S) S {
 398         return s[:len(s):len(s)]
 399 }
 400
 401 // Rotation algorithm explanation:
 402 //
 403 // rotate left by 2
 404 // start with
 405 //   0123456789
 406 // split up like this
 407 //   01 234567 89
 408 // swap first 2 and last 2
 409 //   89 234567 01
 410 // join first parts
 411 //   89234567 01
 412 // recursively rotate first left part by 2
 413 //   23456789 01
 414 // join at the end
 415 //   2345678901
 416 //
 417 // rotate left by 8
 418 // start with
 419 //   0123456789
 420 // split up like this
 421 //   01 234567 89
 422 // swap first 2 and last 2
 423 //   89 234567 01
 424 // join last parts
 425 //   89 23456701
 426 // recursively rotate second part left by 6
 427 //   89 01234567
 428 // join at the end
 429 //   8901234567
 430
 431 // TODO: There are other rotate algorithms.
 432 // This algorithm has the desirable property that it moves each element exactly twice.
 433 // The triple-reverse algorithm is simpler and more cache friendly, but takes more writes.
 434 // The follow-cycles algorithm can be 1-write but it is not very cache friendly.
 435
 436 // rotateLeft rotates b left by n spaces.
 437 // s_final[i] = s_orig[i+r], wrapping around.
 438 func rotateLeft[E any](s []E, r int) {
 439         for r != 0 && r != len(s) {
 440                 if r*2 <= len(s) {
 441                         swap(s[:r], s[len(s)-r:])
 442                         s = s[:len(s)-r]
 443                 } else {
 444                         swap(s[:len(s)-r], s[r:])
 445                         s, r = s[len(s)-r:], r*2-len(s)
 446                 }
 447         }
 448 }
 449 func rotateRight[E any](s []E, r int) {
 450         rotateLeft(s, len(s)-r)
 451 }
 452
 453 // swap swaps the contents of x and y. x and y must be equal length and disjoint.
 454 func swap[E any](x, y []E) {
 455         for i := 0; i < len(x); i++ {
 456                 x[i], y[i] = y[i], x[i]
 457         }
 458 }
 459
 460 // overlaps reports whether the memory ranges a[0:len(a)] and b[0:len(b)] overlap.
 461 func overlaps[E any](a, b []E) bool {
 462         if len(a) == 0 || len(b) == 0 {
 463                 return false
 464         }
 465         elemSize := unsafe.Sizeof(a[0])
 466         if elemSize == 0 {
 467                 return false
 468         }
 469         // TODO: use a runtime/unsafe facility once one becomes available. See issue 12445.
 470         // Also see crypto/internal/alias/alias.go:AnyOverlap
 471         return uintptr(unsafe.Pointer(&a[0])) <= uintptr(unsafe.Pointer(&b[len(b)-1]))+(elemSize-1) &&
 472                 uintptr(unsafe.Pointer(&b[0])) <= uintptr(unsafe.Pointer(&a[len(a)-1]))+(elemSize-1)
 473 }
 474
 475 // startIdx returns the index in haystack where the needle starts.
 476 // prerequisite: the needle must be aliased entirely inside the haystack.
 477 func startIdx[E any](haystack, needle []E) int {
 478         p := &needle[0]
 479         for i := range haystack {
 480                 if p == &haystack[i] {
 481                         return i
 482                 }
 483         }
 484         // TODO: what if the overlap is by a non-integral number of Es?
 485         panic("needle not found")
 486 }
 487
 488 // Reverse reverses the elements of the slice in place.
 489 func Reverse[S ~[]E, E any](s S) {
 490         for i, j := 0, len(s)-1; i < j; i, j = i+1, j-1 {
 491                 s[i], s[j] = s[j], s[i]
 492         }
 493 }
 494
 495 // Concat returns a new slice concatenating the passed in slices.
 496 func Concat[S ~[]E, E any](slices ...S) S {
 497         size := 0
 498         for _, s := range slices {
 499                 size += len(s)
 500                 if size < 0 {
 501                         panic("len out of range")
 502                 }
 503         }
 504         newslice := Grow[S](nil, size)
 505         for _, s := range slices {
 506                 newslice = append(newslice, s...)
 507         }
 508         return newslice
 509 }