Hi Maxi,
The AffineConstraints class does not support using mixed precision and I am
not aware of any alternatives that can solve this problem directly.
There may be a better solution, but one way to solve this problem would be
to set up two different AffineConstraints objects (one for each float
I tried to implement a simple mixed-precision algorithm in deal.II, which
creates my matrix using float values, and my vectors using double values.
Unfortunately, the function distribute_local_to_global(cell_matrix,
cell_rhs, local_dof_indices, system_matrix, system_rhs) does not support
mixed
