[petsc-users] question about MatCreateRedundantMatrix

[Jose E. Roman]

Michael,

In my previous email I should have checked it better. The CISS solver works indeed with dense matrices:

$ mpiexec -n 2 ./ex2 -n 30 -eps_type ciss -terse -rg_type ellipse -rg_ellipse_center 1.175 -rg_ellipse_radius 0.075 -eps_ciss_partitions 2 -mat_type dense 

2-D Laplacian Eigenproblem, N=900 (30x30 grid)

 Solution method: ciss

 Number of requested eigenvalues: 1
 Found 15 eigenvalues, all of them computed up to the required tolerance:
     1.10416, 1.10416, 1.10455, 1.10455, 1.12947, 1.12947, 1.13426, 1.13426, 
     1.16015, 1.16015, 1.19338, 1.19338, 1.21093, 1.21093, 1.24413


There might be something different in the way matrices are initialized in your code. Send me a simple example that reproduces the problem and I will track it down.

Sorry for the confusion.
Jose



> El 19 sept 2019, a las 6:20, hong--- via petsc-users <petsc-users at mcs.anl.gov> escribió:
> 
> Michael,
> We have support of MatCreateRedundantMatrix for dense matrices. For example, petsc/src/mat/examples/tests/ex9.c:
> mpiexec -n 4 ./ex9 -mat_type dense -view_mat -nsubcomms 2
> 
> Hong
> 
> On Wed, Sep 18, 2019 at 5:40 PM Povolotskyi, Mykhailo via petsc-users <petsc-users at mcs.anl.gov> wrote:
> Dear Petsc developers,
> 
> I found that MatCreateRedundantMatrix does not support dense matrices.
> 
> This causes the following problem: I cannot use CISS eigensolver from 
> SLEPC with dense matrices with parallelization over quadrature points.
> 
> Is it possible for you to add this support?
> 
> Thank you,
> 
> Michael.
> 
> 
> p.s. I apologize if you received this e-mail twice, I sent if first from 
> a different address.
>