[FreeFem++] fonction indicatrice

Bonjour,

Je débute avec FreeFem et je suis confronté au problème suivant :

j'ai créé un domaine (rectangulaire) sur lequel j'ai généré un maillage.
Je cherche maintenant à construire une fonction (func) qui soit la somme de fonctions indicatrices valant 1 sur des disques contenus dans ce domaine. Mais pour autant, ces disques en question n'apparaissent pas du tout dans la définition de mon domaine et de mon maillage!
J'ai essayé plusieurs choses mais à chaque fois sans succès.
Cette fonction serait en fait utilisée dans la définition de mon problème sous forme faible comme un terme de source par exemple.

Si quelqu'un pouvait m'aider ce serait super.
Merci!