<div>I let the SLEPc code run for ~45 min. when it terminated with the same values I was getting using DSYGV in LAPACK.<br></div><div><br></div><div>If I write code to directly call LAPACK (i.e., DSYGV), it uses ~3.93 min. What's up with this?</div>
<div><br></div><div>---John</div><div><br></div><br><div class="gmail_quote">On Wed, Jul 27, 2011 at 5:29 AM, John Chludzinski <span dir="ltr"><<a href="mailto:jchludzinski@gmail.com">jchludzinski@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div>
<font face="courier new,monospace">I'm trying to create a dense matrices from values I'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, &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'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><font color="#888888"><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"><br></font></div><div><font face="courier new,monospace"><br></font></div>
</font></blockquote></div><br>