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
------------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It is the best place to buy or sell services for
just about anything Open Source.
http://p.sf.net/sfu/Xq1LFB
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
