diff --git a/CHANGELOG.rst b/CHANGELOG.rst
index a941f13..a0e4363 100644
--- a/CHANGELOG.rst
+++ b/CHANGELOG.rst
@@ -3,6 +3,12 @@
 Next Release
 -----
 
+1.9.1
+-----
+* expressions for Constraint in the matrix interface are now assembled from only non-zero
+  coefficients, speeding up JSON conversion as a side effect
+* fixed a typo in the matrix interface
+
 1.9.0
 -----
 * enable Python 3.14
diff --git a/src/optlang/matrix_interface.py b/src/optlang/matrix_interface.py
index 4ad9b01..999438e 100644
--- a/src/optlang/matrix_interface.py
+++ b/src/optlang/matrix_interface.py
@@ -452,7 +452,11 @@ def _get_expression(self):
         if self.problem is not None:
             variables = self.problem._variables
             all_coefs = self.problem.problem.constraint_coefs
-            coefs = [(v, all_coefs.get((self.name, v.name), 0.0)) for v in variables]
+            coefs = [
+                (variables[vname], coef)
+                for (cname, vname), coef in all_coefs.items()
+                if cname == self.name
+            ]
             expression = add([mul((symbolics.Real(co), v)) for (v, co) in coefs])
             self._expression = expression
         return self._expression
@@ -823,7 +827,7 @@ def objective(self, value):
         (
             offset,
             linear_coefficients,
-            quadratic_coeffients,
+            quadratic_coefficients,
         ) = parse_optimization_expression(value, quadratic=True, expression=expression)
         self._objective_offset = offset
         if linear_coefficients:
@@ -831,7 +835,7 @@ def objective(self, value):
                 v.name: float(c) for v, c in linear_coefficients.items()
             }
 
-        for key, coef in quadratic_coeffients.items():
+        for key, coef in quadratic_coefficients.items():
             if len(key) == 1:
                 var = next(iter(key))
                 self.problem.obj_quadratic_coefs[(var.name, var.name)] = float(coef)