On 6/11/19 1:06 AM, SebG wrote: > > actually I was not aware of this class. But this class is unfortunately only > implemented in 2D. I am not exactly looking for an ellipsoid but for a sphere > with a peturbed radius, where r = R * (1 + e * Y^m_n(\theta, \phi) ). Here R > is the reference radius, e > 0 is a scalar and Y^m_n is a real-valued > spherical harmonic. > > To implement this, I thought that I could map a point from the perturbed > sphere to the reference sphere. Then the updated point is computed on the > reference sphere and mapped back to the perturbed one. > > As a first step, I tried this approach based on the step-53 tutorial. But I > realized that closing the sphere at the poles and in azimuthal direction does > not work. Then I thought that it might easier to program the example of an > ellipsoid from this discussion. > > It would interessting to know why my approach of mapping back to a reference > sphere is not working with the updated version of SphericalManifold. But > maybe > this is worth a new discussion.
I don't know what the issue is, but the way I would write this is a concatenation of two operations. The first would be the transformation from reference coordinates to the sphere (done by the push_forward() function of the SphericalManifold) and the second the scaling in radial direction (done by your own code). step-53 does this same kind of 2-step procedure, but instead of using the WGS84 transformation used there, you'd just defer to SphericalManifold::push_forward. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/bda8d017-7dc4-1413-6ba2-8b1cf10d8dfd%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
