diff --git a/test/test_shell_resultants.jl b/test/test_shell_resultants.jl
index 147f9b5..b0584e4 100644
--- a/test/test_shell_resultants.jl
+++ b/test/test_shell_resultants.jl
@@ -8,6 +8,7 @@ module mcompostat0
 
 using LinearAlgebra
 using FinEtools
+using FinEtools.AlgoBaseModule: solve_blocked!
 using FinEtoolsDeforLinear
 using FinEtoolsFlexStructures.FESetShellT3Module: FESetShellT3
 using FinEtoolsFlexStructures.FESetShellQ4Module: FESetShellQ4
@@ -107,8 +108,7 @@ function test_homogeneous(visualize = true)
     F = distribloads(lfemm, geom0, dchi, fi, 3)
 
     # Solve
-    U = K \ F
-    scattersysvec!(dchi, U[:])
+    solve_blocked!(dchi, K, F)
     result = dchi.values[nl, 3][1]
     # @info "Solution: $(result), $(round(result/analyt_sol*100, digits = 4))%"
 
@@ -289,8 +289,7 @@ function test_composite(visualize = true)
     F = distribloads(lfemm, geom0, dchi, fi, 3)
 
     # Solve
-    U = K \ F
-    scattersysvec!(dchi, U[:])
+    solve_blocked!(dchi, K, F)
     result = dchi.values[nl, 3][1]
     # @info "Solution: $(result), $(round(result/analyt_sol*100, digits = 4))%"