From 12567167546e45226d0da6b6a002c317534ad862 Mon Sep 17 00:00:00 2001
From: Ben Knight <55677727+bknight1@users.noreply.github.com>
Date: Wed, 15 Jan 2025 11:45:09 +0800
Subject: [PATCH] fix pressure variable (#278)

Update petsc_generic_snes_solvers.pyx

Fixes #277 to is not yet optimised. More work to come on that.
---
 src/underworld3/cython/petsc_generic_snes_solvers.pyx | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/underworld3/cython/petsc_generic_snes_solvers.pyx b/src/underworld3/cython/petsc_generic_snes_solvers.pyx
index 887038a6..b4b9a47d 100644
--- a/src/underworld3/cython/petsc_generic_snes_solvers.pyx
+++ b/src/underworld3/cython/petsc_generic_snes_solvers.pyx
@@ -2721,7 +2721,7 @@ class SNES_Stokes_SaddlePt(SolverBaseClass):
                 dof = section.getFieldDof(p, field)
                 if dof > 0:
                     offset = section.getFieldOffset(p, field)
-                    if not unconstrained:
+                    if not unconstrained and self.Unknowns.p.continuous:
                         indices.append(offset)
                     else:
                         cind = section.getFieldConstraintIndices(p, field)