Re: [deal.II] How to set material id with MPI

2019-08-29 Thread Wolfgang Bangerth
On 8/29/19 6:31 PM, Phạm Ngọc Kiên wrote: > > When I run the codes in my computer, it takes  a lot of time for p4est to > load > the grid. > The loading grid  step is more time consuming than solving the system of > equations with a mesh containing about 100,000 cells. It *shouldn't* take that

Re: [deal.II] How to set material id with MPI

2019-08-29 Thread Phạm Ngọc Kiên
Dear all, I think we have two ways to do this. The first one is the way Prof. Wolfgang Bangerth suggested. The second one is to load the grid in a Triangulation in all processor, then we set the material id before copying parallel::distributed::Triangulation from the Triangulation. When I run the

Re: [deal.II] Solving a TrilinosWrappers::BlockSparseMatrix using Trilinos iterative solver (step-57)

2019-08-29 Thread Daniel Arndt
Bruno, I am currently trying to work on a "HPC-ready with Trilinos" version of > step-57 for steady-state solution of the Navier-Stokes equation since I am > curious about compairing it to our stabilized solver. I think that it > could give much better results for steady-state. > The approach deta

Re: [deal.II] Implement 'w.n = u.n' as a boundary condition, where 'u' is a vector valued FE function/solution from a phase/PDE and 'n' is the unit outward normal in a distributed setting

2019-08-29 Thread Daniel Arndt
Bhanu, I have looked at 'VectorTools::compute_nonzero_normal_flux_constraints'. > But what I have here to compute the 'inhomogeneity' is not a closed form > function, but a discrete FE solution computed from another PDE(in my case > Stokes). The 'u' in 'w.n = u.n' is a 'LA::MPI::Vector' FE solutio

[deal.II] Solving a TrilinosWrappers::BlockSparseMatrix using Trilinos iterative solver (step-57)

2019-08-29 Thread Bruno Blais
Hello, I am currently trying to work on a "HPC-ready with Trilinos" version of step-57 for steady-state solution of the Navier-Stokes equation since I am curious about compairing it to our stabilized solver. I think that it could give much better results for steady-state. The approach detailed