Skip to content
Snippets Groups Projects
Commit fff290cd authored by Laura Lindner's avatar Laura Lindner
Browse files

fixed evaluating timestep for splitting read

parent 60eee024
No related branches found
No related tags found
1 merge request!273fixed hopefully problem with computing displacement by depending mechassamble...
Pipeline #307052 failed
...@@ -241,22 +241,23 @@ def WorkPrecisonDiagram(case): ...@@ -241,22 +241,23 @@ def WorkPrecisonDiagram(case):
output_file.write('\\end{figure}\n') output_file.write('\\end{figure}\n')
if __name__=="__main__": if __name__=="__main__":
language='d' language='d'
plotSchwinger()
#AllValueTabs() #AllValueTabs()
#plotPotential() #plotPotential()
diffcase ='tact' diffcase ='tact'
#DifferencePerPointTabs(diffcase) #DifferencePerPointTabs(diffcase)
diffcase='VDiffInt' diffcase='VDiffInt'
# DifferencePerPointTabs(diffcase) #DifferencePerPointTabs(diffcase)
#activationTimeExtrapolateTab() #activationTimeExtrapolateTab()
#tactValuesAndErrorPerPoint() #tactValuesAndErrorPerPoint()
#plotTactErrorPerPoint() plotTactErrorPerPoint()
#ErrorExtraVInt() ErrorExtraVInt()
#WriteDurationPerTimeStepTab() #WriteDurationPerTimeStepTab()
WorkPrecisonDiagram('') WorkPrecisonDiagram('')
WorkPrecisonDiagram('Grob') WorkPrecisonDiagram('Grob')
WorkPrecisonDiagram('tact') WorkPrecisonDiagram('tact')
WorkPrecisonDiagram('tactGrob') WorkPrecisonDiagram('tactGrob')
#plotSchwinger() #plotSchwinger()
# plotExtraError() plotExtraError()
plotSchwinger()
...@@ -597,7 +597,7 @@ class Plot: ...@@ -597,7 +597,7 @@ class Plot:
plt.close() 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) tactRef=comp.getTextFullExtrapolate(self.lL[-1],self.tL[-1],1,p,exP,threshold)
for i in range(len(nL)): for i in range(len(nL)):
tactError=[] tactError=[]
......
...@@ -648,7 +648,7 @@ def extractDataFromLog(path,fn,option): ...@@ -648,7 +648,7 @@ def extractDataFromLog(path,fn,option):
if getModel(path,fn)=='Coupled': if getModel(path,fn)=='Coupled':
elphytimeStep=getElpyhDeltaTime(path,fn) elphytimeStep=getElpyhDeltaTime(path,fn)
else: else:
(T,elphytimeStep)=getTimeDataFromLog(path+'log/'+fn+'.log') (T,elphytimeStep)=getTimeDataFromLog(path+'log/log_'+fn)
#To add Iext for t=0 #To add Iext for t=0
#steps[current_step][12] = '0.0' #steps[current_step][12] = '0.0'
for l in f: for l in f:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment