Re: [deal.II] How to refine mesh in different field

Hi Chen, step-46 gives you a good example on how to perform grid refinement when coupling different kinds of equations in different parts of the domain. I would highly suggest reading this t

Re: [deal.II] On mesh refinement & solution transfer with Raviart-Thomas

Charlie, I don't have time in the next few days to look into the exact cause, but I have a suspicion and if you are up for some digging, here are a few pointers. When you transfer the solution from a cell to its children, then SolutionTransfer what are called the "prolongation matrices" that

[deal.II] On mesh refinement & solution transfer with Raviart-Thomas

Hi there, I've been playing with mesh refinement and solution transfer with the Raviart-Thomas field, but found a very strange thing: in areas where the mesh is coarsened, the solution becomes wrong. In a very simple 2D test, I first created a uniform RT field, i.e. (1, 0), with the help of `V

[deal.II] CUDA without nVidia

Dear all, Have anyone managed to compile and *run* CUDA code on non-nVidia device? I've found HIPify: https://github.com/ROCm-Developer-Tools/HIPIFY that translates CUDA code to HIP code that can be then compiled and ran on both AMD and nVidia GPUs. The thing is, I'm looking for a new laptop.

Re: [deal.II] How to refine mesh in different field

Jest overwrite error estimate in one domain by zeros: unsigned int cell_n = 0; for (const auto & cell : triangulation.active_cell_iterators ){ if (cell-> material_id() == 0 ) estimated_error_per_cell(i) = 0; ++i; } sobota, 3 kwietnia 2021 o 16:48:48 UTC+2 hkch...@gmail.com napisał

[deal.II] How to refine mesh in different field

Dear deal.II community, In a fluid-structure interaction simulation, Computational domain was divided into fluid domain and solid domain, for example, cell->material_id() == 0 and cell->material_id() == 1(the way that deal.ii presents which field the cell belong to) respectively, Now I want to ref