<div dir="ltr">On Tue, Jul 2, 2013 at 12:19 PM, Michael Povolotskyi <span dir="ltr"><<a href="mailto:mpovolot@purdue.edu" target="_blank">mpovolot@purdue.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    <div>On 07/02/2013 12:46 PM, Hong Zhang
      wrote:<br>
    </div>
    <blockquote type="cite">
      <div dir="ltr">
        <div class="gmail_extra">Michael :<br>
          <div class="gmail_quote">
            <blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Dear
              Petsc developers,<br>
              I'm gradually moving from the version 3.2 to the version
              3.4.<br>
              I had to skip the version 3.3 because of the bug in
              MatMatMult that has been reported to you.<br>
              <br>
              What I see is that the MatMatMult works well but there is
              a difference:<br>
              in the version 3.2 if a sparse matrix was multiplied by a
              dense matrix the result was dense.<br>
              In the version 3.4 the result is sparse.<br>
            </blockquote>
            <div><br>
            </div>
            <div> A sparse matrix multiplied by a dense matrix results
              in a dense matrix for all petsc versions.</div>
            <div>What makes you conclude "In the version 3.4
              the result is sparse"?</div>
            <div><br>
            </div>
            <div>Hong</div>
            <blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span><font color="#888888"><br>
                  <br>
                </font></span></blockquote>
          </div>
          <br>
        </div>
      </div>
    </blockquote>
    I use MatMatMult with MAT_INITIAL_MATRIX<br>
    Then I call MatGetType(A, &type) for the product matrix<br>
    Then I have an if statement:<br>
    if (string(type) == "seqdense" || string(type) ==  "mpidense")</div></blockquote><div><br></div><div style>What is the type? M->hdr.type_name</div><div style><br></div><div style>   Matt</div><div>
 </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><span class="HOEnZb"><font color="#888888"><br>
    Michael.<br>
  </font></span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>