[petsc-users] MatCreateSeqDense ?

[Jose E. Roman]

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
> 
> 
>