public static int GetIdxOfFirstBlockWith(MultigridMapping map, bool ZwoSpec)
{
int maxLen = Math.Max(map.GetSubblkLen(0)[0], map.GetSubblkLen(1)[0]);
int minLen = Math.Min(map.GetSubblkLen(0)[0], map.GetSubblkLen(1)[0]);
int crit = ZwoSpec ? maxLen : minLen;
for (int iCell = 0; iCell < map.LocalNoOfBlocks; iCell++)
{
int iBlock = iCell + map.AggGrid.CellPartitioning.i0;
if (map.GetBlockLen(iBlock) == crit)
{
return(iCell);
}
}
return(-1);
}
public static int[] GimmeAllBlocksWithSpec(MultigridMapping map, int DOF)
{
int Bi0 = map.AggGrid.CellPartitioning.i0;
List <int> hits = new List <int>();
for (int i = 0; i < map.LocalNoOfBlocks; i++)
{
int iBlock = Bi0 + i;
int type = map.GetBlockType(iBlock);
int BlockDOF = map.GetSubblkLen(type)[0];
if (BlockDOF == DOF)
{
hits.Add(iBlock);
}
}
int[] recvcnt = hits.Count().MPIAllGather();
Debug.Assert(recvcnt.Length == map.MpiSize);
int[] VecOfAllhits = hits.ToArray().MPIAllGatherv(recvcnt);
return(VecOfAllhits);
}
