[Developers] More informative Hessian info in ADMB-RE

[Hans Julius Skaug]

Thanks,
 
I was planning to print the min and max block sizes in cases
where they are not equal. With your code snippet I can find
the min and max.
 
Hans

________________________________

Fra: developers-bounces at admb-project.org på vegne av dave fournier
Sendt: lø 15.08.2009 14:41
Til: developers at admb-project.org
Emne: Re: [Developers] More informative Hessian info in ADMB-RE





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 <http://otter-rsch.com/> 
_______________________________________________
Developers mailing list
Developers at admb-project.org
http://lists.admb-project.org/mailman/listinfo/developers