An efficient regularization method for a large scale ill-posed geothermal problem

Numerical results, for 2D domains, show that the algorithm works well and the inverse problem can be solved accurately with a realistic noise level in the surface data.