dummy instead of Cii compress

test/compress.jl

@@ -3,7 +3,7 @@
     n_big = 100
     n_lit = 2
 
-    A = 1e-5*abs.(randn(n_big))
+    A = 1e-4*abs.(randn(n_big))
     V = randn(n_big,n_lit)
     W = DiagWoodbury(A,V)
 
@@ -12,31 +12,28 @@
     W1 = read_krybury(fname)
     @test W == W1
 
-    function Cii_precomp_mult(matList)
+    function dummy_precomp_mult(matList)
         return []
     end
 
-    function Cii_fxn_mult(matList,precompList,x)
-        Ctotinv = matList[1]
+    function dummy_fxn_mult(matList,precompList,x)
+        A = matList[1]
         Vi = matList[2]
-        arg2 = Vi*(Vi'*x)
-        return arg2 - Vi*(Vi'*(Ctotinv*arg2))
+        return A*x + Vi*(Vi'*x)
     end
 
-    function Cii_precomp_diag(matList)
-        Ctotinv = matList[1]
-        Vi = matList[2]
-        return [Vi'*(Ctotinv*Vi)]
+    function dummy_precomp_diag(matList)
+        return []
     end
 
-    function Cii_diag_map(matList,precompList)
+    function dummy_diag_map(matList,precompList)
+        A = matList[1]
         Vi = matList[2]
-        arg1 = precompList[1]
-        return dropdims(sum(Vi.^2,dims=2),dims=2).-dropdims(sum(Vi'.*(arg1*Vi'),dims=1),dims=1)
+        return diag(A) + dropdims(sum(Vi.^2,dims=2),dims=2)
     end
 
-    BMat = LowRankMultMat([Diagonal(A),V],Cii_precomp_mult,Cii_fxn_mult);
-    BMatDiag = LowRankDiagMat([Diagonal(A),V],Cii_precomp_diag,Cii_diag_map);
+    BMat = LowRankMultMat([Diagonal(A),V],dummy_precomp_mult,dummy_fxn_mult);
+    BMatDiag = LowRankDiagMat([Diagonal(A),V],dummy_precomp_diag,dummy_diag_map);
 
     out = kryburyCompress(BMat,BMatDiag,n_big)