[Libmesh-users] Using complex DenseMatrix with Real compiled libmesh library

[Mathias N. <mne...@tu...>]

Hello!

Since my system matrices are "Real" valued I'm using a real compiled 
libmesh library. But on element level I need complex matrices. For this 
reason I made explicit instantiations of the class DenseMatrix<Complex> 
( This is only done when using the preprocessor directive "#ifdef 
USE_COMPLEX_NUMBERS" ).

To get my code compiled, due to the complex instantiation of class 
DenseMatrix, I have to uncomment in file "dense_matric.C", in method 
"DenseMatrix<T>::_cholesky_decompose ()" the part

#ifndef USE_COMPLEX_NUMBERS
  if (A(i,j) <= 0.0)
    {
    std::cerr << "Error! Can only use Cholesky decomposition "
              << "with symmetric positive definite matrices."
              << std::endl;
    libmesh_error();
    }
#endif

This makes sense since "A(i,j) <= 0.0" isn't defined using a complex and 
a double value.

Maybe there can be used a alternative code line in future so that it is 
possible to use complex dense matrices.

Thanks in advance,
Mathias