add README example

SciML · Apr 5, 2020 · 2a89e1f · 2a89e1f
1 parent f29d70a
commit 2a89e1f
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 0 deletions.
diff --git a/test/examples.jl b/test/examples.jl
@@ -0,0 +1,16 @@
+fcalls = 0
+function f(dx,x)
+  global fcalls += 1
+  for i in 2:length(x)-1
+    dx[i] = x[i-1] - 2x[i] + x[i+1]
+  end
+  dx[1] = -2x[1] + x[2]
+  dx[end] = x[end-1] - 2x[end]
+  nothing
+end
+
+using SparsityDetection, SparseArrays
+input = rand(10)
+output = similar(input)
+sparsity_pattern = sparsity!(f,output,input)
+jac = Float64.(sparse(sparsity_pattern))
diff --git a/test/runtests.jl b/test/runtests.jl
@@ -5,4 +5,5 @@ include("common.jl")
 @testset "Exploration" begin include("ifsandbuts.jl") end
 @testset "Hessian sparsity" begin include("hessian.jl") end
 @testset "Paraboloid example" begin include("paraboloid.jl") end
+@safetestset "Simple examples" begin include("examples.jl") end
 @safetestset "Global PDE example" begin include("global_PDE.jl") end