The lagrangian package

[Tags:bsd3, library]

Numerically solve convex lagrange multiplier problems with conjugate gradient descent.

Convexity is key, otherwise the descent algorithm can return the wrong answer.

Convexity can be tested by assuring that the hessian of the lagrangian is positive definite over region the function is defined in.

I have provided test that the hessian is positive definite at a point, which is something, but not enough to ensure that the whole function is convex.

Be that as it may, if you know what the your lagrangian is convex you can use solve to find the minimum.

For example, find the maximum entropy with the constraint that the probabilities add up to one.

solve 0.00001 (negate . sum . map (x -> x * log x), [(sum, 1)]) 3

Gives the answer ([0.33, 0.33, 0.33], [-0.09])

The first elements of the result pair are the arguments for the objective function at the minimum. The second elements are the lagrange multipliers.

Properties

Versions 0.1.0.0, 0.2.0.0, 0.2.0.1, 0.2.0.2, 0.3.0.0, 0.3.0.1, 0.4.0.0, 0.4.0.1, 0.5.0.0, 0.6.0.0, 0.6.0.1
Dependencies ad (==3.4.*), base (==4.6.*), hmatrix (==0.14.*), nonlinear-optimization (==0.3.*), vector (==0.10.*) [details]
License BSD3
Author Jonathan Fischoff
Maintainer jonathangfischoff@gmail.com
Stability Unknown
Category Math
Home page http://github.com/jfischoff/lagrangian
Uploaded Wed Mar 6 05:56:57 UTC 2013 by JonathanFischoff
Distributions NixOS:0.6.0.1
Downloads 1882 total (52 in the last 30 days)
Votes
0 []
Status Docs uploaded by user
Build status unknown [no reports yet]

Modules

[Index]

Downloads

Maintainer's Corner

For package maintainers and hackage trustees