Allow SuperLU_DIST to solve with repeated right-hand sides (#4074)

* Make copy of MatrixCSR values

* clang-format

* Failing test - to fix.

* Header in place - need to write implementation.

* Fix

* Done scalepermstruct and lustruct for double type.

* Looks clean.

* Works with two solves.

* Re-use LU factorisation works - bit more boiler plate and then finished

* Pass A by ref

* C++ builds - should probably split out dull stuff to own file.

* Builds

* Tests pass.

* Refactor

* Tidy

* SolveStructDeleter - cleanup should happen at end.

* Fix bug.

* Works - TODO: Add throw on doing second solve without correct factor mode set.

* Tidy.

* Removing unused structs

Removed unused struct definitions for dSOLVE, sSOLVE, and zSOLVE.

* Simplify logging messages

* Tidy and docs

* Bit more checking

* Tweaks

* Only call if solve initialised

* Add setter for A (advanced use).

* Wrap SuperLUDistMatrix - useful for advanced use reusing factors.

* Wrappers.

* More wrapping.

* Wrap set_A. Quite cheap - setter?

* More testing - need splitting into separate tests.

* Add dtype, improve documentation.

* Some more fine tuning

* Tweak

* Better warning message.

* Finish tests.

* Fix docs

* Headers not needed anymore

* Simplify last part of test

Author: jhale