Customization

Medusa library support users defining custom basis types, weights, operators and more, as long as they conform to the prescribed interfaces, given in the Concepts page. Here we show this on an example of the biharmonic equation.

Custom operators: Biharmonic equation

Biharmonic operator $\nabla^4$ is not included in Medusa by default, but it is easy to define.

We solve the problem

$ \begin{align}

 \nabla^4 u &= f      &&\text{in } \Omega, \\
 u &= g_d           &&\text{on } \partial \Omega,\\
 \frac{\partial u}{\partial \vec n} &= g_n           &&\text{on } \partial \Omega, 

\end{align} $ where $u = TODO$ and $f$, $g_d$ and $g_n$ are computed from $u$.

TODO: JureMB