[petsc-users] MatCreateSeqDense ?

[John Chludzinski]

$ 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
> >
> >
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110727/bd4e3646/attachment.htm>