$ time ./ex7.exe -f1 k.dat -f2 m.dat -eps_gen_hermitian -eps_type lapack -eps_smallest_real > x.out 2>&1<br><br>real 19m4.487s<br>user 18m19.650s<br><div>sys 0m1.762s<br></div><div><br></div><div>---John</div>
<div><br></div><br><div class="gmail_quote">On Wed, Jul 27, 2011 at 7:03 AM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Try running with -eps_gen_hermitian (since ex7 does not assume that the problem is symmetric).<br>
Jose<br>
<br>
<br>
<br>
El 27/07/2011, a las 12:51, John Chludzinski escribió:<br>
<div><div class="h5"><br>
> I let the SLEPc code run for ~45 min. when it terminated with the same values I was getting using DSYGV in LAPACK.<br>
><br>
> If I write code to directly call LAPACK (i.e., DSYGV), it uses ~3.93 min. What's up with this?<br>
><br>
> ---John<br>
><br>
><br>
> On Wed, Jul 27, 2011 at 5:29 AM, John Chludzinski <<a href="mailto:jchludzinski@gmail.com">jchludzinski@gmail.com</a>> wrote:<br>
> I'm trying to create a dense matrices from values I'm reading from (binary) files. I tried the following code:<br>
><br>
> Mat A;<br>
> int n = SIZE; //4002<br>
> double *K = (double *)calloc( sizeof(double), SIZE*SIZE );<br>
> ...<br>
> MatCreateSeqDense(PETSC_COMM_SELF, n, n, K, &A);<br>
> MatView(A,PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD));<br>
> ierr = PetscFinalize();CHKERRQ(ierr);<br>
><br>
> NOTE:*** I'm converting K to the FORTRAN column major from the C row major order before I call MatCreateSeqDense(...).<br>
><br>
> 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:<br>
><br>
> ./ex7.exe -f1 k.dat -f2 m.dat -eps_type lapack -eps_smallest_real<br>
><br>
> Am I creating the proper PETSc binary (canonical) format for my 2 matrices?<br>
><br>
> ---John<br>
><br>
><br>
><br>
<br>
</div></div></blockquote></div><br>