Compose shape-to-surface and surface-to-volume maps with the complete derivative chain rule. More...

#include <composition.hpp>

Collaboration diagram for serif::geometry::deformation::ComposedDeformation< SurfaceMap, ExtensionMap >:

Classes
struct	CurvatureWorkspace
	Nonlinear actions retain the tangent and the two contributions to the shape-force derivative. More...

struct	Workspace

Public Types
using	InputSpace = typename SurfaceMap::InputSpace

using	OutputSpace = typename ExtensionMap::OutputSpace

Public Member Functions
std::size_t	input_size () const noexcept
	Keep solver coordinates in the surface model's native parameter space.

std::size_t	output_size () const noexcept
	Expose the displacement field size consumed by the existing geometry mapper.

Workspace	make_workspace () const
	Reserve every intermediate once, outside repeated Newton and Krylov actions.

DeformationStatus	apply (std::span< const double > input, std::span< double > output, Workspace &workspace) const
	A fixed composition uses the same primal action as its Jacobian.

DeformationStatus	apply_transpose (std::span< const double > dual, std::span< double > output, Workspace &workspace) const
	Reverse a fixed composition using one intermediate and no artificial base state.

DeformationStatus	build_displacement (std::span< const double > parameters, std::span< double > displacement, Workspace &workspace) const
	Move the surface first, then extend that same vector trace into the volume.

DeformationStatus	apply_jacobian (std::span< const double > parameters, std::span< const double > direction, std::span< double > displacement_direction, Workspace &workspace) const
	Evaluate the tangent of each nonlinear factor at the accepted surface displacement.

DeformationStatus	apply_jacobian_transpose (std::span< const double > parameters, std::span< const double > displacement_dual, std::span< double > parameter_dual, Workspace &workspace) const
	Reverse both factors in the forward composition, preserving their work pairing.

DeformationStatus	apply_pullback_derivative (std::span< const double > parameters, std::span< const double > direction, std::span< const double > displacement_dual, std::span< double > parameter_dual, Workspace &workspace) const
	Retain both the extension curvature and the derivative of the surface transpose.

const SurfaceMap &	surface_map () const noexcept
	Expose immutable surface mathematics for reference metrics and diagnostics.

Static Public Attributes
static constexpr bool	is_linear = is_linear_deformation<SurfaceMap> && is_linear_deformation<ExtensionMap>

static constexpr bool	uses_linear_actions = LinearDeformationMap<SurfaceMap> && LinearDeformationMap<ExtensionMap>

Private Member Functions
	ComposedDeformation (std::shared_ptr< const SurfaceMap > surface, std::shared_ptr< const ExtensionMap > extension)
	Only the checked composition factory may establish factor compatibility.

Private Attributes
std::shared_ptr< const SurfaceMap >	m_surface

std::shared_ptr< const ExtensionMap >	m_extension

Friends
template<DeformationMap First, DeformationMap Second> requires std::same_as<typename First::OutputSpace, typename Second::InputSpace>
std::expected< ComposedDeformation< First, Second >, DeformationFailure >	compose_deformation (std::shared_ptr< const First >, std::shared_ptr< const Second >)

Detailed Description

template<DeformationMap SurfaceMap, DeformationMap ExtensionMap>
requires std::same_as<typename SurfaceMap::OutputSpace, typename ExtensionMap::InputSpace>
class serif::geometry::deformation::ComposedDeformation< SurfaceMap, ExtensionMap >

Compose shape-to-surface and surface-to-volume maps with the complete derivative chain rule.

Note: Numerical factors are immutable and owned. Workspaces belong to callers and must not overlap public inputs or outputs; this class has no preparation state.