[ADMB Users] Matrix functions for ADMB-RE

[dave fournier]

Turns out that there is a stupid mismatch between the solve for 
dvariables and the
solve for df1b2variables.  Easiest workaround is to change sgn to a 
dvariable in the
tpl as in



dvariable ld;
dvariable sgn;
df1b2vector tmp=solve(varcov,diff,ld,sgn);

nLL+=.5*(nobs*log(2*M_PI)+ld+ diff*diff);


and then to add this section at the bottom of your tpl file.



GLOBALS_SECTION

   #include <admodel.h>
   #include <df1b2fun.h>
   #include <adrndeff.h>

   df1b2vector solve(const df1b2matrix& aa,const df1b2vector& z,
     const df1b2variable & ln_unsigned_det,double& sign);


   df1b2vector solve(const df1b2matrix& aa,const df1b2vector& z,
     const df1b2variable & ld,df1b2variable& sign)
   {
     double sgn;
     return solve(aa,z,ld,sgn);
   }