Revision 1486
trunk/src/linad99/fvar_m20.cpp (revision 1486)  

25  25 
\n\n The implementation of this algorithm was inspired by 
26  26 
"Numerical Recipes in C", 2nd edition, 
27  27 
Press, Teukolsky, Vetterling, Flannery, chapter 2 
28  
29 
Edited by Steve Martell on Dec 30, 2013 to address a bug 

30 
where the lower bound of the square dvar_matrix was greater than 1. 

28  31 
*/ 
29  32 
dvariable det(const dvar_matrix& aa) 
30  33 
{ 
...  ...  
129  132 
} 
130  133 
} 
131  134 
double det=d; 
132 
part_prod(1)=d*bb(1,1); 

135  
136 
// SM Bug 129, issue appears to be at this line 

137 
// part_prod is declared above as dvector(lb,ub) 

138 
// cout<<"Bug 129 ="<<part_prod(lb)<<endl; 

139 
// part_prod(1)=d*bb(1,1); // replaced this line with: 

140 
part_prod(lb) = d*bb(lb,lb); 

141 
// cout<<"Ok got this far; det = "<<det<<endl; 

133  142 
for (j=lb+1;j<=ub;j++) 
134  143 
{ 
135  144 
part_prod(j)=part_prod(j1)*bb(j,j); 
trunk/contrib/statslib/statsLib.h (revision 1486)  

92  92 
dvar_vector plogis( const dvector& x, const prevariable& location, const prevariable& scale ); 
93  93 
dvar_vector plogis( const dvar_vector& x, const prevariable& location, const prevariable& scale ); 
94  94 
// template <typename T1, typename T2> 
95 
// T1 plogis( const T1& x, const T2& location, const T2& scale ); 

95 
// T1 plogis( const T1& x, const T2& locatiogn, const T2& scale );


96  96  
97  97 
// Exponential logistic distribution 
98  98 
dvar_vector eplogis(const dvar_vector& x, const dvariable& alpha, const dvariable& beta, const dvariable& gamma); 
Also available in: Unified diff