//-----------------------------------------------------------------------
// -----------------------------------------------------------------------
// file="MatrixTest"
// Project: StgSharp
// AuthorGroup: Nitload
// Copyright (c) Nitload. All rights reserved.
//     
// 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.
//     
// -----------------------------------------------------------------------
// -----------------------------------------------------------------------
using StgSharp.Mathematics.Memory;
using StgSharp.Mathematics.Numeric;
using System;
using System.Collections.Generic;
using System.Linq;
using System.Numerics;
using System.Text;
using System.Threading.Tasks;

using T = float;

namespace StgSharp
{
    internal static partial class Program
    {

        private const long size = 1024L * 1024L * 1024L * 4L;

        private static int sideLength = 1000;

        internal static unsafe void TestMatrix(
        )
        {
            // TestRender();
            /**/
            HybridLayerSegregatedFitAllocator hlsf;
            unchecked {
                hlsf = new((nuint)size, 32, size);
            }
            Matrix<T> m = Matrix<T>.Create(sideLength, sideLength, MatrixLayout.DenseRectangle,
                                           hlsf);

            // MatrixCompute.Fill(m, 1.0f);
            Matrix<T> n = Matrix<T>.Create(sideLength, sideLength, MatrixLayout.DenseRectangle,
                                           hlsf);

            // MatrixCompute.Fill(n, 1.0f);
            Matrix<T> p = Matrix<T>.Create(sideLength, sideLength, MatrixLayout.DenseRectangle,
                                           hlsf);

            // MatrixCompute.Fill(p, 0.0f);

            /*
            for (int j = 0; j < 20; j++)
            {
                for (int k = 0; k < 10; k++) {
                    Console.Write($"{m[j, k]:0.0E+00}\t");
                }
                Console.WriteLine();
            }

            /**/
            Console.WriteLine("===================================================");

            // MatrixCompute.Mul(n, m, p);
            /**/
            for (int i = 0; i < 20; i++)
            {
                if (i % 4 == 0) {
                    Console.WriteLine();
                }
                for (int j = 0; j < 20; j++)
                {
                    string extra = j % 4 == 3 ? "   " : ",";
                    T value = p[j, i];

                    // value = T.Abs(value) < 1e-6 ? 0 : value;

                    // Console.Write($"{p.GetElementOffset(j,i)}{extra}");

                    Console.Write($"{value:000}{extra}");
                }
                Console.WriteLine();
            }
            Console.WriteLine("===================================================");

            /**/
        }

    }
}