GridapSolvers.MultilevelTools

generate_level_parts(root_parts::AbstractArray,num_procs_x_level::Vector{<:Integer})

From a root communicator root_parts, generate a sequence of nested subcommunicators with sizes given by num_procs_x_level.

HierarchicalArray{T,A,B} <: AbstractVector{T}

Array of hierarchical (nested) distributed objects. Each level might live in a different subcommunicator. If a processor does not belong to subcommunicator ranks[i], then array[i] is nothing.

However, it assumes: - The subcommunicators are nested, so that ranks[i] contains ranks[i+1]. - The first subcommunicator does not have empty parts.

Base.map(f::Function,args::Vararg{HierarchicalArray,N}) where N

Maps a function to a set of HierarchicalArrays. The function is applied only in the subcommunicators where the processor belongs to.

with_level(f::Function,a::HierarchicalArray,lev::Integer;default=nothing)

Applies a function to the lev-th level of a HierarchicalArray. If the processor does not belong to the subcommunicator of the lev-th level, then default is returned.

const ModelHierarchy = HierarchicalArray{<:ModelHierarchyLevel}

A ModelHierarchy is a hierarchical array of ModelHierarchyLevel objects. It stores the adapted/redistributed models and the corresponding subcommunicators.

For convenience, implements some of the API of DiscreteModel.

Single level for a ModelHierarchy.

Note that model_red and red_glue might be of type Nothing whenever there is no redistribution in a given level.

ref_glue is of type Nothing on the coarsest level.

CartesianModelHierarchy(
  ranks::AbstractVector{<:Integer},
  np_per_level::Vector{<:NTuple{D,<:Integer}},
  domain::Tuple,
  nc::NTuple{D,<:Integer};
  nrefs::Union{<:Integer,Vector{<:Integer},Vector{<:NTuple{D,<:Integer}},NTuple{D,<:Integer}},
  map::Function = identity,
  isperiodic::NTuple{D,Bool} = Tuple(fill(false,D)),
  add_labels!::Function = (labels -> nothing),
) where D

Returns a ModelHierarchy with a Cartesian model as coarsest level. The i-th level will be distributed among np_per_level[i] processors. Two consecutive levels are refined by a factor of nrefs[i].

Parameters:

• ranks: Initial communicator. Will be used to generate subcommunicators.
• domain: Tuple containing the domain limits.
• nc: Tuple containing the number of cells in each direction for the coarsest model.
• np_per_level: Vector containing the number of processors we want to distribute each level into. Requires a tuple np = (np_1,...,np_d) for each level, then each level will be distributed among prod(np) processors with np_i processors in the i-th direction.
• nrefs: Vector containing the refinement factor for each level. Has nlevs-1 entries, and each entry can either be an integer (homogeneous refinement) or a tuple with D integers (inhomogeneous refinement).

P4estCartesianModelHierarchy(
  ranks,np_per_level,domain,nc::NTuple{D,<:Integer};
  num_refs_coarse::Integer = 0,
  add_labels!::Function = (labels -> nothing),
  map::Function = identity,
  isperiodic::NTuple{D,Bool} = Tuple(fill(false,D))
) where D

Returns a ModelHierarchy with a Cartesian model as coarsest level, using GridapP4est.jl. The i-th level will be distributed among np_per_level[i] processors. The seed model is given by cmodel = CartesianDiscreteModel(domain,nc).

const FESpaceHierarchy = HierarchicalArray{<:FESpaceHierarchyLevel}

A FESpaceHierarchy is a hierarchical array of FESpaceHierarchyLevel objects. It stores the adapted/redistributed fe spaces and the corresponding subcommunicators.

For convenience, implements some of the API of FESpace.

abstract type LocalProjectionMap{T} <: Map end

struct ReffeProjectionMap{T} <: LocalProjectionMap{T}
  op      :: Operation{T}
  reffe   :: Tuple{<:ReferenceFEName,Any,Any}
  qdegree :: Int
end

Map that projects a field/field-basis onto another local reference space given by a ReferenceFE.

Example:

model = CartesianDiscreteModel((0,1,0,1),(2,2))

reffe_h1 = ReferenceFE(QUAD,lagrangian,Float64,1,space=:Q)
reffe_l2 = ReferenceFE(QUAD,lagrangian,Float64,1,space=:P)
U = FESpace(model,reffe_h1)
u_h1 = interpolate(f,U)

Ω = Triangulation(model)
dΩ = Measure(Ω,2)

Π = LocalProjectionMap(reffe_l2)
u_l2 = Π(u_h1,dΩ)

struct SpaceProjectionMap{T} <: LocalProjectionMap{T}
  op      :: Operation{T}
  space   :: A
  qdegree :: Int
end

Map that projects a CellField onto another FESpace. Necessary when the arrival space has constraints (e.g. boundary conditions) that need to be taken into account.