get solve() working again, add VoltageDerivativeConstraints
refactor solve, fix logging and printing