Don't use "time" to measure performance; instead use -log_summary or PetscGetTime for the interesting part of the computation. In this case, computing the residuals will take a lot of time. Jose
El 27/07/2011, a las 19:28, John Chludzinski escribi?: > $ time ./ex7.exe -f1 k.dat -f2 m.dat -eps_gen_hermitian -eps_type lapack > -eps_smallest_real > x.out 2>&1 > > real 19m4.487s > user 18m19.650s > sys 0m1.762s > > ---John > > > On Wed, Jul 27, 2011 at 7:03 AM, Jose E. Roman <jroman at dsic.upv.es> wrote: > Try running with -eps_gen_hermitian (since ex7 does not assume that the > problem is symmetric). > Jose > > > > El 27/07/2011, a las 12:51, John Chludzinski escribi?: > > > I let the SLEPc code run for ~45 min. when it terminated with the same > > values I was getting using DSYGV in LAPACK. > > > > If I write code to directly call LAPACK (i.e., DSYGV), it uses ~3.93 min. > > What's up with this? > > > > ---John > > > > > > On Wed, Jul 27, 2011 at 5:29 AM, John Chludzinski <jchludzinski at > > gmail.com> wrote: > > I'm trying to create a dense matrices from values I'm reading from (binary) > > files. I tried the following code: > > > > Mat A; > > int n = SIZE; //4002 > > double *K = (double *)calloc( sizeof(double), SIZE*SIZE ); > > ... > > MatCreateSeqDense(PETSC_COMM_SELF, n, n, K, &A); > > MatView(A,PETSC_VIEWER_BINARY_(PETSC_COMM_WORLD)); > > ierr = PetscFinalize();CHKERRQ(ierr); > > > > NOTE:*** I'm converting K to the FORTRAN column major from the C row major > > order before I call MatCreateSeqDense(...). > > > > 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: > > > > ./ex7.exe -f1 k.dat -f2 m.dat -eps_type lapack -eps_smallest_real > > > > Am I creating the proper PETSc binary (canonical) format for my 2 matrices? > > > > ---John > > > > > > > >
