[Developers] (extract_) diagonal

Arni Magnusson arnima at hafro.is
Mon May 31 08:08:37 PDT 2010


The current manuals for ADMB (p. 15-9) and AUTODIF (p. 4-11) describe a 
function 'extract_diagonal', but it seems to be implemented as 'diagonal'. 
To verify this, I used a slightly modified version of simple.tpl:

DATA_SECTION
   init_int n
   init_vector X(1,n)
   init_vector Y(1,n)
PARAMETER_SECTION
   init_number b0
   init_number b1
   vector Yfit(1,n)
   matrix foo(1,3,1,3)
   objective_function_value RSS
PRELIMINARY_CALCS_SECTION
   foo(1,1)=11; foo(1,2)=12; foo(1,3)=13;
   foo(2,1)=21; foo(2,2)=22; foo(2,3)=23;
   foo(3,1)=31; foo(3,2)=32; foo(3,3)=33;
PROCEDURE_SECTION
   Yfit = b0 + b1*X;
   RSS = norm2(Y-Yfit);
REPORT_SECTION
   report << diagonal(foo) << endl;

If you change 'diagonal' to 'extract_diagonal' the model does not compile. 
The function name should either be changed in the manuals, or in the 
source code.

Arni


More information about the Developers mailing list