Re: [deal.II] Questions on Step-18

On 6/24/21 6:32 PM, Michael Li wrote: Andrew, thanks for confirming that. The missing 1/2 does not affect the demonstration of functionalities of deal.II but it may change the results. If this is wrong in the deal.II code, would you be interested in writing a patch to correct this? It should b

Re: [deal.II] Re: Subdivided cylinder and boundary ids

On 6/25/21 12:37 PM, Marc Fehling wrote: # Instead of an absolute threshold, you can introduce a relative one, for example, `epsilon * cell->diameter()` or something similar. FWIW, in most places, this is what we do. Best W. --

[deal.II] Re: Subdivided cylinder and boundary ids

By the way, this function has just been introduced in April with #12003 . Suggestions for improvement on this new feature are welcome! Best, Marc On Friday, June 25, 2021 at 12:37:31 PM UTC-6 Marc Fehling wrote: > Hi Alex, > > > I scrolled through g

[deal.II] Re: Subdivided cylinder and boundary ids

Hi Alex, > I scrolled through grid_generator.cc and found that subdivided_cylinder has a hardcoded epsilon of 1e-5 for determining whether a face is at the boundary. Of course, with a system of my size, all faces will end up being considered to be at the boundary. > > Because it not only sets t

RE: [deal.II] Questions on Step-18

Hi Jean-Paul,Thank you so much for the deliberate suggestions! It seems the code-gallery example is a good starting point for me to implement St-Venant Kirchhoff hyperelastic model with the help of “push-forward” transformations. I need to dig into this example to understand how the geometric stiff

[deal.II] Subdivided cylinder and boundary ids

I have been trying to construct a system that is small (a total length of 1e-5 m) with the subdivided cylinder grid generator. However, I have been having trouble with the boundary indicators. I wanted to set my own boundaries, and because there is no colorize option (it defaults to setting the