Skip to content

[bugfix] changed value-type to double in vector and matrix in JacobiNQuadratureRule1D

Simon Praetorius requested to merge simon.praetorius/dune-geometry:issue/quadrature_rules_ct into master

Summary

changed value type to double in vector and matrix in eigenvalue computations in JacobiNQuadratureRule1D

Details

When instantiating JacobiNQuadratureRule1D with ct=float, the function DynamicMatrixHelp::eigenValuesNonSym is not found since float and double are mixed in the argument types. The reason is that the matrix was filled with ct and the vector was filled with complex<double>.

The bug is introduced in !127 (merged)

Edited by Simon Praetorius

Merge request reports