Skip to content

[bugfix] complex values in FCGSolvers

Simon Praetorius requested to merge simon.praetorius/dune-istl:issue/iterative_solver_complex into master

Summary

Fix complex<double> field type in RestartedFCGSolver and CompleteFCGSolver

Details

Passing vector type with field_type=complex<double> results in compiler errors due to assignment of complex to double. This is fixed by changing types to field_type instead of real_type similar to CGSolver and PCGSolver.

A test with explicit template instantiation of complex type is added to catch this error.

Edited by Simon Praetorius

Merge request reports