From f45da0a51410301377be997f5721b2321f9f15c9 Mon Sep 17 00:00:00 2001
From: welthulk <g1403.2019@udoschmitz.de>
Date: Sun, 14 Apr 2024 16:43:00 +0200
Subject: [PATCH] hotfix issue#48

---
 Project.toml          |  2 +-
 docs/src/changelog.md |  5 +++++
 src/Sparlectra.jl     |  4 ++--
 src/jacobian.jl       | 10 +++++-----
 4 files changed, 13 insertions(+), 8 deletions(-)

diff --git a/Project.toml b/Project.toml
index 5552e69..042b3de 100644
--- a/Project.toml
+++ b/Project.toml
@@ -2,7 +2,7 @@ name = "Sparlectra"
 uuid = "31ce9bba-fd9d-44a1-b005-f5f509afda38"
 authors = ["Udo Schmitz"]
 description = "load flow calculation using newton-raphson"
-version = "0.4.17"
+version = "0.4.18"
 
 [deps]
 BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
diff --git a/docs/src/changelog.md b/docs/src/changelog.md
index 0cc4fd5..3c492e1 100644
--- a/docs/src/changelog.md
+++ b/docs/src/changelog.md
@@ -1,4 +1,9 @@
 # Change Log
+## Version 0.4.18 (2024-04-14)
+
+### Bug Fixes
+ - closes issue#48 "printout jacobian runs into error"
+
 ## Version 0.4.17 (2024-04-14)
 ### New Features
  - added testcase for importing Matpower files
diff --git a/src/Sparlectra.jl b/src/Sparlectra.jl
index e328297..305a064 100644
--- a/src/Sparlectra.jl
+++ b/src/Sparlectra.jl
@@ -1,5 +1,5 @@
 """
-    Sparlectra 0.4.17
+    Sparlectra 0.4.18
 
 Sparlectra is a Julia package for the calculation of electrical networks. It is designed to be used in the context of power system analysis and optimization. 
 
@@ -27,7 +27,7 @@ using LinearAlgebra, SparseArrays, Printf, Logging
 
 # resource data types for working with Sparlectra
 const Wurzel3 = 1.7320508075688772
-const SparlectraVersion = VersionNumber("0.4.17")
+const SparlectraVersion = VersionNumber("0.4.18")
 abstract type AbstractBranch end
 
 export
diff --git a/src/jacobian.jl b/src/jacobian.jl
index b9a7c70..a380b47 100644
--- a/src/jacobian.jl
+++ b/src/jacobian.jl
@@ -676,16 +676,16 @@ function calcNewtonRaphson!(Y::AbstractMatrix{ComplexF64}, nodes::Vector{Node},
     end
 
     if verbose > 3
-      vm_pu = round(nodes[i]._vm_pu, digits = 3)
-      va_deg = round(nodes[i]._va_deg, digits = 3)
+      vm_pu = round(nodes[idx]._vm_pu, digits = 3)
+      va_deg = round(nodes[idx]._va_deg, digits = 3)
       if va_deg > 180
         va_deg -= 360
         va_deg = round(va_deg, digits = 3)
       end
-      if i == lastNode
-        println("(Bus=$(i), vm=$(vm_pu), va=$(va_deg))")
+      if idx == lastNode
+        println("(Bus=$(idx), vm=$(vm_pu), va=$(va_deg))")
       else
-        print("(Bus=$(i), vm=$(vm_pu), va=$(va_deg)), ")
+        print("(Bus=$(idx), vm=$(vm_pu), va=$(va_deg)), ")
       end
     end
   end