Re: [deal.II] On using block diagonal preconditioner for the mixed Poisson problem

On 8/27/19 2:23 PM, Charlie Z wrote: > > I find a conflict here is that:: > 1) TrilinosWrappers::PreconditionAMG can only be initialized with a > matrix, but not a ApproximateSchurComplement; > 2) step-20 uses ApproximateSchurComplement because we actually don't > want to/can't assemble the

[deal.II] On using block diagonal preconditioner for the mixed Poisson problem

Dear all, I'm trying to use a block-diagonal preconditioner for the mixed Poisson problem from step-20, but haven't made it yet. The diagonal preconditioner of interest is from Powell & Silvester,2003

[deal.II] triangulation.add_periodicity takes an enormous amount of time

Dear all, I am trying to read from a mesh file and impose periodic boundary condition. This works out very well for small geometries in three spatial dimensions. However, if the read meshes tend to have more elements, the amount of time needed to create the mesh rises drastically. The time is

[deal.II] On using block-diagonal preconditioner for the mixed Poisson problem from step-20

Dear all, I'm trying to use a block-diagonal preconditioner for the mixed Poisson problem from step-20, but haven't made it yet. The diagonal preconditioner of interest is from Powell & Silvester,2003 . Specifically, for the matrix

Re: [deal.II] Re: Where can I find member function of cell?

Hi, Bruno, Thanks a lot! It really helps! On Tue, Aug 27, 2019 at 7:28 AM Bruno Turcksin wrote: > Yuesu, > > For a Triangulation, the list is here: > https://dealii.org/current/doxygen/deal.II/classCellAccessor-members.html > (CellAccessor) > For a DofHandler, the list is here: >

[deal.II] Re: Where can I find member function of cell?

Yuesu, For a Triangulation, the list is here:https://dealii.org/current/doxygen/deal.II/classCellAccessor-members.html (CellAccessor) For a DofHandler, the list is here: https://dealii.org/current/doxygen/deal.II/classDoFCellAccessor-members.html (DofCellAccessor) Best, Bruno On Monday,

[deal.II] VectorTools::Compress() fails after mesh_loop and distribute_local_to_global

I am trying to assemble my residual/RHS by using mesh_loop: evaluation_point = present_solution; extension_vector = newton_update; extension_vector *= alpha; evaluation_point += extension_vector; residual.reinit(locally_owned_dofs, mpi_communicator);

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug! On Tuesday, August 27, 2019 at 3:41:11 AM UTC+2, Doug wrote: > > Thank you very much for the quick fix! Looking forward to pull this once > it goes through all the checks. > The patch has been merged. Let me know if this fix does the trick for you. We introduced a test named