Some parts of the code expect an empty set of active elements
if all of the elements are active.

The previous implementation was looking for one deactivated
subdomain, but there could be multiple subdomains inactive
at the same time. Now all deactivated subdomains are
searched for active/inactive elements.

Allowing not only a zero value to be applied to the inactive
nodes but also heterogeneous, time-dependent values.

Using elements centers for the condition.

Creating the curve for time interval s.t. the whole material ids group
is deactivated, as it was for the time interval case.

The set of deactivated nodes is now constructed on the fly using the
active elements set passed from the process variable.

Code simplification and one less dereferencing at run time.

Same thing is checked after the find_if.

The name reflects the special nature of the nodes,
namely that these are nodes not touching the active
elements.

Use const references where necessary.

In PV::createBoundaryConditionsForDeactivatedSubDomains method the boundary
conditions must be created only for the number of components of the specified
variable instead for all components of all variables.

The fmt library happily accepts std::string without
any conversions.

Include spdlog via BaseLib/Logging.h, which
provides macro replacements for
ERR/DBUG/WARN/INFO...

Fixes multiple copyright messages in the doxygen docu.

ag '\\file' -L --ignore=ThirdParty *Lib |xargs  sed -i 's/\(.*\)\\copyright/\1\\file\n\1\\copyright/'

See previous commit for details.

std::vector<std::size_t> _ids_of_active_elements.

due to its std::unique_ptr type instance.

Use the specified patch mesh instead of the entire bulk mesh to define
the source term.