List changes here:

The Scythe_Optimize module changed so that the function sent to BFGS
(and all its support functions) requires 3 Matrix arguments: starting
parameter values, a response Matrix, and a covariate Matrix.


Ability to select all row and column elements with the following
usage:
A(_,int)
A(int,_)
To accomplish this a struct called all_elements was added to
Scythe_Double_Matrix.h, and two new overloaded Matrix () operator
definitions were added.
 

c() function(s) added. c() takes a user defined number (between 2 and
26) of doubles as arguments and concatenates them into a column vector.


sort() and sortc() functions added to sort all elements of a Matrix
and all columns of a Matrix respectively. An implementation of
shellsort based on that in Sedgewick, 1992. Algorithms in C++ is used
for both.


log factorial function (logfact()) added. Returns the natural log of
the factorial of the input argument.

Poisson mass function (dpois()) and the natural log of the Poisson 
mass function (logdpois()) added.

function to generate scalar pseudo-random Poisson deviates (rpois()) 
added.

PolySpline class added. (This is still very much under development)