Create TensorBase4x4x4.cs

src/Mathematics.NET/LinearAlgebra/TensorBase4x4x4.cs

@@ -0,0 +1,131 @@
+// <copyright file="TensorBase4x4x4.cs" company="Mathematics.NET">
+// Mathematics.NET
+// https://github.com/HamletTanyavong/Mathematics.NET
+//
+// MIT License
+//
+// Copyright (c) 2023 Hamlet Tanyavong
+//
+// Permission is hereby granted, free of charge, to any person obtaining a copy
+// of this software and associated documentation files (the "Software"), to deal
+// in the Software without restriction, including without limitation the rights
+// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+// copies of the Software, and to permit persons to whom the Software is
+// furnished to do so, subject to the following conditions:
+//
+// The above copyright notice and this permission notice shall be included in all
+// copies or substantial portions of the Software.
+//
+// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+// SOFTWARE.
+// </copyright>
+
+using System.Diagnostics.CodeAnalysis;
+using System.Runtime.CompilerServices;
+using System.Runtime.InteropServices;
+using Mathematics.NET.LinearAlgebra.Abstractions;
+
+namespace Mathematics.NET.LinearAlgebra;
+
+[StructLayout(LayoutKind.Sequential)]
+public struct TensorBase4x4x4<T>
+    : IThreeDimensionalArrayRepresentable<TensorBase4x4x4<T>, T>,
+      IHypercubicArray
+    where T : IComplex<T>
+{
+    public const int Components = 64;
+    public const int E1Components = 16;
+    public const int E2Components = 16;
+    public const int E3Components = 16;
+
+    public Matrix4x4<T> X1;
+    public Matrix4x4<T> X2;
+    public Matrix4x4<T> X3;
+    public Matrix4x4<T> X4;
+
+    //
+    // Constants
+    //
+
+    static int IArrayRepresentable<TensorBase4x4x4<T>, T>.Components => Components;
+    static int IThreeDimensionalArrayRepresentable<TensorBase4x4x4<T>, T>.E1Components => E1Components;
+    static int IThreeDimensionalArrayRepresentable<TensorBase4x4x4<T>, T>.E2Components => E2Components;
+    static int IThreeDimensionalArrayRepresentable<TensorBase4x4x4<T>, T>.E3Components => E3Components;
+
+    //
+    // Indexer
+    //
+
+    public T this[int i, int j, int k]
+    {
+        [MethodImpl(MethodImplOptions.AggressiveInlining)]
+        readonly get
+        {
+            if ((uint)i >= 4)
+            {
+                throw new IndexOutOfRangeException();
+            }
+            return Unsafe.Add(ref Unsafe.AsRef(in X1), i)[j, k];
+        }
+
+        [MethodImpl(MethodImplOptions.AggressiveInlining)]
+        set
+        {
+            if ((uint)i >= 4)
+            {
+                throw new NotImplementedException();
+            }
+            Unsafe.Add(ref X1, i)[j, k] = value;
+        }
+    }
+
+    //
+    // Equality
+    //
+
+    public static bool operator ==(TensorBase4x4x4<T> left, TensorBase4x4x4<T> right)
+        => left.X1 == right.X1 && left.X2 == right.X2 && left.X3 == right.X3 && left.X4 == right.X4;
+
+    public static bool operator !=(TensorBase4x4x4<T> left, TensorBase4x4x4<T> right)
+        => left.X1 != right.X1 || left.X2 != right.X2 || left.X3 != right.X3 || left.X4 != right.X4;
+
+    public override bool Equals([NotNullWhen(true)] object? obj) => obj is TensorBase4x4x4<T> other && Equals(other);
+
+    public bool Equals(TensorBase4x4x4<T> value)
+        => X1.Equals(value.X1) && X2.Equals(value.X2) && X3.Equals(value.X3) && X4.Equals(value.X4);
+
+    public override int GetHashCode() => HashCode.Combine(X1, X2, X3, X4);
+
+    //
+    // Formatting
+    //
+
+    public string ToString(string? format, IFormatProvider? provider) => ToArray3D().ToDisplayString(format, provider);
+
+    //
+    // Methods
+    //
+
+    public T[,,] ToArray3D()
+    {
+        T[,,] result = new T[4, 4, 4];
+
+        for (int i = 0; i < 4; i++)
+        {
+            for (int j = 0; j < 4; j++)
+            {
+                for (int k = 0; k < 4; k++)
+                {
+                    result[i, j, k] = this[i, j, k];
+                }
+            }
+        }
+
+        return result;
+    }
+}