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

status

parent 8bac75bd
No related branches found
No related tags found
2 merge requests!181actual status of monodomain,!179added possibility to use different time integration schemes for ion...
...@@ -130,7 +130,7 @@ def table8(): ...@@ -130,7 +130,7 @@ def table8():
l0Path=pathBV+'SI0004/vtu/AT_l0j0.vtu' l0Path=pathBV+'SI0004/vtu/AT_l0j0.vtu'
#plotTact: #plotTact:
plot=Plot(lL,tL,sdt,endTime,1,listEvaluationPoints) plot=Plot(lL,tL,sdt,endTime,1,listEvaluationPoints)
plot.MeanTactSpace(SI,l0Path) #plot.MeanTactSpace(SI,l0Path)
#tables tact #tables tact
L=4 L=4
J=4 J=4
......
...@@ -399,7 +399,7 @@ def RMSRefError(refDict,path,pathLevel0): ...@@ -399,7 +399,7 @@ def RMSRefError(refDict,path,pathLevel0):
def computeMeanTact(p,pathLevel0): def computeMeanTact(p,pathLevel0):
evalPoints,ateP=read.readVTU(pathLevel0) evalPoints,ateP=read.readVTU(pathLevel0)
points, at = read.readVTU(p) points, at = read.readVTU(p)
unactivatedPoints=0 print('nach readVTU')
if len(points)==0: if len(points)==0:
return 0.0 return 0.0
......
...@@ -316,6 +316,7 @@ class Plot: ...@@ -316,6 +316,7 @@ class Plot:
tactList=[] tactList=[]
jVector=[] jVector=[]
for j in self.tL: for j in self.tL:
print(l,j)
tact=comp.computeMeanTact(eP+'vtu/AT_l'+str(l)+'j'+str(j)+'.vtu',l0Path) tact=comp.computeMeanTact(eP+'vtu/AT_l'+str(l)+'j'+str(j)+'.vtu',l0Path)
if tact > 0.0 : if tact > 0.0 :
tactList.append(tact) tactList.append(tact)
......
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