arXiv:2507.16964v2 Announce Type: replace Abstract: Solving partial differential equations (PDEs) on complex domains can present significant computational challenges. The Diffuse Domain Method (DDM) is an alternative that reformulates the partial differential equations on a larger, simpler domain. The original geometry is embedded into the problem by representing it with a phase-field function. This paper introduces DDFEM, an extensible Python package to provide a framework for transforming PDEs into a Diffuse Domain formulation. We aim to make the application of a variety of different Diffuse Domain approaches more accessible and straightforward to use. The DDFEM package includes features to intuitively define complex domains by combining signed distance functions and provides a number of Diffuse Domain transformers for general second-order evolution equations. In addition, we present a new approach for combining multiple boundary conditions of different types on distinct boundary segments. This is achieved by applying a normalised weighting, derived from the component signed distance functions, to combine the additional boundary terms in the Diffuse Domain formulations without requiring additional phase fields. The domain definition and Diffuse Domain transformation provided by our package are designed to be compatible with a wide range of existing finite element solvers without requiring code alterations. Both the original (non-linear) PDEs provided by the user and the resulting transformed PDEs on the extended domain are defined using the Unified Form Language (UFL). Here we focus on the software implementation and demonstrate its effectiveness across a range of problems.