From 4fde98912ea95ca787891ec6db0ada5419a4427f Mon Sep 17 00:00:00 2001
From: Emanuele Sgroi <emanuele.sgroi@gmail.com>
Date: Thu, 12 Sep 2024 12:18:55 +0200
Subject: [PATCH] Update VNumber.m2

We made a minor change to the function "soc"
---
 M2/Macaulay2/packages/VNumber.m2 | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/M2/Macaulay2/packages/VNumber.m2 b/M2/Macaulay2/packages/VNumber.m2
index 482670b4adf..4ff9acf19dc 100644
--- a/M2/Macaulay2/packages/VNumber.m2
+++ b/M2/Macaulay2/packages/VNumber.m2
@@ -3,7 +3,7 @@
 newPackage(
     "VNumber",
     Version => "1.0",
-    Date => "October 24, 2023",
+    Date => "September 12, 2024",
     Authors => {{Name => "Antonino Ficarra", Email => "antficarra@unime.it", HomePage => "https://www.researchgate.net/profile/Antonino-Ficarra"},
         {Name => "Emanuele Sgroi", Email => "emasgroi@unime.it", HomePage => "https://www.researchgate.net/profile/Emanuele-Sgroi"}},
     Headline => "compute v-number of homogeneous ideals and v-function of monomial ideals",
@@ -376,12 +376,12 @@ soc (Ideal, Ideal) := (I, P) -> (
   XP := ring I;
   if not(isMember(true, (stableMax(I)/(Q->P==Q)))) then XP = product(select(stablePrimes(I), Q -> (Q!=P and isSubset(P,Q)==true)));
   f := reesMap(I);
-  R := source f/(sub(I, source f) + ker f);
+  R := source f/ker f;
   IR := sub(I,R);
   PR := sub(P,R);
-  QR := if not(XP===ring I) then sub(P + saturate XP, R) else ideal(1_R);
+  QR := if not(XP===ring I) then sub(XP, R) else ideal(1_R);
   A := quotient(IR, PR);
-  B := quotient(IR, QR);
+  B := intersect(quotient(IR, PR), saturate(IR, QR));
   A/B
 );