Don't ccall arb_mat_solve_lu_precomp

Nemocas · Oct 23, 2024 · c05ecf2 · c05ecf2
1 parent 5d18487
commit c05ecf2
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 6 deletions.
diff --git a/src/arb/RealMat.jl b/src/arb/RealMat.jl
@@ -584,9 +584,7 @@ function Solve._can_solve_internal_no_check(::Solve.LUTrait, C::Solve.SolveCtx{R
   end
 
   x = similar(b, ncols(C), ncols(b))
-  ccall((:arb_mat_solve_lu_precomp, libflint), Nothing,
-        (Ref{RealMatrix}, Ptr{Int}, Ref{RealMatrix}, Ref{RealMatrix}, Int),
-        x, inv(p).d .- 1, LU, b, precision(Balls))
+  _solve_lu_precomp!(x, p, LU, b)
 
   if side === :left
     x = transpose(x)

diff --git a/src/arb/arb_mat.jl b/src/arb/arb_mat.jl
@@ -611,9 +611,7 @@ function Solve._can_solve_internal_no_check(::Solve.LUTrait, C::Solve.SolveCtx{A
   end
 
   x = similar(b, ncols(C), ncols(b))
-  ccall((:arb_mat_solve_lu_precomp, libflint), Nothing,
-        (Ref{ArbMatrix}, Ptr{Int}, Ref{ArbMatrix}, Ref{ArbMatrix}, Int),
-        x, inv(p).d .- 1, LU, b, precision(base_ring(LU)))
+  _solve_lu_precomp!(x, p, LU, b)
 
   if side === :left
     x = transpose(x)