def func2():
temp = sp.zeros(len(geteval(core, 'v')), 1)
v = geteval(core, 'v')
Mvv = geteval(core, 'Mvv')
Mvy = geteval(core, 'Mvy')
temp += matvecprod(core.I(''), sp.Matrix(v))
temp -= matvecprod(Mvv, sp.Matrix(core.f()))
temp -= matvecprod(Mvy, sp.Matrix(core.u))
return list(temp)
def func():
temp = sp.zeros(len(geteval(core, 'v' + n1)), 1)
v = geteval(core, 'v' + n1)
Mvvl = geteval(core, 'Mv' + n1 + 'vl')
Mvvnl = geteval(core, 'Mv' + n1 + 'vnl')
Mvy = geteval(core, 'Mv' + n1 + 'y')
temp += matvecprod(core.I(n1), sp.Matrix(v))
temp -= matvecprod(Mvvl, sp.Matrix(core.fl()))
temp -= matvecprod(Mvvnl, sp.Matrix(core.fnl()))
temp -= matvecprod(Mvy, sp.Matrix(core.u))
return list(temp)
def func2():
F = geteval(core, 'tempF')
vl = geteval(core, 'vl')
JacFl = jacobian(F, vl)
return list(
regularize_dims(sp.Matrix(F)) -
matvecprod(JacFl, sp.Matrix(vl)))
def func():
Fn1 = geteval(core, 'tempF' + n1)
vl = geteval(core, 'vl')
JacFn1l = geteval(
core,
'jacF' + n1 + 'l',
)
return list(
regularize_dims(sp.Matrix(Fn1)) -
matvecprod(JacFn1l, sp.Matrix(vl)))
def func():
Fa = geteval(core, 'tempF'+a)
vl = geteval(core, 'vl')
JacFal = geteval(core, 'jactempF'+a+'l',)
return list(regularize_dims(sp.Matrix(Fa)) -
matvecprod(JacFal, sp.Matrix(vl)))