我正在处理DAEs问题，出于准确性考虑，我打算使用具有3个搭配点的离散化方案（例如，我正在使用10个有限元）。然后，我只想使用有限元范围内的解来计算给定的度量，因为我正在考虑只有在那儿才有用于比较的度量。
考虑所有点的测量是
def _FIMmatrix(m, i, j):
return m.gFIM[i,j] == sum(sum(m.sen[k,1,i,t]*m.sen[k,1,j,t] for t in m.t) \
for k in m.e)/len(m.t)/len(m.e)
m.FIMmatrix = Constraint(m.p, m.p, rule=_FIMmatrix)
欢迎任何提示。我认为在一种策略中，本地计数器也应等于mod（k，m）eq 1，其中mod是模数，m是配置点的数量，但在Pyomo中找不到模数运算符。
Vasco Luna