bugfix_addition

[mongibellili]

to test
in sol interpolation, mliqss can end simulation before ft because closest nextstatetime can be ft+a and previous simt<ft

investigate the 1000 here!

sol1 = solInterpolated(sol, i, sol.ft / 1000.0)

[mongibellili]

document cycle detection mechanism access from outside (by user)
document inputstep logic
document eventstep logic

[mongibellili]

allow p[5] += 10.0 instead of p[5] =p[5]+ 10.0

--> solved after v1.0.4 (in qssProblemDiscreteHelper) by:
args[j].head in [:(=),:+=, :-=, :*=, :/=]

[mongibellili]

a, b, c, d, I = p leads to wrong dependencies because the parameters do not get changed in the equations and then dD would not get the dependency on "p"

--> solved by adding in main the following:

elseif rhs==:p && lhs isa Expr && lhs.head==:tuple #params: fill dict of param
                for k=1: length(lhs.args)
                    param[lhs.args[k]]=:(p[$k])
                end

[mongibellili]

check under quantizer, in function computenextEventTime if the commented code is needed
(oldsignValue[j,1] != sign(ZCFun[0])) # && abs(oldsignValue[j,2]) >1e-9*absQ

--> solved: the commented code is needed (the buck example)

[mongibellili]

bug: t inside event is of type Taylor:
either transform whole expression or change t to t[0]

--> solved: t became t[0] using this code in disreteHelper:

function changeEventExprToFirstValue(ex::Expr)
  newEx=postwalk(ex) do a  
      if (a==:t)|| (a isa Expr && a.head == :ref && a.args[1]==:q) 
          outerRef=Expr(:ref)
          push!(outerRef.args,a)
          push!(outerRef.args,:(0))
          a=outerRef
      end
      return a
  end
  newEx
end

[mongibellili]

function randNoise()
return randn()*10.0
end
randnoise=randNoise()
ERROR: LoadError: MethodError: no method matching copy(::Symbol)

[mongibellili]

1. allow hidden/interior assignments inside events:

   if y > 0  # this is the event condition
      if x>0    # this is just a check.
        Cd=2.0
      end
   end

1. Solved by assigns = find_assignments(posEvExp) inside qssProblemDiscreteHelper.jl

2. allow multiple assignments inside events:

   if y > 0
        Cd, ρ = Cd_air, ρ_air
   end

[mongibellili]

jac_mode= warning is not valid for explicit mehods

[mongibellili]

stateVarName = odeex.args[1].args[3]
discrParamName = length(odeex.args[1].args)==4 ? odeex.args[1].args[4] : :p

this gives an error if user enters f(u, p, t) and not f(du,u, p, t)