now it turned out this is a bitch. I've spend quite some time trying to get the constraints working in an reliable way, but cant find a solution that works in more than 50% of the test cases. I checked my math with free available solvers (namely octave), turns out my functions are working. This means I have to dig deep and change basic solver algorithms, making the fork needed earlier than I though.
I have some ideas what I can do, starting with changing the linesearch to hold the wolf condition, but this will require some serious reading and therefore quite some time. I will keep you informed of my progress.
Greetings,
Stefan