Dear Paul,<br><br>Thank you very much for your consideration.<br>Actually, if the codes have been finished for complex number-based computation in PETSc, I can do some testings.<br>It will be ok for me. If lots of codes need to be written regarding what I said, I am not sure whether I have enough time. <br>
After all, I am not familiar with PETSc and CUDA.<br><br>Best,<br>Yujie<br><br><div class="gmail_quote">On Mon, Jan 23, 2012 at 5:09 PM, Paul Mullowney <span dir="ltr">&lt;<a href="mailto:paulm@txcorp.com">paulm@txcorp.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>

  
    
  
  <div bgcolor="#ffffff" text="#000000">
    Yujie,<br>
    <br>
    I&#39;m developing GPU code for PETSc that will, at some point, support
    complex numbers. However, I cannot give you a time frame as to when
    that will be complete as I have several other tasks that are higher
    priority right now.<span class="HOEnZb"><font color="#888888"><br>
    <br>
    -Paul</font></span><div><div class="h5"><br>
    <blockquote type="cite">Dear Barry,<br>
      <br>
      For complex number-based implementation, will the functions in
      Vec, Mat, KSP, and PC having been realized for real number work at
      least at the running level?<br>
      Do I need to do more coding for them?<br>
      <br>
      To my knowledge, Victor Minden finished most of the coding work
      for GPU implementation in PETSc with your and Matt&#39;s help.<br>
      However, it seems there is not people to further realize and
      optimize other GPU functions in PETSc.<br>
      Do you think GPU-based computation is not important or has low
      efficiency compared to current CPU-based implementation in PETSc?<br>
      <br>
      <br>
      Thank you very much,<br>
      Yujie<br>
       <br>
      <br>
      <div class="gmail_quote">On Mon, Jan 23, 2012 at 3:46 PM, Barry
        Smith <span dir="ltr">&lt;<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>&gt;</span>
        wrote:<br>
        <blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
            Yujie,<br>
          <br>
             As I said we have never tried it. You are welcome to try it
          following the directions I gave. You must realize that we do
          not have an army of people to support GPUs so you are largely
          on your on in terms of getting things done and need to use
          your initiative in trying things yourself and working through
          the problems that may arise.<br>
          <span><font color="#888888"><br>
                Barry<br>
            </font></span>
          <div>
            <div><br>
              On Jan 23, 2012, at 3:36 PM, recrusader wrote:<br>
              <br>
              &gt; Dear Paul,<br>
              &gt;<br>
              &gt; Please find the changelog in the following link for
              CUSP 0.2.0<br>
              &gt; <a href="http://code.google.com/p/cusp-library/source/browse/CHANGELOG" target="_blank">http://code.google.com/p/cusp-library/source/browse/CHANGELOG</a><br>
              &gt;<br>
              &gt; They has added cusp:complex class to support complex
              number-based operation.<br>
              &gt;<br>
              &gt; does it not work for PETSc?<br>
              &gt;<br>
              &gt; Thanks a lot,<br>
              &gt; Yujie<br>
              &gt;<br>
              &gt;<br>
              &gt; On Mon, Jan 23, 2012 at 3:28 PM, Paul Mullowney &lt;<a href="mailto:paulm@txcorp.com" target="_blank">paulm@txcorp.com</a>&gt;
              wrote:<br>
              &gt; I would claim that Petsc does not support GPU
              capabilities for complex numbers right now.<br>
              &gt;<br>
              &gt; The CUSP library is templated over the scalar type
              (i.e. float or double), however I don&#39;t think it supports
              complex numbers. Although, I could be wrong on this.<br>
              &gt;<br>
              &gt; It could be somewhat straightforward to move from
              CUSP to CUSPARSE (yes, these are different Nvidia
              libraries). I believe CUSPARSE supports SpMV for complex
              types. It also supports triangular solve for complex
              types.<br>
              &gt;<br>
              &gt; -Paul<br>
              &gt;<br>
              &gt;<br>
              &gt;     This is an installation issue, please send all
              installation issues to <a href="mailto:petsc-maint@mcs.anl.gov" target="_blank">petsc-maint@mcs.anl.gov</a>
              with the entire configure.log<br>
              &gt;<br>
              &gt;<br>
              &gt;<br>
              &gt;<br>
              &gt;    You can try editing config/PETSc/packages/cuda.py
              and removing the two lines<br>
              &gt;<br>
              &gt;     if self.scalartypes.scalartype == &#39;complex&#39;:<br>
              &gt;       raise RuntimeError(&#39;Must use real numbers with
              CUDA&#39;)<br>
              &gt;<br>
              &gt;    As it says we have never tested for complex so I
              do not know how far it is from working.<br>
              &gt;<br>
              &gt;     Barry<br>
              &gt;<br>
              &gt; On Jan 23, 2012, at 1:24 PM, recrusader wrote:<br>
              &gt;<br>
              &gt; Dear PETSc Developers,<br>
              &gt;<br>
              &gt; I am compiling PETSc-dev using GPU in complex number
              mode.<br>
              &gt; However, when I configure PETSc-dev, I get the
              following errors<br>
              &gt; &quot;<br>
              &gt;
*******************************************************************************<br>
              &gt;         UNABLE to CONFIGURE with GIVEN OPTIONS  
               (see configure.log<br>
              &gt; for details):<br>
              &gt;
-------------------------------------------------------------------------------<br>
              &gt; Cannot use cuda with complex numbers it is not coded
              for this capability<br>
              &gt;
*******************************************************************************<br>
              &gt; &quot;<br>
              &gt; I have checked CUSP. The latest version can support
              complex number.<br>
              &gt; What should I do for CUDA?<br>
              &gt; Thank you very much.<br>
              &gt;<br>
              &gt; Best,<br>
              &gt; Yujie<br>
              &gt;<br>
              &gt;<br>
              <br>
            </div>
          </div>
        </blockquote>
      </div>
      <br>
    </blockquote>
    <br>
  </div></div></div>

</blockquote></div><br>