Implementation of Point Source

Junxiong Jia, 2020.12.02

Introduction

For solving inverse problems, we often meet point source when tacking adjoint equations, e.g., the following equation

where are some measurement points, is usually some function obtained from forward equation. In the following, we explain that how to implement the point source term on the right-hand side of the above equation.

Measurement Matrix

Usually, for a function in ( is some open bounded domain), it can be decomposed by some fixed orthonormal basis . In the finite element method, we project functions on a finite number of basis as follow

And once we specify the coefficients , the function will be specified. We introduce the following measurement matrix

The matrix has the following property

Implement of the Point Source

Now we can give specific suggestions for implementing the point source. For discretization, we usually introduce the following quantities

with . Actually, we can reduce by some simple calculations as follow:

With these calculations, we find that

Denote , and (coefficients of the function ), we can reduce the above formula as follows:

which is convenient for implementations by some numerical algorithm, e.g., finite element method.