Implementing Sequential Quadratic Programming about spec HOT 1 OPEN

I am happy to announce I have a version of SQP that at least "seems to work".

I have removed all the Lagrangian multipliers in dMA, and moved them to form constraints of the problem.

The test case is InputFiles/Verification/helicity/Taylor/helicity20_SQP.sp

In the input file, I put LBeltrami=3. This means that first the SQP is called, then the result is used to initialize a Newton's method. The screen will output two lines of numbers, with the first in each one being the energy and second one being the helicity (constraint to 20). The first line is the result of SQP and second is the Newton's method. One can verify that in helicity20_SQP.sp.end, mu is indeed 3.13, the result of the Taylor relaxation noted earlier in #59 .

The results are

SQP: 59.420569927889716 19.999995061296641
Newton: 59.420585419870079 19.999999999999993

which seems to match

from spec.