// Warning: Some assembly references could not be resolved automatically. This might lead to incorrect decompilation of some parts,
// for ex. property getter/setter access. To get optimal decompilation results, please manually add the missing references to the list of loaded assemblies.
// EdenAutoMorpherScript, Version=0.0.0.0, Culture=neutral, PublicKeyToken=null
// Eden.AutoMorpher.PcaUtil
using System.Collections.Generic;
using UnityEngine;

public class PcaUtil
{
    public RegionStats ComputeRegionStats(IList<Vector3> points)
    {
        //IL_001c: Unknown result type (might be due to invalid IL or missing references)
        //IL_0021: Unknown result type (might be due to invalid IL or missing references)
        //IL_0027: Unknown result type (might be due to invalid IL or missing references)
        //IL_002b: Unknown result type (might be due to invalid IL or missing references)
        //IL_0030: Unknown result type (might be due to invalid IL or missing references)
        //IL_0035: Unknown result type (might be due to invalid IL or missing references)
        //IL_0041: Unknown result type (might be due to invalid IL or missing references)
        //IL_0044: Unknown result type (might be due to invalid IL or missing references)
        //IL_0049: Unknown result type (might be due to invalid IL or missing references)
        //IL_007e: Unknown result type (might be due to invalid IL or missing references)
        //IL_0083: Unknown result type (might be due to invalid IL or missing references)
        //IL_0084: Unknown result type (might be due to invalid IL or missing references)
        //IL_0089: Unknown result type (might be due to invalid IL or missing references)
        //IL_008c: Unknown result type (might be due to invalid IL or missing references)
        //IL_0093: Unknown result type (might be due to invalid IL or missing references)
        //IL_009f: Unknown result type (might be due to invalid IL or missing references)
        //IL_00a6: Unknown result type (might be due to invalid IL or missing references)
        //IL_00b3: Unknown result type (might be due to invalid IL or missing references)
        //IL_00ba: Unknown result type (might be due to invalid IL or missing references)
        //IL_00c7: Unknown result type (might be due to invalid IL or missing references)
        //IL_00ce: Unknown result type (might be due to invalid IL or missing references)
        //IL_00db: Unknown result type (might be due to invalid IL or missing references)
        //IL_00e2: Unknown result type (might be due to invalid IL or missing references)
        //IL_00ef: Unknown result type (might be due to invalid IL or missing references)
        //IL_00f6: Unknown result type (might be due to invalid IL or missing references)
        //IL_017e: Unknown result type (might be due to invalid IL or missing references)
        //IL_0183: Unknown result type (might be due to invalid IL or missing references)
        //IL_01a2: Unknown result type (might be due to invalid IL or missing references)
        //IL_01a7: Unknown result type (might be due to invalid IL or missing references)
        //IL_01a8: Unknown result type (might be due to invalid IL or missing references)
        //IL_01ad: Unknown result type (might be due to invalid IL or missing references)
        //IL_0209: Unknown result type (might be due to invalid IL or missing references)
        //IL_020a: Unknown result type (might be due to invalid IL or missing references)
        //IL_0211: Unknown result type (might be due to invalid IL or missing references)
        //IL_0213: Unknown result type (might be due to invalid IL or missing references)
        //IL_01ca: Unknown result type (might be due to invalid IL or missing references)
        //IL_01cb: Unknown result type (might be due to invalid IL or missing references)
        //IL_01cf: Unknown result type (might be due to invalid IL or missing references)
        //IL_01d4: Unknown result type (might be due to invalid IL or missing references)
        //IL_01d9: Unknown result type (might be due to invalid IL or missing references)
        //IL_01de: Unknown result type (might be due to invalid IL or missing references)
        //IL_01e3: Unknown result type (might be due to invalid IL or missing references)
        //IL_01e5: Unknown result type (might be due to invalid IL or missing references)
        //IL_01ea: Unknown result type (might be due to invalid IL or missing references)
        RegionStats result = default(RegionStats);
        if (points == null || points.Count == 0)
        {
            return result;
        }
        int count = points.Count;
        Vector3 val = Vector3.zero;
        for (int i = 0; i < count; i++)
        {
            val += points[i];
        }
        val /= (float)count;
        float num = 0f;
        float num2 = 0f;
        float num3 = 0f;
        float num4 = 0f;
        float num5 = 0f;
        float num6 = 0f;
        for (int j = 0; j < count; j++)
        {
            Vector3 val2 = points[j] - val;
            num += val2.x * val2.x;
            num2 += val2.x * val2.y;
            num3 += val2.x * val2.z;
            num4 += val2.y * val2.y;
            num5 += val2.y * val2.z;
            num6 += val2.z * val2.z;
        }
        float num7 = 1f / (float)count;
        num *= num7;
        num2 *= num7;
        num3 *= num7;
        num4 *= num7;
        num5 *= num7;
        num6 *= num7;
        this.JacobiEigenDecomposition3x3(num, num2, num3, num4, num5, num6, out var eigenValues, out var eigenVectors);
        int num8 = 0;
        if (eigenValues[1] > eigenValues[num8])
        {
            num8 = 1;
        }
        if (eigenValues[2] > eigenValues[num8])
        {
            num8 = 2;
        }
        Vector3 normalized = eigenVectors[num8].normalized;
        float num9 = float.PositiveInfinity;
        float num10 = float.NegativeInfinity;
        float num11 = 0f;
        for (int k = 0; k < count; k++)
        {
            float num12 = Vector3.Dot(points[k] - val, normalized);
            if (num12 < num9)
            {
                num9 = num12;
            }
            if (num12 > num10)
            {
                num10 = num12;
            }
            Vector3 val3 = val + normalized * num12;
            Vector3 val4 = points[k] - val3;
            float magnitude = val4.magnitude;
            num11 += magnitude;
        }
        result.center = val;
        result.principalAxis = normalized;
        result.length = num10 - num9;
        result.avgRadius = num11 / (float)count;
        return result;
    }

