<div> 

<font face="courier new,monospace">I&#39;m trying to create a dense matrices from values I&#39;m reading from (binary) files.  I tried the following code:</font></div><div><font face="courier new,monospace"><br></font></div>
<div><font face="courier new,monospace">Mat A;</font><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">int n = SIZE; //4002</font></div><div><font face="courier new,monospace">double *K = (double *)calloc( sizeof(double), SIZE*SIZE );</font><font face="courier new,monospace">
</font></div><div><font face="courier new,monospace">...</font><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">M</font><font face="courier new,monospace">atCreateSeqDense(PETSC_COMM_SELF, n, n, K, &amp;A);</font><font face="courier new,monospace"> </font></div>
<div><font face="courier new,monospace">MatView(A,PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD));</font><font face="courier new,monospace"><br></font><font face="courier new,monospace">  ierr = PetscFinalize();CHKERRQ(ierr);</font><font face="courier new,monospace"><br>
</font></div><div><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">NOTE:*** I&#39;m converting K to the FORTRAN column major from the C row major order before I call </font>
<font face="courier new,monospace">M</font><font face="courier new,monospace">atCreateSeqDense(...).</font>
<font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">This appears to work but when I try to use the 2 matrices I thus created with SLEPc ex7 (generalized eigenvalue problem) it never terminates, using:</font><font face="courier new,monospace"><br>
</font></div><div><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">./ex7.exe -f1 k.dat -f2 m.dat -eps_type lapack -eps_smallest_real</font><font face="courier new,monospace"><br>
</font></div><div><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace">Am I creating the proper PETSc binary (canonical) format for my 2 matrices?<br></font></div><div><font face="courier new,monospace"><br>
</font></div><div><font face="courier new,monospace">---John<br></font></div><div><font face="courier new,monospace"></font><font face="courier new,monospace"><br></font></div><div><font face="courier new,monospace"><br></font></div>