Santiago,
The shift '-eps_target -2e-3+1.01i' is very close to the eigenvalues. What
happens if you pick a target little away from your eigenvalues?
I suspect mumps encounters a zero pivot during numerical factorization. There
are options to handle it, but I need matrices A and B to investigate.
Yeah this doesn't look good for MUMPS but isn't for sure the problem either.
The valgrind output should be sent to the MUMPS developers.
Hong,
Can you send this to the MUMPS developers and see what they say?
Thanks
Barry
> On Oct 30, 2018, at 2:04 PM, Santiago
Dear users,
Following a suggestion from Matthew Knepley I’ve been trying to apply
fieldsplit/gamg for my set of PDEs but I’m still encountering issues despite
various tests. pc_gamg simply won’t start.
Note that direct solvers always yield the correct, physical result.
Removing the fieldsplit
And you can also always use -mg_levels_esteig_ksp_type cg
On Tue, Oct 30, 2018 at 2:51 PM Mark Adams wrote:
> The last two cases were cut off but that is OK.
>
> Lets just do -pc_gamg_square_graph 0 and please add/do -ksp_view -
> mg_levels_ksp_type richardson -pc_mg_levels 5 -gamg_est_ksp_type
The last two cases were cut off but that is OK.
Lets just do -pc_gamg_square_graph 0 and please add/do -ksp_view -
mg_levels_ksp_type richardson -pc_mg_levels 5 -gamg_est_ksp_type cg
Send the grep GAMG again, and the ksp_view output.
And then do this again with: -pc_gamg_agg_nsmooths
I think
On Tue, Oct 30, 2018 at 1:18 PM Santiago Andres Triana via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi petsc-users,
>
> I am solving a generalized eigenvalue problem using ex7 in
> $SLEPC_DIR/src/eps/examples/tutorials/. I provide the A and B matrices.
> The program runs fine, with correct
Can you run the code on the "failing" machine using valgrind?
https://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
Barry
> On Oct 30, 2018, at 12:10 PM, Santiago Andres Triana via petsc-users
> wrote:
>
> Hi petsc-users,
>
> I am solving a generalized eigenvalue problem
Nicolas,
Smoothed aggregation is fine with shells. see the original SA paper (
https://link.springer.com/article/10.1007/BF02238511).
The rotational modes, which are the non-trivial modes that must be
supplied, are used in the interpolation.
Mark
On Tue, Oct 30, 2018 at 5:22 AM Karin wrote:
Thank you, Satish.
I apologize for my undetailed question. The makefile problem usually happens
due to different environment variables, so I find the administrator of the
cluster and solve the problem.
Another question is about function DMDACreate3d, const PetscInt lx[] needs to
be defined if I
