<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <div class="moz-cite-prefix">On 07/02/2013 12:46 PM, Hong Zhang
      wrote:<br>
    </div>
    <blockquote
cite="mid:CAGCphBtJyjni7DDAX3Q+pMUCH56U82tjaccgCvCTOkYaK5TREg@mail.gmail.com"
      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 style="">What makes you conclude "In the version 3.4
              the result is sparse"?</div>
            <div style=""><br>
            </div>
            <div style="">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
                class=""><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")<br>
    <br>
    Michael.<br>
  </body>
</html>