rcce/bt/solve_subs.c

213       double pivot, coeff;
219       pivot = 1.00e0/lhs(1,1);
220       lhs(1,2) = lhs(1,2)*pivot;
221       lhs(1,3) = lhs(1,3)*pivot;
222       lhs(1,4) = lhs(1,4)*pivot;
223       lhs(1,5) = lhs(1,5)*pivot;
224       c(1,1) = c(1,1)*pivot;
225       c(1,2) = c(1,2)*pivot;
226       c(1,3) = c(1,3)*pivot;
227       c(1,4) = c(1,4)*pivot;
228       c(1,5) = c(1,5)*pivot;
229       r(1)   = r(1)  *pivot;
280       pivot = 1.00e0/lhs(2,2);
281       lhs(2,3) = lhs(2,3)*pivot;
282       lhs(2,4) = lhs(2,4)*pivot;
283       lhs(2,5) = lhs(2,5)*pivot;
284       c(2,1) = c(2,1)*pivot;
285       c(2,2) = c(2,2)*pivot;
286       c(2,3) = c(2,3)*pivot;
287       c(2,4) = c(2,4)*pivot;
288       c(2,5) = c(2,5)*pivot;
289       r(2)   = r(2)  *pivot;
336       pivot = 1.00e0/lhs(3,3);
337       lhs(3,4) = lhs(3,4)*pivot;
338       lhs(3,5) = lhs(3,5)*pivot;
339       c(3,1) = c(3,1)*pivot;
340       c(3,2) = c(3,2)*pivot;
341       c(3,3) = c(3,3)*pivot;
342       c(3,4) = c(3,4)*pivot;
343       c(3,5) = c(3,5)*pivot;
344       r(3)   = r(3)  *pivot;
387       pivot = 1.00e0/lhs(4,4);
388       lhs(4,5) = lhs(4,5)*pivot;
389       c(4,1) = c(4,1)*pivot;
390       c(4,2) = c(4,2)*pivot;
391       c(4,3) = c(4,3)*pivot;
392       c(4,4) = c(4,4)*pivot;
393       c(4,5) = c(4,5)*pivot;
394       r(4)   = r(4)  *pivot;
433       pivot = 1.00e0/lhs(5,5);
434       c(5,1) = c(5,1)*pivot;
435       c(5,2) = c(5,2)*pivot;
436       c(5,3) = c(5,3)*pivot;
437       c(5,4) = c(5,4)*pivot;
438       c(5,5) = c(5,5)*pivot;
439       r(5)   = r(5)  *pivot;
491       double pivot, coeff;
498       pivot = 1.00e0/lhs(1,1);
499       lhs(1,2) = lhs(1,2)*pivot;
500       lhs(1,3) = lhs(1,3)*pivot;
501       lhs(1,4) = lhs(1,4)*pivot;
502       lhs(1,5) = lhs(1,5)*pivot;
503       r(1)   = r(1)  *pivot;
534       pivot = 1.00e0/lhs(2,2);
535       lhs(2,3) = lhs(2,3)*pivot;
536       lhs(2,4) = lhs(2,4)*pivot;
537       lhs(2,5) = lhs(2,5)*pivot;
538       r(2)   = r(2)  *pivot;
565       pivot = 1.00e0/lhs(3,3);
566       lhs(3,4) = lhs(3,4)*pivot;
567       lhs(3,5) = lhs(3,5)*pivot;
568       r(3)   = r(3)  *pivot;
591       pivot = 1.00e0/lhs(4,4);
592       lhs(4,5) = lhs(4,5)*pivot;
593       r(4)   = r(4)  *pivot;
612       pivot = 1.00e0/lhs(5,5);
613       r(5)   = r(5)  *pivot;