scribbles for cg method

AlgebraicJulia · Sep 25, 2024 · 029b7cf · 029b7cf
1 parent 0654a71
commit 029b7cf
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 0 deletions.
diff --git a/Project.toml b/Project.toml
@@ -14,9 +14,11 @@ Catlab = "134e5e36-593f-5add-ad60-77f754baafbe"
 DataMigrations = "0c4ad18d-0c49-4bc2-90d5-5bca8f00d6ae"
 FileIO = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
 GeometryBasics = "5c1252a2-5f33-56bf-86c9-59e7332b4326"
+IterativeSolvers = "42fd0dbc-a981-5370-80f2-aaf504508153"
 JSON = "682c06a0-de6a-54ab-a142-c8b1cf79cde6"
 LazyArrays = "5078a376-72f3-5289-bfd5-ec5146d43c02"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+LinearMaps = "7a12625a-238d-50fd-b39a-03d52299707e"
 Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
 MeshIO = "7269a6da-0436-5bbc-96c2-40638cbb6118"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"

diff --git a/test/Operators.jl b/test/Operators.jl
@@ -11,6 +11,8 @@ using Random
 using GeometryBasics: Point2, Point3
 using StaticArrays: SVector
 using Statistics: mean, var
+using IterativeSolvers
+using LinearMaps
 
 Point2D = Point2{Float64}
 Point3D = Point3{Float64}
@@ -358,6 +360,16 @@ function plot_dual0form(sd, f0)
   f
 end
 
+x = rand(1113)
+Δᵣ = Δ(rect)
+@test IterativeSolvers.cg(Δᵣ,x) ≈ Δᵣ \ x
+
+
+
+
+
+
+
 function euler_equation_test(X♯, sd)
   interior_tris = setdiff(triangles(sd), boundary_inds(Val{2}, sd))