The "Principal Axis" Optimizer described here works particuarly well on my problem. Now I do know some constraints about my problem and would like to implement them to make it faster. Even though it is an unconstrained optimizer it performs better without specified constraints than the other constrained optimizers with specified constraints. Mathematica forbids to specify constraints for this particular method. Is there a way around that?