[aspect-devel] Cylindrical coordinates

[Maxwell Rudolph]

If you decide to pursue this it might make sense to also implement the 2D spherical annulus simultaneously. 

> On May 9, 2018, at 10:16 PM, Jonathan Perry-Houts <jperryh2 at uoregon.edu> wrote:
> 
> 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
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel