Dear Barry,<br><br>What are the difficult things if PETSc is revised with matrices of complex and real numbers? It should be more flexible for a general scientific toolbox. I am curious almost all the packages don't support both simultaneously. Thanks a lot.<br>
<br>Regards,<br>Yujie <br><br><div class="gmail_quote">On Fri, Jan 22, 2010 at 2:35 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
The way we handle complex numbered linear systems in PETSc is to compile all of PETSc with complex numbers and then just use the solvers on those complex numbers. The current drawback to this is that PETSc can only be built with support for complex numbers or for real numbers. We cannot build a PETSc where some matrices are complex and some are real.<br>
We don't have any interest in solving complex systems as larger real systems.<br><font color="#888888">
<br>
Barry</font><div><div></div><div class="h5"><br>
<br>
On Jan 22, 2010, at 11:04 AM, Yujie wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear PETSc Developers,<br>
<br>
Recently, I am trying to find some complex number-based solvers and preconditioners. However, it is difficult to find a general framework to include some solvers and preconditioners. Trilinos is developing a package, komplex, to use the real-number-based solver to solve complex number -based problem. I don't know whether PETSc wants to develop such the function for complex number-based problem. I think it will significantly increase the application range of PETSc. After all, in PETSc, lots of solvers and preconditioners have been developed. Thanks a lot.<br>
<br>
Regards,<br>
Yujie<br>
</blockquote>
<br>
</div></div></blockquote></div><br>