<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div dir="ltr">
<div dir="ltr">Michael,<br>
</div>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
--------------------------------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object type<br>
[0]PETSC ERROR: Mat type seqdense<br>
[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">
http://www.mcs.anl.gov/petsc/documentation/faq.html</a> <br>
for trouble shooting.<br>
[0]PETSC ERROR: Petsc Release Version 3.8.4, Mar, 24, 2018<br>
</blockquote>
<div><br>
</div>
<div>This is an old version of Petsc. Can you update to the latest Petsc release?</div>
<div>Hong</div>
<div><br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
On 09/19/2019 04:55 AM, Jose E. Roman wrote:<br>
> Michael,<br>
><br>
> In my previous email I should have checked it better. The CISS solver works indeed with dense matrices:<br>
><br>
> $ 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<br>
><br>
> 2-D Laplacian Eigenproblem, N=900 (30x30 grid)<br>
><br>
> Solution method: ciss<br>
><br>
> Number of requested eigenvalues: 1<br>
> Found 15 eigenvalues, all of them computed up to the required tolerance:<br>
> 1.10416, 1.10416, 1.10455, 1.10455, 1.12947, 1.12947, 1.13426, 1.13426,<br>
> 1.16015, 1.16015, 1.19338, 1.19338, 1.21093, 1.21093, 1.24413<br>
><br>
><br>
> 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.<br>
><br>
> Sorry for the confusion.<br>
> Jose<br>
><br>
><br>
><br>
>> El 19 sept 2019, a las 6:20, hong--- via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> escribió:<br>
>><br>
>> Michael,<br>
>> We have support of MatCreateRedundantMatrix for dense matrices. For example, petsc/src/mat/examples/tests/ex9.c:<br>
>> mpiexec -n 4 ./ex9 -mat_type dense -view_mat -nsubcomms 2<br>
>><br>
>> Hong<br>
>><br>
>> On Wed, Sep 18, 2019 at 5:40 PM Povolotskyi, Mykhailo via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
>> Dear Petsc developers,<br>
>><br>
>> I found that MatCreateRedundantMatrix does not support dense matrices.<br>
>><br>
>> This causes the following problem: I cannot use CISS eigensolver from<br>
>> SLEPC with dense matrices with parallelization over quadrature points.<br>
>><br>
>> Is it possible for you to add this support?<br>
>><br>
>> Thank you,<br>
>><br>
>> Michael.<br>
>><br>
>><br>
>> p.s. I apologize if you received this e-mail twice, I sent if first from<br>
>> a different address.<br>
>><br>
<br>
</blockquote>
</div>
</div>
</body>
</html>