[Developers] More informative Hessian info in ADMB-RE

[dave fournier]

The block diagonal Hessian

is in

 d3_array * block_diagonal_hessian;

in the structure usually pointed to by lapprox if it is allocated.

Note that there is no requirement that the size of each block be the
same in more complicated models.

So
 (*(lapprox->block_diagonal_hessian))(i)

is the i'th block

and

 (*(lapprox->block_diagonal_hessian))(i).indexmax()
  -(*(lapprox->block_diagonal_hessian))(i).indexmin()+1 is the size of
the i'th block


-- 
David A. Fournier
P.O. Box 2040,
Sidney, B.C. V8l 3S3
Canada
Phone/FAX 250-655-3364
http://otter-rsch.com