Add MacaulayNullspace constructor

JuliaAlgebra · Jul 2, 2024 · 69dc288 · 69dc288
1 parent 745c911
commit 69dc288
Showing 1 changed file with 11 additions and 0 deletions.
diff --git a/src/null.jl b/src/null.jl
@@ -32,6 +32,17 @@ function Base.getindex(
     )
 end
 
+function MacaulayNullspace(
+    ν::MomentMatrix,
+    rank_check::RankCheck,
+    ldlt::LowRankLDLTAlgorithm = SVDLDLT(),
+)
+    M = value_matrix(ν)
+    chol = low_rank_ldlt(M, ldlt, rank_check)
+    @assert size(chol.L, 1) == LinearAlgebra.checksquare(M)
+    return MacaulayNullspace(chol.L, ν.basis, accuracy(chol))
+end
+
 abstract type MacaulayNullspaceSolver end
 
 function solve(null::MacaulayNullspace, solver::MacaulayNullspaceSolver)