Victor A. P. Magri
Victor A. P. Magri
Hi @Heatherms27, there aren't plans now for adding support of complex numbers to preconditioners in the IJ interface. Sorry about that! Currently, many compilation warnings are issued when `--enable-complex` is...
Hi Jed, I'm not sure about your first question. Regarding the second, have you tried the unknown-based approach of BoomerAMG? In this approach, interpolation is computed only within the same...
Thanks for the instructions, @jedbrown! I'm able to run your test case, but the number of iterations that I'm getting is different than yours. Could you double-check if my command-line...
@jedbrown, I see that PETSc changes its default options for BoomerAMG when using the device, see lines [lines 1950-1962](https://petsc.org/main/src/ksp/pc/impls/hypre/hypre.c.html#:~:text=1950%3A%20%23if%20defined(PETSC_HAVE_HYPRE_DEVICE)). Alternatively, you could match the options used for GPUs on the...
Jed, this is the closest GPU-friendly configuration, in terms of convergence, to the best CPU-friendly configuration of BoomerAMG for your problem: ``` -mg_coarse_pc_hypre_boomeramg_coarsen_type pmis -mg_coarse_pc_hypre_boomeramg_interp_type ext+i -mg_coarse_pc_hypre_boomeramg_no_CF -mg_coarse_pc_hypre_boomeramg_P_max 6 -mg_coarse_pc_hypre_boomeramg_print_statistics...
That makes sense, @liruipeng! Thanks! Perhaps, I was not very careful in making the suggestion to change the defaults... Please, disregard that @jedbrown. Still, as Ruipeng said, I believe it...
Sounds good, @jedbrown. Let us know the results you get. Thanks!
Hi @yufeimi, It seems that the issue comes from PETSc calling `MPI_Comm_dup` multiple times during execution of your program, [see this line]( https://gitlab.com/petsc/petsc/-/blob/main/src/mat/impls/hypre/mhypre.c#L2130). Supposing that a hypre matrix is always...
Hi @rfalgout, I found out that in one of my SSAMG tests, the value of `gamma` is negative and thus `gamma > HYPRE_REAL_MIN` is false but not because gamma is...
Hello Ruipeng and Rob, I agree `alpha` and `gamma` will be positive for an SPD `inv(M)` operator and non-negative for the SPSD case. Since SSAMG shows a negative `gamma` for...