add dt_fixed
active penalization working!! some slight difference to matlab solution. parallel.