La correttezza della soluzione molto spesso controllata verificando che il principio di conservazione della massa venga rispettato. Per questo si deve calcolare la massa di acqua che entra nel dominio e verificare che questa sia uguale alla massa d'acqua che esce dal dominio stesso.
La massa d'acqua entrata o uscita dai nodi di Dirichlet viene
calcolata semplicemente moltiplicando per il vettore soluzione
la matrice di rigidezza ,
come calcolata prima dell'imposizione delle condizioni al contorno.
Per fare ciø bisogna salvare i coefficienti diagonali della matrice
prima che essi siano cambiati dall'imposizione delle condizioni al
contorno.
La massa d'acqua entrata o uscita dai nodi di Neumann data direttamente dal valore della portata in quei nodi.
Il bilancio di massa viene eseguito sommando algebricamente i valori ottenuti per ciascun nodo di Dirichlet e di Neumann. Il risultato deve essere teoricamente nullo. In pratica il bilancio di massa si considera corretto se tale risultato è uguale ad una frazione piccola della massa entrante o uscente.