Home | History | Annotate | Download | only in db_vlvm

Lines Matching refs:ss

179     double ss[6];
182 ss[0]=db_SquareSum7(A);
183 ss[1]=db_SquareSum7(A+7);
184 ss[2]=db_SquareSum7(A+14);
185 ss[3]=db_SquareSum7(A+21);
186 ss[4]=db_SquareSum7(A+28);
187 ss[5]=db_SquareSum7(A+35);
189 ss[1]-=db_OrthogonalizePair7(A+7 ,A,ss[0]);
190 ss[2]-=db_OrthogonalizePair7(A+14,A,ss[0]);
191 ss[3]-=db_OrthogonalizePair7(A+21,A,ss[0]);
192 ss[4]-=db_OrthogonalizePair7(A+28,A,ss[0]);
193 ss[5]-=db_OrthogonalizePair7(A+35,A,ss[0]);
195 /*Pivot on largest ss (could also be done on ss/(original_ss))*/
196 i=db_MaxIndex5(ss+1);
197 db_OrthogonalizationSwap7(A+7,i,ss+1);
199 ss[2]-=db_OrthogonalizePair7(A+14,A+7,ss[1]);
200 ss[3]-=db_OrthogonalizePair7(A+21,A+7,ss[1]);
201 ss[4]-=db_OrthogonalizePair7(A+28,A+7,ss[1]);
202 ss[5]-=db_OrthogonalizePair7(A+35,A+7,ss[1]);
204 i=db_MaxIndex4(ss+2);
205 db_OrthogonalizationSwap7(A+14,i,ss+2);
207 ss[3]-=db_OrthogonalizePair7(A+21,A+14,ss[2]);
208 ss[4]-=db_OrthogonalizePair7(A+28,A+14,ss[2]);
209 ss[5]-=db_OrthogonalizePair7(A+35,A+14,ss[2]);
211 i=db_MaxIndex3(ss+3);
212 db_OrthogonalizationSwap7(A+21,i,ss+3);
214 ss[4]-=db_OrthogonalizePair7(A+28,A+21,ss[3]);
215 ss[5]-=db_OrthogonalizePair7(A+35,A+21,ss[3]);
217 i=db_MaxIndex2(ss+4);
218 db_OrthogonalizationSwap7(A+28,i,ss+4);
220 ss[5]-=db_OrthogonalizePair7(A+35,A+28,ss[4]);
224 db_MultiplyScalar7(A ,db_SafeSqrtReciprocal(ss[0]));
225 db_MultiplyScalar7(A+7 ,db_SafeSqrtReciprocal(ss[1]));
226 db_MultiplyScalar7(A+14,db_SafeSqrtReciprocal(ss[2]));
227 db_MultiplyScalar7(A+21,db_SafeSqrtReciprocal(ss[3]));
228 db_MultiplyScalar7(A+28,db_SafeSqrtReciprocal(ss[4]));
229 db_MultiplyScalar7(A+35,db_SafeSqrtReciprocal(ss[5]));
236 double ss[8];
239 ss[0]=db_SquareSum9(A);
240 ss[1]=db_SquareSum9(A+9);
241 ss[2]=db_SquareSum9(A+18);
242 ss[3]=db_SquareSum9(A+27);
243 ss[4]=db_SquareSum9(A+36);
244 ss[5]=db_SquareSum9(A+45);
245 ss[6]=db_SquareSum9(A+54);
246 ss[7]=db_SquareSum9(A+63);
248 ss[1]-=db_OrthogonalizePair9(A+9 ,A,ss[0]);
249 ss[2]-=db_OrthogonalizePair9(A+18,A,ss[0]);
250 ss[3]-=db_OrthogonalizePair9(A+27,A,ss[0]);
251 ss[4]-=db_OrthogonalizePair9(A+36,A,ss[0]);
252 ss[5]-=db_OrthogonalizePair9(A+45,A,ss[0]);
253 ss[6]-=db_OrthogonalizePair9(A+54,A,ss[0]);
254 ss[7]-=db_OrthogonalizePair9(A+63,A,ss[0]);
256 /*Pivot on largest ss (could also be done on ss/(original_ss))*/
257 i=db_MaxIndex7(ss+1);
258 db_OrthogonalizationSwap9(A+9,i,ss+1);
260 ss[2]-=db_OrthogonalizePair9(A+18,A+9,ss[1]);
261 ss[3]-=db_OrthogonalizePair9(A+27,A+9,ss[1]);
262 ss[4]-=db_OrthogonalizePair9(A+36,A+9,ss[1]);
263 ss[5]-=db_OrthogonalizePair9(A+45,A+9,ss[1]);
264 ss[6]-=db_OrthogonalizePair9(A+54,A+9,ss[1]);
265 ss[7]-=db_OrthogonalizePair9(A+63,A+9,ss[1]);
267 i=db_MaxIndex6(ss+2);
268 db_OrthogonalizationSwap9(A+18,i,ss+2);
270 ss[3]-=db_OrthogonalizePair9(A+27,A+18,ss[2]);
271 ss[4]-=db_OrthogonalizePair9(A+36,A+18,ss[2]);
272 ss[5]-=db_OrthogonalizePair9(A+45,A+18,ss[2]);
273 ss[6]-=db_OrthogonalizePair9(A+54,A+18,ss[2]);
274 ss[7]-=db_OrthogonalizePair9(A+63,A+18,ss[2]);
276 i=db_MaxIndex5(ss+3);
277 db_OrthogonalizationSwap9(A+27,i,ss+3);
279 ss[4]-=db_OrthogonalizePair9(A+36,A+27,ss[3]);
280 ss[5]-=db_OrthogonalizePair9(A+45,A+27,ss[3]);
281 ss[6]-=db_OrthogonalizePair9(A+54,A+27,ss[3]);
282 ss[7]-=db_OrthogonalizePair9(A+63,A+27,ss[3]);
284 i=db_MaxIndex4(ss+4);
285 db_OrthogonalizationSwap9(A+36,i,ss+4);
287 ss[5]-=db_OrthogonalizePair9(A+45,A+36,ss[4]);
288 ss[6]-=db_OrthogonalizePair9(A+54,A+36,ss[4]);
289 ss[7]-=db_OrthogonalizePair9(A+63,A+36,ss[4]);
291 i=db_MaxIndex3(ss+5);
292 db_OrthogonalizationSwap9(A+45,i,ss+5);
294 ss[6]-=db_OrthogonalizePair9(A+54,A+45,ss[5]);
295 ss[7]-=db_OrthogonalizePair9(A+63,A+45,ss[5]);
297 i=db_MaxIndex2(ss+6);
298 db_OrthogonalizationSwap9(A+54,i,ss+6);
300 ss[7]-=db_OrthogonalizePair9(A+63,A+54,ss[6]);
304 db_MultiplyScalar9(A ,db_SafeSqrtReciprocal(ss[0]));
305 db_MultiplyScalar9(A+9 ,db_SafeSqrtReciprocal(ss[1]));
306 db_MultiplyScalar9(A+18,db_SafeSqrtReciprocal(ss[2]));
307 db_MultiplyScalar9(A+27,db_SafeSqrtReciprocal(ss[3]));
308 db_MultiplyScalar9(A+36,db_SafeSqrtReciprocal(ss[4]));
309 db_MultiplyScalar9(A+45,db_SafeSqrtReciprocal(ss[5]));
310 db_MultiplyScalar9(A+54,db_SafeSqrtReciprocal(ss[6]));
311 db_MultiplyScalar9(A+63,db_SafeSqrtReciprocal(ss[7]));