bugfixes: active penalization working, comparison with fortran okay.