Hello FreeFem community,

I started using FreeFEM. I find it very easy to understand and use. I started liking it and considering to integrate it into my research workflow. I am looking for some clarification regarding the application of boundary conditions.

I read in the FreeFEM documentation that Dirichlet boundary conditions are imposed using the penalty method using a very high value (10^30) for the penalty constant.

My experience indicates that such an approach results in ill-conditioned matrices and poses many issues in the modal analysis (eigenvalue problems) and buckling analysis due to the sensitivity of results to the penalty constant.

I am wondering you have any plans for implementing the elimination approach for imposing Dirichlet BCs in FreeFEM in the future.

Thanks,

Chenna