Some refactorings.

move F out of the loop

node_xs as 2-array of global dim vectors. rename

eliminate node var unused

Extract i from Df

Extract Df out of comp loop

Extract Df and f out of node loop

Extract l and l_original_norm out of f, Df

local rhs stride

Block operations. Local x, node_local_indices

the latter is the key for all node 0/1 specific operations

Remove unused version of findElementOrError with iterators
and error string.

Removing unique_ptr and constness from the members allowing removing the
constructors.

Defaulting ProcessVariable move ctor.

The previous method was based on material ids, but now
it is gone there must be another method to find out
whether an element is active or not.

A deactivated subdomain consists of a single set of elements
(with inner and outer nodes) instead of multiple meshes, one
for each material id.

In the former implementation the artificially created
boundary between two deactivated subdomains was not handled
correctly resulting in singular global matrices.

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.