Dear Jed,<br><br>According to the codes you added, I added them in PETSc 3.0.0p8 I am using. However, when I compile it, I get the following errors<br>&quot;libfast in: /home/yujie/codes/petsc-3.0.0-p8_complex/src/mat/impls/dense/seq/ftn-custom<br>
libfast in: /home/yujie/codes/petsc-3.0.0-p8_complex/src/mat/impls/dense/mpi<br>mpidense.c:1583: error: invalid conversion from &#39;PetscErrorCode (*)(_p_Mat*)&#39; to &#39;PetscErrorCode (*)(_p_Mat*, _p_Mat*, _p_Mat*)&#39;<br>
mpidense.c:1583: error: invalid conversion from &#39;PetscErrorCode (*)(_p_Mat*)&#39; to &#39;PetscErrorCode (*)(_p_Mat*, PetscInt, const PetscScalar*)&#39;<br>/usr/bin/ar: mpidense.o: No such file or directory<br>&quot;<br>
I use &quot;C++&quot; mode for compilation. In PETsc I am using , the codes are around Line 1583:<br><br>&quot;<br>/*100*/0,<br>       0,<br>       0,<br>       MatConjugate_MPIDense, //by Yujie<br>       0,<br>/*105*/0,<br>
       MatRealPart_MPIDense,<br>       MatImaginaryPart_MPIDense,<br>       0,<br>       0<br>};  //here is 1583<br><br>EXTERN_C_BEGIN<br>#undef __FUNCT__<br>#define __FUNCT__ &quot;MatMPIDenseSetPreallocation_MPIDense&quot;<br>
&quot;<br>Thanks a lot.<br><br>Regards,<br>Yujie <br><br><div class="gmail_quote">On Thu, Feb 25, 2010 at 10:53 AM, Jed Brown <span dir="ltr">&lt;<a href="mailto:jed@59a2.org">jed@59a2.org</a>&gt;</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;">
<div class="im">On Thu, 25 Feb 2010 12:43:35 -0600, Yujie &lt;<a href="mailto:recrusader@gmail.com">recrusader@gmail.com</a>&gt; wrote:<br>
&gt; Dear Jed,<br>
&gt;<br>
&gt; Should You add the judgemnet to the complex number?<br>
<br>
</div>It should be fine.  PetscConj, PetscRealPart, PetscImaginaryPart are<br>
defined without PETSC_USE_COMPLEX, and the loop is completely compiled<br>
out at high optimization levels.<br>
<font color="#888888"><br>
Jed<br>
</font></blockquote></div><br>