diff --git a/src/polysolve/linear/FEMSolver.cpp b/src/polysolve/linear/FEMSolver.cpp
index 5d1f8f8f..edbcd909 100644
--- a/src/polysolve/linear/FEMSolver.cpp
+++ b/src/polysolve/linear/FEMSolver.cpp
@@ -103,7 +103,6 @@ namespace polysolve::linear
         const bool remove_zero_cols,
         const bool skip_last_cols)
     {
-        assert(!solver.is_dense());
         // Let Γ be the set of Dirichlet dofs.
         // To implement nonzero Dirichlet boundary conditions, we seek to replace
         // the linear system Au = f with a new system Ãx = g, where
@@ -342,4 +341,4 @@ namespace polysolve::linear
         f = g;
     }
 
-} // namespace polysolve::linear
\ No newline at end of file
+} // namespace polysolve::linear