In weakly formulated PDE (or even ODE), we seem to be interested in solutions that satisfy or take desired values at some boundary points of a domain we are interested in. For example, Dirichlet boundary conditions impose the function to take zero or certain value on the boundary of the domain, and so are other types like periodic boundary conditions and Neumann boundary conditions, Robin boundary conditions, etc.,

My question is, can we also impose condition like, not only the function staisfy one of the above mentioned conditions, but also a new condition that function should take a so and so values at so and so points (finite in number) that lie in the interior of the domain.

Writing down as equation, let $P = \{p_i\}$ be a finite set of points in the interior of the domain, and condition is that $f(p_i) = d_i$, $d_i\in\mathbb{R}$ $d_i$ are the desired values the function should take at the points $p_i$

I want to know the questions in the title of this post, for the above mentioned conditions. The above pointwise conditions make sense only for weakly formulated PDE, so I am not interested in strong PDE.

boundary pointsof the domain $\Omega\setminus \{p_i\}$. (Of course, they may or may not satisfy a given regularity property for the boundary) $\endgroup$