> On Sep 20, 2019, at 12:35 AM, Young Hyun Jo wrote:
>
> Thanks for the answer.
> It's really helpful to understand the PETSc library.
> By the way, I just want to ask two more questions not related to the
> multigrid.
>
> 1. Is there a method known as the fastest solver for Poisson's
The DMDA structured grid management in PETSc does not provide the needed
interpolations for doing non-nested multigrid. The algebraic portions of
PETSc's geometric multigrid would work fine for that case if you have your own
way to provide the needed interpolation. I'll note that
You didn't indicate "why" you can't use multiple levels with "130 grids", is
there some error message? Nor do you mention if you have periodic boundary
conditions or are using cell or vertex centered unknowns. All of these things
affect when you can coarsen for multigrain or not.
I think you are fine with DMForest. I just mentioned ForestClaw for
background. It has a bunch of hyperbolic stuff in there that is specialized.
On Thu, Sep 19, 2019 at 8:16 AM Mohammad Hassan
wrote:
> In fact, I would create my base mesh in DMPlex and use DMForest to
> construct the
Hong,
do you have in mind a reason why the newer version should work or is it a
general recommendation?
Which stable version would you recommend to upgrade to?
Thank you,
Michael.
On 09/19/2019 02:22 PM, Zhang, Hong wrote:
Michael,
Michael,
--
[0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: Mat type seqdense
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
for trouble shooting.
[0]PETSC ERROR: Petsc Release
Hello Jose,
I have done the test case to reproduce my error.
1. You will need to download a file "matrix.bin" from the following link
https://www.dropbox.com/s/6y7ro99ou4qr8uy/matrix.bin?dl=0
2. Here is the C++ code I use
#include
#include
#include
#include
#include
#include
#include
In fact, I would create my base mesh in DMPlex and use DMForest to construct
the non-conformal meshing obtained by my own block-based AMR functions.
ForestClaw is based on p4est. However, I may need to implement the AMR
algorithm on DM in p4est library and then convert it to DMPlex. Do you
Note, Forest gives you individual elements at the leaves. Donna Calhoun, a
former Chombo user, has developed a block structured solver on p4est (
https://math.boisestate.edu/~calhoun/ForestClaw/index.html), but I would
imagine that you could just take the Plex that DMForest creates and just
call
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
10 matches
Mail list logo
