I have a 4-dimensional FESpace with approximately 500k-1M degrees of freedom. Given the (partial) translational symmetry my hope is to split the problems in order to speed up the calculation.
The given total matrix to solve is complex and neither symmetric nor positive definite.
Is there a recommendation which solver/strategy I should use.