def simulate(self, p=None, tindex=None, dt=1):
h = pd.Series(data=0, index=self.stress[0].index, name=self.name)
for i in self.stress:
self.npoints = self.stress.index.size
b = self.rfunc.block(p, self.r[i]) # nparam-1 depending on rfunc
h += fftconvolve(self.stress[i], b, 'full')[:self.npoints]
if tindex is not None:
h = h[tindex]
return h
评论列表
文章目录