Improve the algorithm efficiency.

Different approaches :

• julia profiling
• (asynchronous) parallelization
• Linear modelling (especially considering to have or not the artificial constraint x=x_t^k used to get the multiplier)
• Initialization