
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=windows-1256">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi,<br>

I have a valgrind "invalid read" in MatGetBrowsOfAoCols_MPIAIJ in the lines<br>

 rvalues  = gen_from->values; /* holds the length of receiving row */<br>

 svalues  = gen_to->values;   /* holds the length of sending row */<br>

(and others, but that are the first). This function is called indirectly via MatMatMultSymbolic.<br>

All invalid reads apear in MatGetBrowsOfAoCols_MPIAIJ.<br>

The calling sequence for the matrix factors are<br>

  MatCreate(PETSC_COMM_WORLD,A);<br>

  MatSetSizes(A,NumLocRows,NumLocCols,PETSC_DETERMINE,PETSC_DETERMINE);<br>

  MatSetType(A,MATMPIAIJ);<br>

  MatMPIAIJSetPreallocation(A,PETSC_NULL,nnz,PETSC_NULL,onz);<br>

  MatSetValues(...);<br>

  MatAssemblyBegin(...)<br>

  MatAssemblyEnd(...)<br>

and then a call to MatMatMultSymbolic. <br>

I want the result matrix to be created. In the documentation I see, that this happens, if I pass MAT_REUSE_MATRIX as scall, but how can I pass the value, I don't the parameter.<br>

Where can I find an example for MatMatMultSymbolic / MatMatMultNumeric?<br>

Thanks,<br>

    Hannes<br>

</div>

</body>

</html>