Inheritance diagram for aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >:

Public Member Functions
	SphericalManifoldWithTopography (const std::shared_ptr< const InitialTopographyModel::Interface< dim >> topography, const double inner_radius, const double outer_radius)

	SphericalManifoldWithTopography (const SphericalManifoldWithTopography< dim > &)=default

virtual std::unique_ptr< Manifold< dim, dim > >	clone () const override

Point< dim >	push_forward_from_sphere (const Point< dim > &p) const

Point< dim >	pull_back_to_sphere (const Point< dim > &p) const

virtual Point< dim >	get_intermediate_point (const Point< dim > &p1, const Point< dim > &p2, const double w) const override

virtual Tensor< 1, dim >	get_tangent_vector (const Point< dim > &x1, const Point< dim > &x2) const override

virtual Tensor< 1, dim >	normal_vector (const typename Triangulation< dim, dim >::face_iterator &face, const Point< dim > &p) const override

virtual void	get_normals_at_vertices (const typename Triangulation< dim, dim >::face_iterator &face, typename Manifold< dim, dim >::FaceVertexNormals &face_vertex_normals) const override

virtual void	get_new_points (const ArrayView< const Point< dim >> &surrounding_points, const Table< 2, double > &weights, ArrayView< Point< dim >> new_points) const override

virtual Point< dim >	get_new_point (const ArrayView< const Point< dim >> &vertices, const ArrayView< const double > &weights) const override

double	topography_for_point (const Point< dim > &x_y_z) const

Public Member Functions inherited from aspect::SphericalManifold< dim, spacedim >
	SphericalManifold (const Point< spacedim > center=Point< spacedim >())

virtual Point< spacedim >	get_intermediate_point (const Point< spacedim > &p1, const Point< spacedim > &p2, const double w) const override

virtual Tensor< 1, spacedim >	get_tangent_vector (const Point< spacedim > &x1, const Point< spacedim > &x2) const override

virtual Tensor< 1, spacedim >	normal_vector (const typename Triangulation< dim, spacedim >::face_iterator &face, const Point< spacedim > &p) const override

virtual void	get_normals_at_vertices (const typename Triangulation< dim, spacedim >::face_iterator &face, typename Manifold< dim, spacedim >::FaceVertexNormals &face_vertex_normals) const override

virtual void	get_new_points (const ArrayView< const Point< spacedim >> &surrounding_points, const Table< 2, double > &weights, ArrayView< Point< spacedim >> new_points) const override

virtual Point< spacedim >	get_new_point (const ArrayView< const Point< spacedim >> &vertices, const ArrayView< const double > &weights) const override

Public Member Functions inherited from Manifold< dim, dim >
virtual	~Manifold () override=default

void	serialize (Archive &ar, const unsigned int version)

virtual Point< spacedim >	get_new_point (const ArrayView< const Point< spacedim > > &surrounding_points, const ArrayView< const double > &weights) const

virtual Point< spacedim >	get_new_point (const ArrayView< const Point< spacedim > > &surrounding_points, const ArrayView< const double > &weights) const

virtual void	get_new_points (const ArrayView< const Point< spacedim > > &surrounding_points, const Table< 2, double > &weights, ArrayView< Point< spacedim > > new_points) const

virtual void	get_new_points (const ArrayView< const Point< spacedim > > &surrounding_points, const Table< 2, double > &weights, ArrayView< Point< spacedim > > new_points) const

virtual Point< spacedim >	project_to_manifold (const ArrayView< const Point< spacedim > > &surrounding_points, const Point< spacedim > &candidate) const

virtual Point< spacedim >	project_to_manifold (const ArrayView< const Point< spacedim > > &surrounding_points, const Point< spacedim > &candidate) const

virtual Point< spacedim >	get_new_point_on_line (const typename Triangulation< dim, spacedim >::line_iterator &line) const

virtual Point< spacedim >	get_new_point_on_line (const typename Triangulation< dim, spacedim >::line_iterator &line) const

virtual Point< spacedim >	get_new_point_on_quad (const typename Triangulation< dim, spacedim >::quad_iterator &quad) const

virtual Point< spacedim >	get_new_point_on_quad (const typename Triangulation< dim, spacedim >::quad_iterator &quad) const