    private void JacobiEigenDecomposition3x3(float c00, float c01, float c02, float c11, float c12, float c22, out float[] eigenValues, out Vector3[] eigenVectors)
    {
        //IL_032d: Unknown result type (might be due to invalid IL or missing references)
        //IL_0332: Unknown result type (might be due to invalid IL or missing references)
        //IL_0353: Unknown result type (might be due to invalid IL or missing references)
        //IL_0358: Unknown result type (might be due to invalid IL or missing references)
        //IL_0379: Unknown result type (might be due to invalid IL or missing references)
        //IL_037e: Unknown result type (might be due to invalid IL or missing references)
        float[,] array = new float[3, 3]
        {
            { c00, c01, c02 },
            { c01, c11, c12 },
            { c02, c12, c22 }
        };
        float[,] array2 = new float[3, 3]
        {
            { 1f, 0f, 0f },
            { 0f, 1f, 0f },
            { 0f, 0f, 1f }
        };
        for (int i = 0; i < 32; i++)
        {
            int num = 0;
            int num2 = 1;
            float num3 = Mathf.Abs(array[0, 1]);
            float num4 = Mathf.Abs(array[0, 2]);
            if (num4 > num3)
            {
                num3 = num4;
                num = 0;
                num2 = 2;
            }
            float num5 = Mathf.Abs(array[1, 2]);
            if (num5 > num3)
            {
                num3 = num5;
                num = 1;
                num2 = 2;
            }
            if (num3 < 1E-10f)
            {
                break;
            }
            float num6 = array[num, num];
            float num7 = array[num2, num2];
            float num8 = array[num, num2];
            float num9 = 0.5f * Mathf.Atan2(2f * num8, num7 - num6);
            float num10 = Mathf.Cos(num9);
            float num11 = Mathf.Sin(num9);
            for (int j = 0; j < 3; j++)
            {
                if (j != num && j != num2)
                {
                    float num12 = array[j, num];
                    float num13 = array[j, num2];
                    array[j, num] = num10 * num12 - num11 * num13;
                    array[num, j] = array[j, num];
                    array[j, num2] = num11 * num12 + num10 * num13;
                    array[num2, j] = array[j, num2];
                }
            }
            float num14 = num10 * num10 * num6 - 2f * num11 * num10 * num8 + num11 * num11 * num7;
            float num15 = num11 * num11 * num6 + 2f * num11 * num10 * num8 + num10 * num10 * num7;
            array[num, num] = num14;
            array[num2, num2] = num15;
            array[num, num2] = 0f;
            array[num2, num] = 0f;
            for (int k = 0; k < 3; k++)
            {
                float num16 = array2[k, num];
                float num17 = array2[k, num2];
                array2[k, num] = num10 * num16 - num11 * num17;
                array2[k, num2] = num11 * num16 + num10 * num17;
            }
        }
        eigenValues = new float[3];
        eigenVectors = (Vector3[])(object)new Vector3[3];
        eigenValues[0] = array[0, 0];
        eigenValues[1] = array[1, 1];
        eigenValues[2] = array[2, 2];
        eigenVectors[0] = new Vector3(array2[0, 0], array2[1, 0], array2[2, 0]);
        eigenVectors[1] = new Vector3(array2[0, 1], array2[1, 1], array2[2, 1]);
        eigenVectors[2] = new Vector3(array2[0, 2], array2[1, 2], array2[2, 2]);
    }
}