Struct P1HarmonicCoordinateRelation 

pub struct P1HarmonicCoordinateRelation<const D: usize> { /* private fields */ }

The P1 Laplace relation defining one fixed-topology harmonic coordinate field.

The relation is sealed to an immutable reference mesh and an exact two-region cell partition. Solid/interface vertices are driven by the supplied solid displacement, fluid-only exterior vertices are fixed, and the remaining fluid vertices satisfy K_II u_I + K_ID u_D = 0 component by component. Construction assembles that relation but performs no solve.

Implementations

§

impl<const D: usize> P1HarmonicCoordinateRelation<D>

Signature links:usizeP1HarmonicCoordinateRelation

pub fn new( reference_mesh: &SimplicialMesh, fluid_cells: Vec<CellId>, solid_cells: Vec<CellId>, interface_facets: Vec<FacetId>, ) -> Result<P1HarmonicCoordinateRelation<D>, Diagnostic>

Signature links:newSimplicialMeshVecCellIdVecCellIdVecFacetIdResultP1HarmonicCoordinateRelationDiagnostic

Seal the solver-independent coordinate relation on one reference mesh.

Cell and facet inventories must be strictly increasing. Fluid and solid cells must cover every cell exactly once, and interface_facets must equal the complete set of cross-region facets. Solver tolerances and admission policy are intentionally absent from this geometry-level relation.

Errors

Returns EQ0803 unless the mesh is intrinsic 2D/3D, the partition and complete Dirichlet closure are exact and the bounded dense relation can be assembled from positive affine simplices.

pub const fn reference_mesh(&self) -> &SimplicialMesh

Signature links:reference_meshSimplicialMesh

Immutable mesh revision defining topology, reference coordinates, and coefficients.

pub fn fluid_cells(&self) -> &[CellId]

Signature links:fluid_cellsCellId

Fluid cells in canonical mesh order.

pub fn solid_cells(&self) -> &[CellId]

Signature links:solid_cellsCellId

Solid cells in canonical mesh order.

pub fn interface_facets(&self) -> &[FacetId]

Signature links:interface_facetsFacetId

Complete cross-region interface in canonical mesh-facet order.

pub fn solid_vertices(&self) -> &[VertexId]

Signature links:solid_verticesVertexId

Vertices driven exactly by the mesh-wide solid-displacement field.

pub fn driver_vertices(&self) -> &[VertexId]

Signature links:driver_verticesVertexId

Interface vertices supplying the fluid harmonic Dirichlet trace.

pub fn fixed_exterior_vertices(&self) -> &[VertexId]

Signature links:fixed_exterior_verticesVertexId

Fluid-only physical-exterior vertices fixed to reference coordinates.

pub fn fluid_interior_vertices(&self) -> &[VertexId]

Signature links:fluid_interior_verticesVertexId

Fluid vertices governed by the assembled interior Laplace relation.

pub fn interior_stiffness(&self) -> &[f64]

Signature links:interior_stiffnessf64

Row-major K_II in Self::fluid_interior_vertices order.

This solver-neutral lowered operator is exposed read-only so numerical realizations can solve the same relation later certified here.

pub fn driver_stiffness(&self) -> &[f64]

Signature links:driver_stiffnessf64

Row-major K_ID with interior rows and Self::driver_vertices columns.

pub fn driver_rhs_norms(&self) -> &[f64]

Signature links:driver_rhs_normsf64

Euclidean norm of each unit-driver right-hand side in driver-vertex order.

A solver-owning layer can combine these geometry-derived norms with its exact linear plan to derive the residual targets later supplied to Self::validate_current_coordinates.

pub fn validate_current_coordinates( &self, solid_displacement: &[[f64; D]], current_coordinates: &[Vec<f64>], driver_residual_targets: &[f64], ) -> Result<(), Diagnostic>

Signature links:validate_current_coordinatesf64DVecf64f64Result()Diagnostic

Certify candidate coordinates as the admitted harmonic extension.

solid_displacement is mesh-wide and must be exact zero outside the solid closure. Candidate solid coordinates must equal reference + displacement exactly, fixed fluid-exterior coordinates must equal the reference exactly, and fluid-interior coordinates must satisfy the P1 Laplace relation. The residual allowance is the triangle-inequality combination of the supplied per-driver solver targets plus a binary64 reduction roundoff term; no solver-produced coordinate or caller-authored residual is trusted. Targets are an external admission policy, not geometry state: artifact/numerics callers must derive and bind them from the exact realization plan or accepted solve reports.

Errors

Returns EQ0803 for incompatible/non-finite fields or targets, altered Dirichlet coordinates, overflow, or a harmonic residual outside the supplied solver bound.

Trait Implementations

§

impl<const D: usize> Clone for P1HarmonicCoordinateRelation<D>

Signature links:usizeCloneP1HarmonicCoordinateRelation

§

fn clone(&self) -> P1HarmonicCoordinateRelation<D>

Signature links:cloneP1HarmonicCoordinateRelation

Returns a duplicate of the value.

§

fn clone_from(&mut self, source: &Self)

Signature links:unstableSourceclone_from

Performs copy-assignment from source.

§

impl<const D: usize> Debug for P1HarmonicCoordinateRelation<D>

Signature links:usizeDebugP1HarmonicCoordinateRelation

§

fn fmt(&self, f: &mut Formatter<'_>) -> Result<(), Error>

Signature links:fmtFormatterResult()Error

Formats the value using the given formatter.

§

impl<const D: usize> PartialEq for P1HarmonicCoordinateRelation<D>

Signature links:usizePartialEqP1HarmonicCoordinateRelation

§

fn eq(&self, other: &P1HarmonicCoordinateRelation<D>) -> bool

Signature links:eqP1HarmonicCoordinateRelationbool

Tests for self and other values to be equal, and is used by ==.

§

fn ne(&self, other: &Rhs) -> bool

Signature links:unstableSourcene&Rhsbool

Tests for !=. The default implementation is almost always sufficient, and should not be overridden without very good reason.

impl<const D: usize> StructuralPartialEq for P1HarmonicCoordinateRelation<D>