Skip to content

Implement pointobservation - #61

Open
finsberg wants to merge 12 commits into
mainfrom
finberg/pointobservation-by-hand
Open

Implement pointobservation#61
finsberg wants to merge 12 commits into
mainfrom
finberg/pointobservation-by-hand

Conversation

@finsberg

@finsberg finsberg commented Aug 14, 2026

Copy link
Copy Markdown
Member

I this PR we add PointObersevation class which builds the operator $B$, such that

$$ d = B u, \qquad B \in \mathbb{R}^{n_d \times N}, \qquad B_{ij} = \phi_j(x_i), $$

as well as a new function point_obersvation_misfit which build the following functional

$$ J(u) = \frac{1}{2\sigma^2} \lVert W (B u - d) \rVert^2 . $$

Given some poinst where you have your obersverations points, the operatror B is created as

B = dolfinx_adjoint.PointObservation(V,  points)

where V is the relevant finite element space.

AI assistance

I used Claude Code (Claude Sonnet 5) to help profile, implement, and iterate on several of these changes. I reviewed, tested, and take responsibility for the final contribution.

@finsberg
finsberg marked this pull request as ready for review August 19, 2026 08:45
@finsberg
finsberg requested a review from jorgensd August 19, 2026 08:45
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