From 7887c8d17442da6e68e2223309969fc8c3bdcb82 Mon Sep 17 00:00:00 2001
From: Francesc Verdugo <fverdugo@fs0.das6.cs.vu.nl>
Date: Sat, 18 May 2024 06:13:32 +0200
Subject: [PATCH] Sorting ghost col ids in own_ghost block according global ids

---
 src/ksp.jl | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/src/ksp.jl b/src/ksp.jl
index c40d884..f76b466 100644
--- a/src/ksp.jl
+++ b/src/ksp.jl
@@ -263,6 +263,20 @@ function MatCreateMPIAIJWithSplitArrays_args(a::PSparseMatrix,petsc_comm,cols=re
         oj .-= u
         m = own_length(rows)
         n = own_length(cols)
+        for i in 1:(length(oi)-1)
+            pini = oi[i]+1
+            pend = oi[i+1]
+            issorted = true
+            for p in pini:(pend-1)
+               if oj[p+1] < oj[p]
+                   issorted = false
+                   break
+               end
+            end
+            if ! issorted
+                sort!(view(oj,pini:pend))
+            end
+        end
         (petsc_comm,m,n,M,N,i,j,v,oi,oj,ov)
     end
     args = map(setup,partition(a),rows,cols)