diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp b/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
index ba8ab4c497bcff1cc7daddbf4068c855b734b676..ef66ec6f20eebe3cf3a5386dcbbbe0a2711c79ca 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
@@ -465,10 +465,6 @@ void RichardsMechanicsProcess<DisplacementDim>::assembleConcreteProcess(
 {
     DBUG("Assemble the equations for RichardsMechanics");
 
-    // Note: This assembly function is for the Picard nonlinear solver. Since
-    // only the Newton-Raphson method is employed to simulate coupled HM
-    // processes in this class, this function is actually not used so far.
-
     std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
         dof_table = {std::ref(*_local_to_global_index_map)};
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];