- Timestamp:
- 09/24/08 13:07:50 (16 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
bdm/math/libDC.cpp
r168 r170 157 157 double x = 0.0, sum; 158 158 int i,j; 159 vec s(v.length());160 vec S=L*v;161 159 162 160 for ( i=0; i<D.length(); i++ ) { //rows of L … … 164 162 for ( j=0; j<=i; j++ ){sum+=L( i,j )*v( j );} 165 163 x +=D( i )*sum*sum; 166 s(i)=sum;167 164 }; 168 165 return x;