Skip to content

Conversation

@kmokstad
Copy link
Contributor

@kmokstad kmokstad commented Nov 9, 2025

This implements general point searching for unstructured 2D Lagrange meshes, for use in result point definition. It reuses some old Fortran codes of the 90s (the AFEM project) to conduct explicit parametric inversion for 3- and 4-noded shell elements.

@kmokstad
Copy link
Contributor Author

kmokstad commented Nov 9, 2025

jenkins build this with downstreams ifem-shellex=40 please

@kmokstad kmokstad force-pushed the point-search branch 2 times, most recently from fd89fd0 to e428aa1 Compare November 10, 2025 06:39
@kmokstad
Copy link
Contributor Author

jenkins build this with downstreams please

Changed: Use RealArray instead of Vector when no algebraic operations.
to reduce repetive code. The read() call will not fail anyway
(and if it does, it will be caugth by the subsequent checks).
Also merge methods genThreadGroups() and getThreadGroups().
Use loops instead of lambdas in the tests for simplicity.
Fixed: Swap node 3 and 4 when reading 4-noded element topology.
@kmokstad
Copy link
Contributor Author

jenkins build this with downstreams ifem-shellex=41 please

spatial points. The default implementation invokes ASMbase::findPoint()
in a loop over the specified points. It is overridden in ASMu2DLag to
perform parametric inversion for bilinear quads and triangles.
Tighten the comparison with zero tolerance (1e-8 is too loose).
nodal field at specified points defined by element index and
local xi,eta parameters within those element.
Added: private method evalPt() evaluating a nodal field at a point.
@kmokstad
Copy link
Contributor Author

jenkins build this with downstreams ifem-shellex=41 please

Added: virtual method ASMbase::printElmInfo().
Added: Evaluation of 2ndary solution at specified Cartesian points.
Changed: Override evalSolPt() in ASMu2DLag with the triangle support.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant