Add Skew map to BCO

[knelli2]

Proposed changes

Adds the Skew map to the BinaryCompactObject domain. Sets it to None in all input files for now until a control system is written for it.

Part of #4603.

Upgrade instructions

When using the BinaryCompactObject domain creator, you must add the following option to the time depndent options:

  TimeDependentOptions:
    SkewMap:
      InitialValuesY: [0.0, 0.01, 0.0]
      InitialValuesZ: [0.0, 0.01, 0.0]

It has the same capabilities to be None or from a volume file as the other time dependent map options.

Code review checklist

• The code is documented and the documentation renders correctly. Run
  make doc to generate the documentation locally into BUILD_DIR/docs/html.
  Then open index.html.
• The code follows the stylistic and code quality guidelines listed in the
  code review guide.
• The PR lists upgrade instructions and is labeled bugfix or
  new feature if appropriate.

Further comments

Depends on and includes #6446.

[knelli2]

Rebased and no longer dependent

[nilsvu]

If use_rigid_map=false but both rot_scale_trans and skew maps are given, this will ignore the skew map. Is that what you want? Could you make this clearer, maybe with a comment or by throwing errors or something

[knelli2]

Yes this is correct. I'll add a comment

[nilsvu]

So when use_rigid_map=false the skew map can't be applied? Better throw an error then if the skew map was requested.

[knelli2]

Not "can't" be applied, but "won't" be applied. If there is a skew map, but the block we are in is in the outer shell, then use_rigid_map=false and skew_map_.has_value()=true, yet we don't want to apply the skew map

[nilsvu]

Ok that clarifies it, thanks 👍

I think in this code branch and also in the block_has_shape_map branch in grid_to_distorted_map below we're always within the inner cubes, so use_rigid_map should be always true. Maybe assert that instead of the sentinel.

[knelli2]

Added some errors