[aspect-devel] Cylindrical coordinates

[Wolfgang Bangerth]

On 05/10/2018 10:16 AM, Jonathan Perry-Houts wrote:
> 
> Does anybody have a sense of how difficult it would be to implement 2D
> cylindrical coordinates in Aspect? It doesn't look like there's support
> for non-Cartesian derivatives in Deal.II, but maybe somebody has a
> suggestion for hacking something together.

You'd have to rewrite the assemblers to implement the bilinear form that 
describes the equation in cylindrical coordinates. You'd also have to hack all 
of the other places where we do any integrals (nullspace removal, postprocessors).

It's likely going to be a substantial amount of work -- not difficult, just 
going through it all.

Best
  W.

-- 
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bangerth at colostate.edu
                            www: http://www.math.colostate.edu/~bangerth/