use updateindex! instead of += for matrix update ; version bump

Project.toml

@@ -1,7 +1,7 @@
 name = "VoronoiFVM"
 uuid = "82b139dc-5afc-11e9-35da-9b9bdfd336f3"
 authors = ["Juergen Fuhrmann <[email protected]>"]
-version = "0.6.4"
+version = "0.6.5"
 
 [deps]
 BinDeps = "9e28174c-4ba2-5203-b857-d8d62c4213ee"
@@ -23,7 +23,7 @@ BinDeps = "^0.8.0"
 DiffResults = "^1.0.0"
 DocStringExtensions = "^0.8.0"
 ElasticArrays = "^1.0.0"
-ExtendableSparse = "^0.2.0"
+ExtendableSparse = "^0.2.6"
 ForwardDiff = "^0.10.0"
 IterativeSolvers = "^0.8.0"
 SparseDiffTools = "^1.0.0"

docs/src/changes.md

@@ -1,5 +1,24 @@
 # Changes
-## v0.6, Dec 15 2019
+## v0.6.5 Jan 25 2019
+- use updateindex! for matrix, depend on ExtendableSparse 0.2.6
+
+## v0.6.4 2020-01-20
+- Rearranged + commented boundary assembly loop
+- Reworked + renamed some examples
+- Document that unknowns doesn't initialize values.
+
+## v0.6.3 2019-12-21 
+remove xcolptrs call
+Update dependency on ExtendableSparse
+
+## v0.6.2 2019-12-20
+Updated dependency list (Triangulate ^0.4.0)
+
+## v0.6.1, 2019-12-17
+- return "plotted" for being able to  place colormap
+- require Triangulate >= 0.3.0
+
+## v0.6.0, Dec 15 2019
 - Removed Triangle submodule, depend on new Triangulate.jl Triangle wrapper
 - link to source code in examples
 - boundary_dirichlet! etc methods for setting boundary conditions

src/vfvm_solver.jl

@@ -14,7 +14,7 @@ const value=ForwardDiff.value
 # Add value to matrix if it is nonzero
 @inline function _addnz(matrix,i,j,v::Tv,fac) where Tv
     if v!=zero(Tv)
-        matrix[i,j]+=v*fac
+        updateindex!(matrix,+,v*fac,i,j)
     end
 end