From 2a131c07ea7d5cfd8b7d8ad23d64a074b622b06f Mon Sep 17 00:00:00 2001 From: Juergen Fuhrmann Date: Sat, 25 Jan 2020 23:34:34 +0100 Subject: [PATCH] use updateindex! instead of += for matrix update ; version bump --- Project.toml | 4 ++-- docs/src/changes.md | 21 ++++++++++++++++++++- src/vfvm_solver.jl | 2 +- 3 files changed, 23 insertions(+), 4 deletions(-) diff --git a/Project.toml b/Project.toml index 34ddcd848..b77da6e3d 100644 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "VoronoiFVM" uuid = "82b139dc-5afc-11e9-35da-9b9bdfd336f3" authors = ["Juergen Fuhrmann "] -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" diff --git a/docs/src/changes.md b/docs/src/changes.md index 34e7ff78d..8711875bc 100644 --- a/docs/src/changes.md +++ b/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 diff --git a/src/vfvm_solver.jl b/src/vfvm_solver.jl index 0352d47c4..b5a8ccec8 100644 --- a/src/vfvm_solver.jl +++ b/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