/// <summary>
        /// Subdivides an incoming Floor according to x- and y-axis distances and grid rotation.
        /// </summary>
        /// <param name="model">The input model.</param>
        /// <param name="input">The arguments to the execution.</param>
        /// <returns>A SubdivideSlabOutputs instance containing computed results and the model with any new elements.</returns>
        public static SubdivideSlabOutputs Execute(Dictionary <string, Model> inputModels, SubdivideSlabInputs input)
        {
            var allFloors = new List <Floor>();

            inputModels.TryGetValue("Floors", out var flrModel);
            if (flrModel == null || flrModel.AllElementsOfType <Floor>().Count() == 0)
            {
                throw new ArgumentException("No Floors found.");
            }
            allFloors.AddRange(flrModel.AllElementsOfType <Floor>());
            List <SlabSubdivision> subdivisions = new List <SlabSubdivision>();

            for (int i = 0; i < allFloors.Count; i++)
            {
                Floor floor      = allFloors[i];
                var   floorId    = StringExtensions.NumberToString(i);
                var   profile    = floor.Profile;
                var   perimeter  = profile.Perimeter;
                var   voids      = profile.Voids;
                var   openings   = floor.Openings.Select(o => o.Profile);
                var   elevation  = floor.Elevation;
                var   boundaries = new List <Polygon>();
                boundaries.Add(perimeter);
                if (voids != null)
                {
                    boundaries.AddRange(voids);
                }
                if (openings != null)
                {
                    boundaries.AddRange(openings.Select(o => o.Perimeter));
                }
                Transform transform = null;
                if (input.AlignToLongestEdge)
                {
                    var longestEdge = perimeter.Segments().OrderByDescending(p => p.Length()).First();
                    var xAxis       = (longestEdge.End - longestEdge.Start).Unitized();
                    if (perimeter.IsClockWise())
                    {
                        xAxis = xAxis * -1;
                    }
                    transform = new Transform(Vector3.Origin, xAxis, Vector3.ZAxis, 0);
                    // if (perimeter.IsClockWise())
                    // {
                    //     transform.Invert();
                    // }
                }
                var grid = new Grid2d(boundaries, transform);
                if (input.SubdivideAtVoidCorners && voids != null && voids.Count > 0)
                {
                    foreach (var voidCrv in voids)
                    {
                        grid.SplitAtPoints(voidCrv.Vertices);
                    }
                    foreach (var cell in grid.CellsFlat)
                    {
                        cell.U.DivideByApproximateLength(input.Length, EvenDivisionMode.RoundUp);
                        cell.V.DivideByApproximateLength(input.Width, EvenDivisionMode.RoundUp);
                    }
                }
                else
                {
                    grid.U.DivideByApproximateLength(input.Length, EvenDivisionMode.RoundUp);
                    grid.V.DivideByApproximateLength(input.Width, EvenDivisionMode.RoundUp);
                }

                var cells = grid.GetCells();

                for (int i1 = 0; i1 < cells.Count; i1++)
                {
                    var    id        = $"{floorId}-{i1:000}";
                    Grid2d cell      = cells[i1];
                    var    cellCrvs  = cell.GetTrimmedCellGeometry();
                    var    isTrimmed = cell.IsTrimmed();
                    if (cellCrvs != null && cellCrvs.Length > 0)
                    {
                        subdivisions.Add(CreateSlabSubdivision(id, cellCrvs, floor, isTrimmed));
                    }
                }
            }

            var output = new SubdivideSlabOutputs(subdivisions.Count);

            output.Model.AddElements(subdivisions);
            foreach (var subdiv in subdivisions)
            {
                var thicknessXform = new Transform(0, 0, subdiv.Depth);
                var profile        = subdiv.Transform.OfProfile(subdiv.Profile);
                profile = thicknessXform.OfProfile(profile);
                output.Model.AddElement(new ModelCurve(profile.Perimeter));
                if (profile.Voids != null && profile.Voids.Count > 0)
                {
                    output.Model.AddElements(profile.Voids.Select(v => new ModelCurve(v)));
                }
            }
            return(output);
        }