From 07d7e2b142f310ccdf1b41bd66b748897e9161ee Mon Sep 17 00:00:00 2001
From: Sam Yates <halfflat@gmail.com>
Date: Wed, 26 Oct 2016 15:01:32 +0200
Subject: [PATCH] Use correct routine in numeric_rallpack1.jl x0.3
---
validation/ref/numeric/numeric_rallpack1.jl | 9 +++++++--
1 file changed, 7 insertions(+), 2 deletions(-)
diff --git a/validation/ref/numeric/numeric_rallpack1.jl b/validation/ref/numeric/numeric_rallpack1.jl
index d0f5997e..ea9755aa 100644
--- a/validation/ref/numeric/numeric_rallpack1.jl
+++ b/validation/ref/numeric/numeric_rallpack1.jl
@@ -6,8 +6,12 @@ using JSON
using SIUnits.ShortUnits
using PassiveCable
+# This should run the same effective model
+# as rallpack1, but with differing
+# electrical parameters (see below).
+
function run_cable(x_prop, ts)
- # physical properties
+ # Physical properties:
# f is a fudge factor. rM needs to be the same
# the same as in nestmc, where we cannot yet set
@@ -40,6 +44,7 @@ end
function run_rallpack1(x_prop, ts)
return [rallpack1(0.001*x_prop, t/s)*V for t in ts]
end
+
# Generate traces at x=0, x=0.3L, x=L
ts = collect(0s: 0.025ms: 250ms)
@@ -51,7 +56,7 @@ trace = Dict(
:data => Dict(
:time => map(t->t/ms, ts),
symbol("cable.x0.0") => map(v->v/mV, run_cable(0, ts)),
- symbol("cable.x0.3") => map(v->v/mV, run_rallpack1(0.3, ts)),
+ symbol("cable.x0.3") => map(v->v/mV, run_cable(0.3, ts)),
symbol("cable.x1.0") => map(v->v/mV, run_cable(1.0, ts))
)
)
--
GitLab