Lines Matching refs:Lx
64 double *Lx = L->x ; variable
87 y /= Lx [p] ;
90 X [j] = y / Lx [p] ;
94 X [Li [p]] -= Lx [p] * y ;
109 y [0] = X [j] / Lx [p] ;
110 y [1] = (X [j+1] - Lx [p+1] * y [0]) / Lx [q] ;
115 y [1] = X [j+1] - Lx [p+1] * y [0] ;
116 X [j ] = y [0] / Lx [p] ;
117 X [j+1] = y [1] / Lx [q] ;
120 y [1] = X [j+1] - Lx [p+1] * y [0] ;
125 X [Li [p]] -= Lx [p] * y [0] + Lx [q] * y [1] ;
141 y [0] = X [j] / Lx [p] ;
142 y [1] = (X [j+1] - Lx [p+1] * y [0]) / Lx [q] ;
143 y [2] = (X [j+2] - Lx [p+2] * y [0] - Lx [q+1] * y [1]) / Lx [r] ;
149 y [1] = X [j+1] - Lx [p+1] * y [0] ;
150 y [2] = X [j+2] - Lx [p+2] * y [0] - Lx [q+1] * y [1] ;
151 X [j ] = y [0] / Lx [p] ;
152 X [j+1] = y [1] / Lx [q] ;
153 X [j+2] = y [2] / Lx [r] ;
156 y [1] = X [j+1] - Lx [p+1] * y [0] ;
157 y [2] = X [j+2] - Lx [p+2] * y [0] - Lx [q+1] * y [1] ;
163 X [Li [p]] -= Lx [p] * y [0] + Lx [q] * y [1] + Lx [r] * y [2] ;
183 double *Lx = L->x ; variable
208 y [0] /= Lx [p] ;
209 y [1] /= Lx [p] ;
213 X [j][0] = y [0] / Lx [p] ;
214 X [j][1] = y [1] / Lx [p] ;
219 X [i][0] -= Lx [p] * y [0] ;
220 X [i][1] -= Lx [p] * y [1] ;
237 y [0][0] /= Lx [p] ;
238 y [0][1] /= Lx [p] ;
239 y [1][0] = (X [j+1][0] - Lx [p+1] * y [0][0]) / Lx [q] ;
240 y [1][1] = (X [j+1][1] - Lx [p+1] * y [0][1]) / Lx [q] ;
246 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
247 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
248 X [j ][0] = y [0][0] / Lx [p] ;
249 X [j ][1] = y [0][1] / Lx [p] ;
250 X [j+1][0] = y [1][0] / Lx [q] ;
251 X [j+1][1] = y [1][1] / Lx [q] ;
253 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
254 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
261 X [i][0] -= Lx [p] * y [0][0] + Lx [q] * y [1][0] ;
262 X [i][1] -= Lx [p] * y [0][1] + Lx [q] * y [1][1] ;
280 y [0][0] /= Lx [p] ;
281 y [0][1] /= Lx [p] ;
282 y [1][0] = (X [j+1][0] - Lx[p+1] * y[0][0]) / Lx [q] ;
283 y [1][1] = (X [j+1][1] - Lx[p+1] * y[0][1]) / Lx [q] ;
284 y [2][0] = (X [j+2][0] - Lx[p+2] * y[0][0] - Lx[q+1]*y[1][0])/Lx[r];
285 y [2][1] = (X [j+2][1] - Lx[p+2] * y[0][1] - Lx[q+1]*y[1][1])/Lx[r];
293 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
294 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
295 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
296 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
297 X [j ][0] = y [0][0] / Lx [p] ;
298 X [j ][1] = y [0][1] / Lx [p] ;
299 X [j+1][0] = y [1][0] / Lx [q] ;
300 X [j+1][1] = y [1][1] / Lx [q] ;
301 X [j+2][0] = y [2][0] / Lx [r] ;
302 X [j+2][1] = y [2][1] / Lx [r] ;
304 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
305 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
306 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
307 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
316 X[i][0] -= Lx[p] * y[0][0] + Lx[q] * y[1][0] + Lx[r] * y[2][0] ;
317 X[i][1] -= Lx[p] * y[0][1] + Lx[q] * y[1][1] + Lx[r] * y[2][1] ;
337 double *Lx = L->x ; variable
363 y [0] /= Lx [p] ;
364 y [1] /= Lx [p] ;
365 y [2] /= Lx [p] ;
370 X [j][0] = y [0] / Lx [p] ;
371 X [j][1] = y [1] / Lx [p] ;
372 X [j][2] = y [2] / Lx [p] ;
377 double lx = Lx [p] ;
398 y [0][0] /= Lx [p] ;
399 y [0][1] /= Lx [p] ;
400 y [0][2] /= Lx [p] ;
401 y [1][0] = (X [j+1][0] - Lx [p+1] * y [0][0]) / Lx [q] ;
402 y [1][1] = (X [j+1][1] - Lx [p+1] * y [0][1]) / Lx [q] ;
403 y [1][2] = (X [j+1][2] - Lx [p+1] * y [0][2]) / Lx [q] ;
411 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
412 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
413 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
414 X [j ][0] = y [0][0] / Lx [p] ;
415 X [j ][1] = y [0][1] / Lx [p] ;
416 X [j ][2] = y [0][2] / Lx [p] ;
417 X [j+1][0] = y [1][0] / Lx [q] ;
418 X [j+1][1] = y [1][1] / Lx [q] ;
419 X [j+1][2] = y [1][2] / Lx [q] ;
421 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
422 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
423 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
432 lx [0] = Lx [p] ;
433 lx [1] = Lx [q] ;
455 y [0][0] /= Lx [p] ;
456 y [0][1] /= Lx [p] ;
457 y [0][2] /= Lx [p] ;
458 y [1][0] = (X [j+1][0] - Lx[p+1] * y[0][0]) / Lx [q] ;
459 y [1][1] = (X [j+1][1] - Lx[p+1] * y[0][1]) / Lx [q] ;
460 y [1][2] = (X [j+1][2] - Lx[p+1] * y[0][2]) / Lx [q] ;
461 y [2][0] = (X [j+2][0] - Lx[p+2] * y[0][0] - Lx[q+1]*y[1][0])/Lx[r];
462 y [2][1] = (X [j+2][1] - Lx[p+2] * y[0][1] - Lx[q+1]*y[1][1])/Lx[r];
463 y [2][2] = (X [j+2][2] - Lx[p+2] * y[0][2] - Lx[q+1]*y[1][2])/Lx[r];
474 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
475 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
476 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
477 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
478 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
479 y [2][2] = X [j+2][2] - Lx [p+2] * y [0][2] - Lx [q+1] * y [1][2] ;
480 X [j ][0] = y [0][0] / Lx [p] ;
481 X [j ][1] = y [0][1] / Lx [p] ;
482 X [j ][2] = y [0][2] / Lx [p] ;
483 X [j+1][0] = y [1][0] / Lx [q] ;
484 X [j+1][1] = y [1][1] / Lx [q] ;
485 X [j+1][2] = y [1][2] / Lx [q] ;
486 X [j+2][0] = y [2][0] / Lx [r] ;
487 X [j+2][1] = y [2][1] / Lx [r] ;
488 X [j+2][2] = y [2][2] / Lx [r] ;
490 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
491 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
492 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
493 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
494 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
495 y [2][2] = X [j+2][2] - Lx [p+2] * y [0][2] - Lx [q+1] * y [1][2] ;
507 lx [0] = Lx [p] ;
508 lx [1] = Lx [q] ;
509 lx [2] = Lx [r] ;
532 double *Lx = L->x ; variable
559 y [0] /= Lx [p] ;
560 y [1] /= Lx [p] ;
561 y [2] /= Lx [p] ;
562 y [3] /= Lx [p] ;
568 X [j][0] = y [0] / Lx [p] ;
569 X [j][1] = y [1] / Lx [p] ;
570 X [j][2] = y [2] / Lx [p] ;
571 X [j][3] = y [3] / Lx [p] ;
576 double lx = Lx [p] ;
599 y [0][0] /= Lx [p] ;
600 y [0][1] /= Lx [p] ;
601 y [0][2] /= Lx [p] ;
602 y [0][3] /= Lx [p] ;
603 y [1][0] = (X [j+1][0] - Lx [p+1] * y [0][0]) / Lx [q] ;
604 y [1][1] = (X [j+1][1] - Lx [p+1] * y [0][1]) / Lx [q] ;
605 y [1][2] = (X [j+1][2] - Lx [p+1] * y [0][2]) / Lx [q] ;
606 y [1][3] = (X [j+1][3] - Lx [p+1] * y [0][3]) / Lx [q] ;
616 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
617 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
618 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
619 y [1][3] = X [j+1][3] - Lx [p+1] * y [0][3] ;
620 X [j ][0] = y [0][0] / Lx [p] ;
621 X [j ][1] = y [0][1] / Lx [p] ;
622 X [j ][2] = y [0][2] / Lx [p] ;
623 X [j ][3] = y [0][3] / Lx [p] ;
624 X [j+1][0] = y [1][0] / Lx [q] ;
625 X [j+1][1] = y [1][1] / Lx [q] ;
626 X [j+1][2] = y [1][2] / Lx [q] ;
627 X [j+1][3] = y [1][3] / Lx [q] ;
629 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
630 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
631 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
632 y [1][3] = X [j+1][3] - Lx [p+1] * y [0][3] ;
642 lx [0] = Lx [p] ;
643 lx [1] = Lx [q] ;
667 y [0][0] /= Lx [p] ;
668 y [0][1] /= Lx [p] ;
669 y [0][2] /= Lx [p] ;
670 y [0][3] /= Lx [p] ;
671 y [1][0] = (X [j+1][0] - Lx[p+1] * y[0][0]) / Lx [q] ;
672 y [1][1] = (X [j+1][1] - Lx[p+1] * y[0][1]) / Lx [q] ;
673 y [1][2] = (X [j+1][2] - Lx[p+1] * y[0][2]) / Lx [q] ;
674 y [1][3] = (X [j+1][3] - Lx[p+1] * y[0][3]) / Lx [q] ;
675 y [2][0] = (X [j+2][0] - Lx[p+2] * y[0][0] - Lx[q+1]*y[1][0])/Lx[r];
676 y [2][1] = (X [j+2][1] - Lx[p+2] * y[0][1] - Lx[q+1]*y[1][1])/Lx[r];
677 y [2][2] = (X [j+2][2] - Lx[p+2] * y[0][2] - Lx[q+1]*y[1][2])/Lx[r];
678 y [2][3] = (X [j+2][3] - Lx[p+2] * y[0][3] - Lx[q+1]*y[1][3])/Lx[r];
692 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
693 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
694 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
695 y [1][3] = X [j+1][3] - Lx [p+1] * y [0][3] ;
696 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
697 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
698 y [2][2] = X [j+2][2] - Lx [p+2] * y [0][2] - Lx [q+1] * y [1][2] ;
699 y [2][3] = X [j+2][3] - Lx [p+2] * y [0][3] - Lx [q+1] * y [1][3] ;
700 X [j ][0] = y [0][0] / Lx [p] ;
701 X [j ][1] = y [0][1] / Lx [p] ;
702 X [j ][2] = y [0][2] / Lx [p] ;
703 X [j ][3] = y [0][3] / Lx [p] ;
704 X [j+1][0] = y [1][0] / Lx [q] ;
705 X [j+1][1] = y [1][1] / Lx [q] ;
706 X [j+1][2] = y [1][2] / Lx [q] ;
707 X [j+1][3] = y [1][3] / Lx [q] ;
708 X [j+2][0] = y [2][0] / Lx [r] ;
709 X [j+2][1] = y [2][1] / Lx [r] ;
710 X [j+2][2] = y [2][2] / Lx [r] ;
711 X [j+2][3] = y [2][3] / Lx [r] ;
713 y [1][0] = X [j+1][0] - Lx [p+1] * y [0][0] ;
714 y [1][1] = X [j+1][1] - Lx [p+1] * y [0][1] ;
715 y [1][2] = X [j+1][2] - Lx [p+1] * y [0][2] ;
716 y [1][3] = X [j+1][3] - Lx [p+1] * y [0][3] ;
717 y [2][0] = X [j+2][0] - Lx [p+2] * y [0][0] - Lx [q+1] * y [1][0] ;
718 y [2][1] = X [j+2][1] - Lx [p+2] * y [0][1] - Lx [q+1] * y [1][1] ;
719 y [2][2] = X [j+2][2] - Lx [p+2] * y [0][2] - Lx [q+1] * y [1][2] ;
720 y [2][3] = X [j+2][3] - Lx [p+2] * y [0][3] - Lx [q+1] * y [1][3] ;
734 lx [0] = Lx [p] ;
735 lx [1] = Lx [q] ;
736 lx [2] = Lx [r] ;
768 double *Lx = L->x ; in LSOLVE() local
828 DIV_REAL (yx,yz,0, yx,yz,0, Lx,p) ; in LSOLVE()
832 DIV_REAL (Xx,Xz,j, yx,yz,0, Lx,p) ; in LSOLVE()
839 MULTSUB (Xx,Xz,i, Lx,Lz,p, yx,yz,0) ; in LSOLVE()