Fix bugs in solid assemblage solver

Author: Kastakin

src/libeq/solver/nr.py

@@ -199,9 +199,9 @@ def _panic_save():
             np.abs(F) < threshold
         ):
             if solids_to_remove.size != 0:
-                x = np.insert(
-                    x, np.clip(solids_to_remove, a_min=0, a_max=x.shape[1]), 0, axis=1
-                )
+                temp = np.zeros((x.shape[0], n_components + n_solids))
+                temp[:, ~np.isin(range(n_components + n_solids), solids_to_remove)] = x
+                x = temp
             return x
 
         if damping:

src/libeq/solver/solids_solver.py

@@ -164,7 +164,7 @@ def _update_solids_set(
 
     """
     adjust_solids = True
-    negative_solid_concentration = c[:, point_log_ks.size :] < 0
+    negative_solid_concentration = c[:, -point_log_ks.size :] < 0
     supersaturated_solid = saturation_index > 1 + 1e-9
 
     # any negative solid concentration remove them from the solids set