Abacus.DoublePrecision.Vector4Tests.TestStaticFn_IsUnit_iii C# (CSharp) Method

        public void TestStaticFn_IsUnit_iii ()
        {
            Double radius = 1;

            Double pi; Maths.Pi(out pi);

            for( Int32 i = 0; i <= 10; ++ i)
            {
                for( Int32 j = 0; j <= 10; ++ j)
                {
                    for( Int32 k = 0; k <= 10; ++ k)
                    {
                        Double theta = 2 * pi * i / 100;
                        Double phi = 2 * pi * j / 100;
                        Double gamma = 2 * pi * k / 100;

                        Double x =
                            Maths.Cos(theta) *
                            Maths.Sin(phi) *
                            Maths.Sin(gamma) * radius;

                        Double y =
                            Maths.Sin(theta) *
                            Maths.Sin(phi) *
                            Maths.Sin(gamma) * radius;

                        Double z =
                            Maths.Cos(phi) *
                            Maths.Sin(gamma) * radius;

                        Double w =
                            Maths.Cos(gamma) * radius;

                        Assert.That(
                            new Vector4(x, y,  z, w).IsUnit(),
                            Is.EqualTo(true));
                    }
                }
            }
        }