From fff290cd56973b5cd85615f3465a4acb4cefc258 Mon Sep 17 00:00:00 2001
From: Laura Pfeiffer <laura.pfeiffer@kit.edu>
Date: Sun, 14 Jul 2024 07:43:20 +0200
Subject: [PATCH] fixed evaluating timestep for splitting read
---
tools/evaluation/thesisBiventricle.py | 11 ++++++-----
tools/evaluation/utility/plotting.py | 2 +-
tools/evaluation/utility/reading.py | 2 +-
3 files changed, 8 insertions(+), 7 deletions(-)
diff --git a/tools/evaluation/thesisBiventricle.py b/tools/evaluation/thesisBiventricle.py
index 62ffb1f94..27833b79f 100644
--- a/tools/evaluation/thesisBiventricle.py
+++ b/tools/evaluation/thesisBiventricle.py
@@ -241,22 +241,23 @@ def WorkPrecisonDiagram(case):
output_file.write('\\end{figure}\n')
if __name__=="__main__":
language='d'
- plotSchwinger()
+
#AllValueTabs()
#plotPotential()
diffcase ='tact'
#DifferencePerPointTabs(diffcase)
diffcase='VDiffInt'
- # DifferencePerPointTabs(diffcase)
+ #DifferencePerPointTabs(diffcase)
#activationTimeExtrapolateTab()
#tactValuesAndErrorPerPoint()
- #plotTactErrorPerPoint()
- #ErrorExtraVInt()
+ plotTactErrorPerPoint()
+ ErrorExtraVInt()
#WriteDurationPerTimeStepTab()
WorkPrecisonDiagram('')
WorkPrecisonDiagram('Grob')
WorkPrecisonDiagram('tact')
WorkPrecisonDiagram('tactGrob')
#plotSchwinger()
- # plotExtraError()
+ plotExtraError()
+ plotSchwinger()
diff --git a/tools/evaluation/utility/plotting.py b/tools/evaluation/utility/plotting.py
index c0500205e..145a9e77c 100644
--- a/tools/evaluation/utility/plotting.py
+++ b/tools/evaluation/utility/plotting.py
@@ -597,7 +597,7 @@ class Plot:
plt.close()
- def tactErrorInTime(self,fL,nL,threshold,p,l,sP):
+ def tactErrorInTime(self,fL,nL,threshold,p,l,sP,exP):
tactRef=comp.getTextFullExtrapolate(self.lL[-1],self.tL[-1],1,p,exP,threshold)
for i in range(len(nL)):
tactError=[]
diff --git a/tools/evaluation/utility/reading.py b/tools/evaluation/utility/reading.py
index c017bbbcd..dac041c0c 100644
--- a/tools/evaluation/utility/reading.py
+++ b/tools/evaluation/utility/reading.py
@@ -648,7 +648,7 @@ def extractDataFromLog(path,fn,option):
if getModel(path,fn)=='Coupled':
elphytimeStep=getElpyhDeltaTime(path,fn)
else:
- (T,elphytimeStep)=getTimeDataFromLog(path+'log/'+fn+'.log')
+ (T,elphytimeStep)=getTimeDataFromLog(path+'log/log_'+fn)
#To add Iext for t=0
#steps[current_step][12] = '0.0'
for l in f:
--
GitLab