virtual Point< spacedim >	get_new_point_on_hex (const typename Triangulation< dim, spacedim >::hex_iterator &hex) const

virtual Point< spacedim >	get_new_point_on_hex (const typename Triangulation< dim, spacedim >::hex_iterator &hex) const

Point< spacedim >	get_new_point_on_face (const typename Triangulation< dim, spacedim >::face_iterator &face) const

Point< spacedim >	get_new_point_on_face (const typename Triangulation< dim, spacedim >::face_iterator &face) const

Point< spacedim >	get_new_point_on_cell (const typename Triangulation< dim, spacedim >::cell_iterator &cell) const

Point< spacedim >	get_new_point_on_cell (const typename Triangulation< dim, spacedim >::cell_iterator &cell) const

virtual void	get_normals_at_vertices (const typename Triangulation< dim, spacedim >::face_iterator &face, FaceVertexNormals &face_vertex_normals) const

virtual void	get_normals_at_vertices (const typename Triangulation< dim, spacedim >::face_iterator &face, FaceVertexNormals &face_vertex_normals) const

unsigned int	n_subscriptions () const

unsigned int	n_subscriptions () const

void	list_subscribers (StreamType &stream) const

void	list_subscribers () const

void	list_subscribers (StreamType &stream) const

void	list_subscribers () const

Public Member Functions inherited from EnableObserverPointer
	EnableObserverPointer ()

	EnableObserverPointer (const EnableObserverPointer &)

	EnableObserverPointer (EnableObserverPointer &&) noexcept

virtual	~EnableObserverPointer ()

EnableObserverPointer &	operator= (const EnableObserverPointer &)

EnableObserverPointer &	operator= (EnableObserverPointer &&) noexcept

void	serialize (Archive &ar, const unsigned int version)

unsigned int	n_subscriptions () const

void	list_subscribers (StreamType &stream) const

void	list_subscribers () const

unsigned int	n_subscriptions () const

void	list_subscribers (StreamType &stream) const

void	list_subscribers () const

Private Attributes
const std::shared_ptr< const InitialTopographyModel::Interface< dim > >	topo

const double	R0

const double	R1

Additional Inherited Members
Public Types inherited from Manifold< dim, dim >
typedef std::array< Tensor< 1, spacedim >, GeometryInfo< dim >::vertices_per_face >	FaceVertexNormals

Static Public Member Functions inherited from Manifold< dim, dim >
static ::ExceptionBase &	ExcInUse (int arg1, std::string arg2, std::string arg3)

static ::ExceptionBase &	ExcNoSubscriber (std::string arg1, std::string arg2)

Static Public Member Functions inherited from EnableObserverPointer
static ::ExceptionBase &	ExcInUse (int arg1, std::string arg2, std::string arg3)

static ::ExceptionBase &	ExcNoSubscriber (std::string arg1, std::string arg2)

Public Attributes inherited from aspect::SphericalManifold< dim, spacedim >
const Point< spacedim >	center

Detailed Description

template<int dim>
class aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >

A description of a manifold that describes a spherical shell with overlaid topography.

Definition at line 39 of file spherical_shell.h.

Constructor & Destructor Documentation

◆ SphericalManifoldWithTopography() [1/2]

template<int dim>

aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::SphericalManifoldWithTopography	(	const std::shared_ptr< const InitialTopographyModel::Interface< dim >>	topography,
		const double	inner_radius,
		const double	outer_radius
	)

Constructor.

◆ SphericalManifoldWithTopography() [2/2]

template<int dim>

aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::SphericalManifoldWithTopography ( const SphericalManifoldWithTopography< dim > & )

default

Copy constructor.

Member Function Documentation

◆ clone()

template<int dim>

virtual std::unique_ptr<Manifold<dim, dim> > aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::clone ( ) const

overridevirtual

Make a clone of this Manifold object.

Reimplemented from aspect::SphericalManifold< dim, spacedim >.

◆ push_forward_from_sphere()

template<int dim>

Point<dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::push_forward_from_sphere ( const Point< dim > & p ) const

Given a point in the undeformed spherical geometry, push it forward to the corresponding point in the sphere with surface topography.

◆ pull_back_to_sphere()

template<int dim>

Point<dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::pull_back_to_sphere ( const Point< dim > & p ) const

Given a point in the deformed spherical geometry with topography, pull it back to the corresponding point in the undeformed sphere.

