using System; using System.Runtime.CompilerServices; using System.Runtime.InteropServices; using System.Threading; using TValue = System.Int32; namespace DCFApixels.DragonECS.Graphs.Internal { internal sealed unsafe class SparseMatrix { public const int MIN_CAPACITY_BITS_OFFSET = 4; public const int MIN_CAPACITY = 1 << MIN_CAPACITY_BITS_OFFSET; private const int CHAIN_LENGTH_THRESHOLD = 5; private const float CHAIN_LENGTH_THRESHOLD_CAPCITY_THRESHOLD = 0.7f; private Basket* _buckets; private Entry* _entries; private int _capacity; private int _count_Threshold; private int _count; private int _freeList; private int _freeCount; private int _modBitMask; #region Properties public int Count { [MethodImpl(MethodImplOptions.AggressiveInlining)] get { return _count; } } public int Capacity { [MethodImpl(MethodImplOptions.AggressiveInlining)] get { return _capacity; } } #endregion #region Constructors [MethodImpl(MethodImplOptions.AggressiveInlining)] public SparseMatrix(int minCapacity = MIN_CAPACITY) { minCapacity = NormalizeCapacity(minCapacity); //_buckets = new UnsafeArray(minCapacity); //_entries = new UnsafeArray(minCapacity, true); _buckets = UnmanagedArrayUtility.New(minCapacity); _entries = UnmanagedArrayUtility.NewAndInit(minCapacity); for (int i = 0; i < minCapacity; i++) { _buckets[i] = Basket.Empty; } _modBitMask = (minCapacity - 1) & 0x7FFFFFFF; _count = 0; _freeList = 0; _freeCount = 0; SetCapacity(minCapacity); } ~SparseMatrix() { UnmanagedArrayUtility.Free(_buckets); UnmanagedArrayUtility.Free(_entries); } #endregion #region Add/TryAdd/Set [MethodImpl(MethodImplOptions.AggressiveInlining)] public void Add(int x, int y, TValue value) { unchecked { long key = KeyUtility.FromXY(x, y); #if DEBUG if (FindEntry(key) >= 0) { Throw.ArgumentException("Has(x, y) is true"); } #endif int targetBucket = (int)key & _modBitMask; AddInternal(key, targetBucket, value); } } [MethodImpl(MethodImplOptions.AggressiveInlining)] public bool TryAdd(int x, int y, TValue value) { long key = KeyUtility.FromXY(x, y); if (FindEntry(key) >= 0) { return false; } int targetBucket = (int)key & _modBitMask; AddInternal(key, targetBucket, value); return true; } [MethodImpl(MethodImplOptions.AggressiveInlining)] public void Set(int x, int y, TValue value) { long key = KeyUtility.FromXY(x, y); int targetBucket = (int)key & _modBitMask; for (int i = _buckets[targetBucket].index; i >= 0; i = _entries[i].next) { if (_entries[i].key == key) { _entries[i].value = value; return; } } AddInternal(key, targetBucket, value); } [MethodImpl(MethodImplOptions.AggressiveInlining)] private void AddInternal(long key, int targetBucket, int value) { unchecked { int index; if (_freeCount == 0) { if (_count >= _capacity) { Resize(); // обновляем под новое значение _modBitMask targetBucket = (int)key & _modBitMask; } //index = Interlocked.Increment(ref _count); index = _count++; } else { index = _freeList; _freeList = _entries[index].next; _freeCount--; } #if DEBUG if (_freeCount < 0) { Throw.UndefinedException(); } #endif ref Basket basket = ref _buckets[targetBucket]; ref Entry entry = ref _entries[index]; entry.next = basket.index; entry.key = key; entry.value = value; Interlocked.Increment(ref basket.count); //basket.count++; basket.index = index; if (basket.count >= CHAIN_LENGTH_THRESHOLD && _count > _count_Threshold) //_count / _capacity >= CHAIN_LENGTH_THRESHOLD_CAPCITY_THRESHOLD) { Resize(); } } } #endregion #region FindEntry/Has [MethodImpl(MethodImplOptions.AggressiveInlining)] private int FindEntry(int x, int y) { long key = KeyUtility.FromXY(x, y); for (int i = _buckets[key & _modBitMask].index; i >= 0; i = _entries[i].next) { if (_entries[i].key == key) { return i; } } return -1; } [MethodImpl(MethodImplOptions.AggressiveInlining)] private int FindEntry(long key) { for (int i = _buckets[key & _modBitMask].index; i >= 0; i = _entries[i].next) { if (_entries[i].key == key) { return i; } } return -1; } [MethodImpl(MethodImplOptions.AggressiveInlining)] public bool HasKey(int x, int y) { return FindEntry(x, y) >= 0; } #endregion #region GetValue/TryGetValue [MethodImpl(MethodImplOptions.AggressiveInlining)] public TValue GetValue(int x, int y) { int index = FindEntry(x, y); #if DEBUG if (index < 0) { Throw.KeyNotFound(); } #endif return _entries[index].value; } [MethodImpl(MethodImplOptions.AggressiveInlining)] public bool TryGetValue(int x, int y, out TValue value) { int index = FindEntry(x, y); if (index < 0) { value = default; return false; } value = _entries[index].value; return true; } #endregion #region TryDel [MethodImpl(MethodImplOptions.AggressiveInlining)] public bool TryDel(int x, int y) { long key = KeyUtility.FromXY(x, y); int targetBucket = (int)key & _modBitMask; ref Basket basket = ref _buckets[targetBucket]; int last = -1; for (int i = basket.index; i >= 0; last = i, i = _entries[i].next) { if (_entries[i].key == key) { if (last < 0) { basket.index = _entries[i].next; } else { _entries[last].next = _entries[i].next; } _entries[i].next = _freeList; _entries[i].key = default; //Key.Null; _entries[i].value = default; _freeList = i; _freeCount++; basket.count--; return true; } } return false; } #endregion #region Clear public void Clear() { if (_count > 0) { for (int i = 0; i < _capacity; i++) { _buckets[i] = Basket.Empty; } for (int i = 0; i < _capacity; i++) { _entries[i] = default; } _count = 0; } } #endregion #region Resize [MethodImpl(MethodImplOptions.NoInlining)] private void Resize() { int newSize = _capacity << 1; Console.WriteLine($"Resize {newSize}"); _modBitMask = (newSize - 1) & 0x7FFFFFFF; //newBuckets create and ini Basket* newBuckets = UnmanagedArrayUtility.New(newSize); for (int i = 0; i < newSize; i++) { newBuckets[i] = Basket.Empty; } //END newBuckets create and ini Entry* newEntries = UnmanagedArrayUtility.ResizeAndInit(_entries, _capacity, newSize); for (int i = 0; i < _count; i++) { if (newEntries[i].key >= 0) { ref Entry entry = ref newEntries[i]; ref Basket basket = ref newBuckets[entry.key & _modBitMask]; entry.next = basket.index; basket.index = i; basket.count++; } } UnmanagedArrayUtility.Free(_buckets); _buckets = newBuckets; _entries = newEntries; SetCapacity(newSize); } private void SetCapacity(int newSize) { _capacity = newSize; _count_Threshold = (int)(_capacity * CHAIN_LENGTH_THRESHOLD_CAPCITY_THRESHOLD); } [MethodImpl(MethodImplOptions.AggressiveInlining)] private static int NormalizeCapacity(int capacity) { int result = MIN_CAPACITY; while (result < capacity) { result <<= 1; } return result; } #endregion #region Utils [StructLayout(LayoutKind.Sequential, Pack = 4)] private struct Entry { public int next; // Index of next entry, -1 if last public long key; public TValue value; public override string ToString() { return key == 0 ? "NULL" : $"{key} {value}"; } } [StructLayout(LayoutKind.Sequential, Pack = 4, Size = 8)] public struct Basket { public static readonly Basket Empty = new Basket(-1, 0); public int index; public int count; [MethodImpl(MethodImplOptions.AggressiveInlining)] public Basket(int index, int count) { this.index = index; this.count = count; } public override string ToString() { return index < 0 ? "NULL" : $"{index} {count}"; } } public static class KeyUtility { [MethodImpl(MethodImplOptions.AggressiveInlining)] public static long FromXY(int x, int y) { unchecked { //long result = ((long)x << 32) | (long)(x ^ y ^ Mixing(y)); //Console.WriteLine($"key {x}-{y}-{result}"); //return result; //return ((long)x << 32) | (long)(x ^ y ^ Mixing(y)); return ((long)x << 32) | ((long)(x << 2) ^ y); } } [MethodImpl(MethodImplOptions.AggressiveInlining)] private static int Mixing(int v) { unchecked { v *= 3571; //v ^= v << 13; //v ^= v >> 17; v ^= v >> 8; return v; } } } //[StructLayout(LayoutKind.Explicit, Pack = 4, Size = 8)] //public readonly struct Key : IEquatable //{ // public static readonly Key Null = new Key(-1, 0); // // [FieldOffset(0)] // public readonly long Full; // [FieldOffset(0)] // public readonly int X; // [FieldOffset(4)] // public readonly int YHash; // // [MethodImpl(MethodImplOptions.AggressiveInlining)] // private Key(int x, int yHash) : this() // { // this.X = x; // this.YHash = yHash; // } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // public static unsafe Key FromXY(int x, int y) // { // unchecked // { // return new Key(x, x ^ y ^ Mixing(y)); // } // } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // private static int Mixing(int x) // { // unchecked // { // x *= 3571; // x ^= x << 13; // x ^= x >> 17; // return x; // } // } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // internal static bool EqualsInFind(Key a, Key b) { return a.X == b.X; } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // public static bool operator ==(Key a, Key b) { return a.X == b.X && a.YHash == b.YHash; } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // public static bool operator !=(Key a, Key b) { return a.X != b.X || a.YHash != b.YHash; } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // public override int GetHashCode() { return YHash; } // [MethodImpl(MethodImplOptions.AggressiveInlining)] // public bool Equals(Key other) { return this == other; } // public override bool Equals(object obj) { return obj is Key && Equals((Key)obj); } // public override string ToString() { return $"({X}, {YHash})"; } //} #endregion } }