[aspect-devel] Cylindrical coordinates

[Jonathan Perry-Houts]

Thanks Wolfgang and Felipe. Sounds like there's no library magic that 
saves me in this case. Thanks for the advice!

-Jonathan

On 05/09/2018 08:11 PM, Wolfgang Bangerth wrote:
> 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.
> 

-- 
Jonathan Perry-Houts
Ph.D. Candidate
Department of Earth Sciences
1272 University of Oregon
Eugene, OR 97403-1272