◆ get_intermediate_point()

template<int dim>

virtual Point<dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::get_intermediate_point	(	const Point< dim > &	p1,
		const Point< dim > &	p2,
		const double	w
	)		const

overridevirtual

Given any two points in space, first project them on the surface of a sphere with unit radius, then connect them with a geodesic and find the intermediate point, and finally rescale the final radius so that the resulting one is the convex combination of the starting radii.

◆ get_tangent_vector()

template<int dim>

virtual Tensor<1, dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::get_tangent_vector	(	const Point< dim > &	x1,
		const Point< dim > &	x2
	)		const

overridevirtual

Compute the derivative of the get_intermediate_point() function with parameter w equal to zero.

◆ normal_vector()

template<int dim>

virtual Tensor<1, dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::normal_vector	(	const typename Triangulation< dim, dim >::face_iterator &	face,
		const Point< dim > &	p
	)		const

overridevirtual

We fudge here, but for a good reason. What the function is supposed to compute is the normal vector to the surface. This should be the normal to the surface with topography, but instead we return the normal to the undeformed surface – i.e., the radial direction. This is, in particular, used to compute no-flux boundary conditions, for which we want to impose a boundary condition that allows for plate-like motion – that is, we need to allow horizontal motion, even if that is not tangential to the surface along the slopes of mountains or ocean trenches. Using the radial direction, i.e., the normal vector to the undeformed surface (= a radial vector) allows for exactly this.

◆ get_normals_at_vertices()

template<int dim>

virtual void aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::get_normals_at_vertices	(	const typename Triangulation< dim, dim >::face_iterator &	face,
		typename Manifold< dim, dim >::FaceVertexNormals &	face_vertex_normals
	)		const

overridevirtual

Compute the normal vectors to the boundary at each vertex.

◆ get_new_points()

template<int dim>

virtual void aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::get_new_points	(	const ArrayView< const Point< dim >> &	surrounding_points,
		const Table< 2, double > &	weights,
		ArrayView< Point< dim >>	new_points
	)		const

overridevirtual

Compute a new set of points that interpolate between the given points surrounding_points. weights is a table with as many columns as surrounding_points.size(). The number of rows in weights must match the length of new_points.

This function is optimized to perform on a collection of new points, by collecting operations that are not dependent on the weights outside of the loop over all new points.

The implementation does not allow for surrounding_points and new_points to point to the same array, so make sure to pass different objects into the function.

◆ get_new_point()

template<int dim>

virtual Point<dim> aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::get_new_point	(	const ArrayView< const Point< dim >> &	vertices,
		const ArrayView< const double > &	weights
	)		const

overridevirtual

Return a point on the spherical manifold which is intermediate with respect to the surrounding points.

◆ topography_for_point()

template<int dim>

double aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::topography_for_point ( const Point< dim > & x_y_z ) const

Return the topography of the surface directly above the point given by the coordinates stored in the argument.

Member Data Documentation

◆ topo

template<int dim>

const std::shared_ptr<const InitialTopographyModel::Interface<dim> > aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::topo

private

A pointer to the topography model.

Definition at line 160 of file spherical_shell.h.

◆ R0

template<int dim>

const double aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::R0

private

Inner and outer radii of the spherical shell.

Definition at line 165 of file spherical_shell.h.

◆ R1

template<int dim>

const double aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >::R1

private

Definition at line 165 of file spherical_shell.h.

The documentation for this class was generated from the following file:

/home/docs/checkouts/readthedocs.org/user_builds/aspect-documentation/checkouts/latest/include/aspect/geometry_model/spherical_shell.h

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<int dim> class aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >

Constructor & Destructor Documentation

◆ SphericalManifoldWithTopography() [1/2]

◆ SphericalManifoldWithTopography() [2/2]

Member Function Documentation

◆ clone()

◆ push_forward_from_sphere()

◆ pull_back_to_sphere()

◆ get_intermediate_point()

◆ get_tangent_vector()

◆ normal_vector()

◆ get_normals_at_vertices()

◆ get_new_points()

◆ get_new_point()

◆ topography_for_point()

Member Data Documentation

◆ topo

◆ R0

◆ R1

template<int dim>
class aspect::GeometryModel::internal::SphericalManifoldWithTopography< dim >