protected override void AdaptMesh(int TimestepNo, out GridCommons newGrid, out GridCorrelation old2NewGrid)
        {
            if (TimestepNo > 3 && TimestepNo % 3 != 0)
            {
                //{
                newGrid     = null;
                old2NewGrid = null;
                return;
            }

            // Check grid changes
            // ==================
            bool                     AnyChange;
            List <int>               CellsToRefineList;
            List <int[]>             Coarsening;
            GridRefinementController gridRefinementController = new GridRefinementController((GridData)this.GridData, LsTrk.Regions.GetCutCellMask());

            if (MPISize > 1)
            {
                List <Tuple <int, BitArray> > cellMaskRefinementLevel = GetCellMaskWithRefinementLevels();
                AnyChange = gridRefinementController.ComputeGridChange(cellMaskRefinementLevel, out CellsToRefineList, out Coarsening);
            }
            else
            {
                AnyChange = gridRefinementController.ComputeGridChange(LevelIndicator, out CellsToRefineList, out Coarsening);
            }
            int NoOfCellsToRefine  = 0;
            int NoOfCellsToCoarsen = 0;

            if (AnyChange)
            {
                int[] glb = (new int[] {
                    CellsToRefineList.Count,
                    Coarsening.Sum(L => L.Length),
                }).MPISum();
                NoOfCellsToRefine  = glb[0];
                NoOfCellsToCoarsen = glb[1];
            }
            int oldJ = this.GridData.CellPartitioning.TotalLength;

            // Update Grid
            // ===========

            if (AnyChange)
            {
                Console.WriteLine("       Refining " + NoOfCellsToRefine + " of " + oldJ + " cells");
                Console.WriteLine("       Coarsening " + NoOfCellsToCoarsen + " of " + oldJ + " cells");


                newGrid = ((GridData)this.GridData).Adapt(CellsToRefineList, Coarsening, out old2NewGrid);
            }
            else
            {
                newGrid     = null;
                old2NewGrid = null;
            }
        }
Esempio n. 2
0
        /// <summary>
        /// control of mesh adaptation
        /// </summary>
        protected override void AdaptMesh(int TimestepNo, out GridCommons newGrid, out GridCorrelation old2NewGrid)
        {
            if (this.Control.AdaptiveMeshRefinement && TimestepNo > 1)
            {
                int oldJ = this.GridData.CellPartitioning.TotalLength;

                double LocNormPow2         = this.ResiualKP1.CoordinateVector.L2NormPow2(); // norm of residual on this processor
                double TotNormPow2         = LocNormPow2.MPISum();                          //                          norm of residual over all processors
                double MeanNormPow2PerCell = TotNormPow2 / oldJ;                            //                    mean norm per cell

                double maxSoFar = 0;
                int    jMax     = -1;
                for (int j = 0; j < oldJ; j++)
                {
                    double CellNorm = Error.Coordinates.GetRow(j).L2NormPow2();

                    if (CellNorm > maxSoFar)
                    {
                        jMax     = j;
                        maxSoFar = CellNorm;
                    }
                }


                int LevelIndicator(int j, int CurrentLevel)
                {
                    double CellNorm = this.ResiualKP1.Coordinates.GetRow(j).L2NormPow2();

                    if (j == jMax)
                    {
                        return(CurrentLevel + 1);
                    }
                    else
                    {
                        return(CurrentLevel);
                    }
                }

                bool AnyChange          = GridRefinementController.ComputeGridChange((GridData)(this.GridData), null, LevelIndicator, out List <int> CellsToRefineList, out List <int[]> Coarsening);
                int  NoOfCellsToRefine  = 0;
                int  NoOfCellsToCoarsen = 0;
                if (AnyChange)
                {
                    int[] glb = (new int[] {
                        CellsToRefineList.Count,
                        Coarsening.Sum(L => L.Length),
                    }).MPISum();

                    NoOfCellsToRefine  = glb[0];
                    NoOfCellsToCoarsen = glb[1];
                }

                // Update Grid
                // ===========

                if (AnyChange)
                {
                    Console.WriteLine("       Refining " + NoOfCellsToRefine + " of " + oldJ + " cells");
                    Console.WriteLine("       Coarsening " + NoOfCellsToCoarsen + " of " + oldJ + " cells");

                    newGrid = ((GridData)(this.GridData)).Adapt(CellsToRefineList, Coarsening, out old2NewGrid);
                }
                else
                {
                    newGrid     = null;
                    old2NewGrid = null;
                }
            }
            else
            {
                newGrid     = null;
                old2NewGrid = null;
            }
        }