Add in toggle to matrixequations

[jlperla]

• [x] We should get support for the sylvester solver choice. I would change https://github.com/HighDimensionalEconLab/DifferentiableStateSpaceModels.jl/blob/main/src/types.jl#L331-L346

To add in somehing like

sylvester_solver = :MatrixEquations

As a default, where we can have it as :MatrixEquations vs. :GeneralizedSylvesterSolver and pick whichever one works the best for now.

• [x] Then I think in the code we just need to change code like https://github.com/HighDimensionalEconLab/DifferentiableStateSpaceModels.jl/blob/main/src/generate_perturbation.jl#L308 into somethingg like

  if settings.sylvester_solver == :GeneralizedSylvesterSolver
  ws = GeneralizedSylvesterWs(n, n, n_x, 2)
  generalized_sylvester_solver!(buff.A, buff.C, c.h_x, buff.E, 2, ws)
  else
  # whatever the old code was
  end

• [x] In the case of the derivatives, I think it is totally fine for us to keep the workspace around - even if we actually execute with MatrixEquations (i.e. https://github.com/HighDimensionalEconLab/DifferentiableStateSpaceModels.jl/blob/main/src/generate_perturbation_derivatives.jl#L111 can stay there even if https://github.com/HighDimensionalEconLab/DifferentiableStateSpaceModels.jl/blob/main/src/generate_perturbation_derivatives.jl#L133 ends up using MatrixEquations If we find later that we mostly want to use the MatrixEquations then we can get fancier.