Lines Matching refs:ay
165 INTEGER :: ax, ay, az, bx, by, bz, cda, cdax, cday, cdaz, coa, coamx, coamy, coamz, coapx, & local
279 DO ay = 2, la
280 az = la - ay
281 sc(coset(0, ay, az), 1) = rap(2)*sc(coset(0, ay - 1, az), 1) + &
282 f2*REAL(ay - 1, dp)*sc(coset(0, ay - 2, az), 1) - &
283 f2*kvec(2)*ss(coset(0, ay - 1, az), 1)
284 ss(coset(0, ay, az), 1) = rap(2)*ss(coset(0, ay - 1, az), 1) + &
285 f2*REAL(ay - 1, dp)*ss(coset(0, ay - 2, az), 1) + &
286 f2*kvec(2)*sc(coset(0, ay - 1, az), 1)
291 DO ay = 0, la - 1
292 az = la - 1 - ay
293 sc(coset(1, ay, az), 1) = rap(1)*sc(coset(0, ay, az), 1) - &
294 f2*kvec(1)*ss(coset(0, ay, az), 1)
295 ss(coset(1, ay, az), 1) = rap(1)*ss(coset(0, ay, az), 1) + &
296 f2*kvec(1)*sc(coset(0, ay, az), 1)
301 DO ay = 0, la - ax
302 az = la - ax - ay
303 sc(coset(ax, ay, az), 1) = rap(1)*sc(coset(ax - 1, ay, az), 1) + &
304 f3*sc(coset(ax - 2, ay, az), 1) - &
305 f2*kvec(1)*ss(coset(ax - 1, ay, az), 1)
306 ss(coset(ax, ay, az), 1) = rap(1)*ss(coset(ax - 1, ay, az), 1) + &
307 f3*ss(coset(ax - 2, ay, az), 1) + &
308 f2*kvec(1)*sc(coset(ax - 1, ay, az), 1)
339 DO ay = 0, la - ax
340 az = la - ax - ay
341 sc(coset(ax, ay, az), 2) = sc(coset(ax + 1, ay, az), 1) - &
342 rab(1)*sc(coset(ax, ay, az), 1)
343 sc(coset(ax, ay, az), 3) = sc(coset(ax, ay + 1, az), 1) - &
344 rab(2)*sc(coset(ax, ay, az), 1)
345 sc(coset(ax, ay, az), 4) = sc(coset(ax, ay, az + 1), 1) - &
346 rab(3)*sc(coset(ax, ay, az), 1)
347 ss(coset(ax, ay, az), 2) = ss(coset(ax + 1, ay, az), 1) - &
348 rab(1)*ss(coset(ax, ay, az), 1)
349 ss(coset(ax, ay, az), 3) = ss(coset(ax, ay + 1, az), 1) - &
350 rab(2)*ss(coset(ax, ay, az), 1)
351 ss(coset(ax, ay, az), 4) = ss(coset(ax, ay, az + 1), 1) - &
352 rab(3)*ss(coset(ax, ay, az), 1)
364 DO ay = 0, la_max - ax
365 fy = f2*REAL(ay, dp)
366 az = la_max - ax - ay
369 sc(coset(ax, ay, az), 2) = rbp(1)*sc(coset(ax, ay, az), 1) - &
370 f2*kvec(1)*ss(coset(ax, ay, az), 1)
371 ss(coset(ax, ay, az), 2) = rbp(1)*ss(coset(ax, ay, az), 1) + &
372 f2*kvec(1)*sc(coset(ax, ay, az), 1)
374 sc(coset(ax, ay, az), 2) = rbp(1)*sc(coset(ax, ay, az), 1) + &
375 fx*sc(coset(ax - 1, ay, az), 1) - &
376 f2*kvec(1)*ss(coset(ax, ay, az), 1)
377 ss(coset(ax, ay, az), 2) = rbp(1)*ss(coset(ax, ay, az), 1) + &
378 fx*ss(coset(ax - 1, ay, az), 1) + &
379 f2*kvec(1)*sc(coset(ax, ay, az), 1)
381 IF (ay == 0) THEN
382 sc(coset(ax, ay, az), 3) = rbp(2)*sc(coset(ax, ay, az), 1) - &
383 f2*kvec(2)*ss(coset(ax, ay, az), 1)
384 ss(coset(ax, ay, az), 3) = rbp(2)*ss(coset(ax, ay, az), 1) + &
385 f2*kvec(2)*sc(coset(ax, ay, az), 1)
387 sc(coset(ax, ay, az), 3) = rbp(2)*sc(coset(ax, ay, az), 1) + &
388 fy*sc(coset(ax, ay - 1, az), 1) - &
389 f2*kvec(2)*ss(coset(ax, ay, az), 1)
390 ss(coset(ax, ay, az), 3) = rbp(2)*ss(coset(ax, ay, az), 1) + &
391 fy*ss(coset(ax, ay - 1, az), 1) + &
392 f2*kvec(2)*sc(coset(ax, ay, az), 1)
395 sc(coset(ax, ay, az), 4) = rbp(3)*sc(coset(ax, ay, az), 1) - &
396 f2*kvec(3)*ss(coset(ax, ay, az), 1)
397 ss(coset(ax, ay, az), 4) = rbp(3)*ss(coset(ax, ay, az), 1) + &
398 f2*kvec(3)*sc(coset(ax, ay, az), 1)
400 sc(coset(ax, ay, az), 4) = rbp(3)*sc(coset(ax, ay, az), 1) + &
401 fz*sc(coset(ax, ay, az - 1), 1) - &
402 f2*kvec(3)*ss(coset(ax, ay, az), 1)
403 ss(coset(ax, ay, az), 4) = rbp(3)*ss(coset(ax, ay, az), 1) + &
404 fz*ss(coset(ax, ay, az - 1), 1) + &
405 f2*kvec(3)*sc(coset(ax, ay, az), 1)
426 DO ay = 0, la - ax
427 az = la - ax - ay
431 sc(coset(ax, ay, az), coset(0, 0, lb)) = &
432 sc(coset(ax, ay, az + 1), coset(0, 0, lb - 1)) - &
433 rab(3)*sc(coset(ax, ay, az), coset(0, 0, lb - 1))
434 ss(coset(ax, ay, az), coset(0, 0, lb)) = &
435 ss(coset(ax, ay, az + 1), coset(0, 0, lb - 1)) - &
436 rab(3)*ss(coset(ax, ay, az), coset(0, 0, lb - 1))
442 sc(coset(ax, ay, az), coset(0, by, bz)) = &
443 sc(coset(ax, ay + 1, az), coset(0, by - 1, bz)) - &
444 rab(2)*sc(coset(ax, ay, az), coset(0, by - 1, bz))
445 ss(coset(ax, ay, az), coset(0, by, bz)) = &
446 ss(coset(ax, ay + 1, az), coset(0, by - 1, bz)) - &
447 rab(2)*ss(coset(ax, ay, az), coset(0, by - 1, bz))
455 sc(coset(ax, ay, az), coset(bx, by, bz)) = &
456 sc(coset(ax + 1, ay, az), coset(bx - 1, by, bz)) - &
457 rab(1)*sc(coset(ax, ay, az), coset(bx - 1, by, bz))
458 ss(coset(ax, ay, az), coset(bx, by, bz)) = &
459 ss(coset(ax + 1, ay, az), coset(bx - 1, by, bz)) - &
460 rab(1)*ss(coset(ax, ay, az), coset(bx - 1, by, bz))
475 DO ay = 0, la_max - ax
476 fy = f2*REAL(ay, dp)
477 az = la_max - ax - ay
485 sc(coset(ax, ay, az), coset(0, 0, lb)) = &
486 rbp(3)*sc(coset(ax, ay, az), coset(0, 0, lb - 1)) + &
487 f3*sc(coset(ax, ay, az), coset(0, 0, lb - 2)) - &
488 f2*kvec(3)*ss(coset(ax, ay, az), coset(0, 0, lb - 1))
489 ss(coset(ax, ay, az), coset(0, 0, lb)) = &
490 rbp(3)*ss(coset(ax, ay, az), coset(0, 0, lb - 1)) + &
491 f3*ss(coset(ax, ay, az), coset(0, 0, lb - 2)) + &
492 f2*kvec(3)*sc(coset(ax, ay, az), coset(0, 0, lb - 1))
494 sc(coset(ax, ay, az), coset(0, 0, lb)) = &
495 rbp(3)*sc(coset(ax, ay, az), coset(0, 0, lb - 1)) + &
496 fz*sc(coset(ax, ay, az - 1), coset(0, 0, lb - 1)) + &
497 f3*sc(coset(ax, ay, az), coset(0, 0, lb - 2)) - &
498 f2*kvec(3)*ss(coset(ax, ay, az), coset(0, 0, lb - 1))
499 ss(coset(ax, ay, az), coset(0, 0, lb)) = &
500 rbp(3)*ss(coset(ax, ay, az), coset(0, 0, lb - 1)) + &
501 fz*ss(coset(ax, ay, az - 1), coset(0, 0, lb - 1)) + &
502 f3*ss(coset(ax, ay, az), coset(0, 0, lb - 2)) + &
503 f2*kvec(3)*sc(coset(ax, ay, az), coset(0, 0, lb - 1))
508 IF (ay == 0) THEN
510 sc(coset(ax, ay, az), coset(0, 1, bz)) = &
511 rbp(2)*sc(coset(ax, ay, az), coset(0, 0, bz)) - &
512 f2*kvec(2)*ss(coset(ax, ay, az), coset(0, 0, bz))
513 ss(coset(ax, ay, az), coset(0, 1, bz)) = &
514 rbp(2)*ss(coset(ax, ay, az), coset(0, 0, bz)) + &
515 f2*kvec(2)*sc(coset(ax, ay, az), coset(0, 0, bz))
519 sc(coset(ax, ay, az), coset(0, by, bz)) = &
520 rbp(2)*sc(coset(ax, ay, az), coset(0, by - 1, bz)) + &
521 f3*sc(coset(ax, ay, az), coset(0, by - 2, bz)) - &
522 f2*kvec(2)*ss(coset(ax, ay, az), coset(0, by - 1, bz))
523 ss(coset(ax, ay, az), coset(0, by, bz)) = &
524 rbp(2)*ss(coset(ax, ay, az), coset(0, by - 1, bz)) + &
525 f3*ss(coset(ax, ay, az), coset(0, by - 2, bz)) + &
526 f2*kvec(2)*sc(coset(ax, ay, az), coset(0, by - 1, bz))
530 sc(coset(ax, ay, az), coset(0, 1, bz)) = &
531 rbp(2)*sc(coset(ax, ay, az), coset(0, 0, bz)) + &
532 fy*sc(coset(ax, ay - 1, az), coset(0, 0, bz)) - &
533 f2*kvec(2)*ss(coset(ax, ay, az), coset(0, 0, bz))
534 ss(coset(ax, ay, az), coset(0, 1, bz)) = &
535 rbp(2)*ss(coset(ax, ay, az), coset(0, 0, bz)) + &
536 fy*ss(coset(ax, ay - 1, az), coset(0, 0, bz)) + &
537 f2*kvec(2)*sc(coset(ax, ay, az), coset(0, 0, bz))
541 sc(coset(ax, ay, az), coset(0, by, bz)) = &
542 rbp(2)*sc(coset(ax, ay, az), coset(0, by - 1, bz)) + &
543 fy*sc(coset(ax, ay - 1, az), coset(0, by - 1, bz)) + &
544 f3*sc(coset(ax, ay, az), coset(0, by - 2, bz)) - &
545 f2*kvec(2)*ss(coset(ax, ay, az), coset(0, by - 1, bz))
546 ss(coset(ax, ay, az), coset(0, by, bz)) = &
547 rbp(2)*ss(coset(ax, ay, az), coset(0, by - 1, bz)) + &
548 fy*ss(coset(ax, ay - 1, az), coset(0, by - 1, bz)) + &
549 f3*ss(coset(ax, ay, az), coset(0, by - 2, bz)) + &
550 f2*kvec(2)*sc(coset(ax, ay, az), coset(0, by - 1, bz))
559 sc(coset(ax, ay, az), coset(1, by, bz)) = &
560 rbp(1)*sc(coset(ax, ay, az), coset(0, by, bz)) - &
561 f2*kvec(1)*ss(coset(ax, ay, az), coset(0, by, bz))
562 ss(coset(ax, ay, az), coset(1, by, bz)) = &
563 rbp(1)*ss(coset(ax, ay, az), coset(0, by, bz)) + &
564 f2*kvec(1)*sc(coset(ax, ay, az), coset(0, by, bz))
570 sc(coset(ax, ay, az), coset(bx, by, bz)) = &
571 rbp(1)*sc(coset(ax, ay, az), &
573 f3*sc(coset(ax, ay, az), coset(bx - 2, by, bz)) - &
574 f2*kvec(1)*ss(coset(ax, ay, az), coset(bx - 1, by, bz))
575 ss(coset(ax, ay, az), coset(bx, by, bz)) = &
576 rbp(1)*ss(coset(ax, ay, az), &
578 f3*ss(coset(ax, ay, az), coset(bx - 2, by, bz)) + &
579 f2*kvec(1)*sc(coset(ax, ay, az), coset(bx - 1, by, bz))
585 sc(coset(ax, ay, az), coset(1, by, bz)) = &
586 rbp(1)*sc(coset(ax, ay, az), coset(0, by, bz)) + &
587 fx*sc(coset(ax - 1, ay, az), coset(0, by, bz)) - &
588 f2*kvec(1)*ss(coset(ax, ay, az), coset(0, by, bz))
589 ss(coset(ax, ay, az), coset(1, by, bz)) = &
590 rbp(1)*ss(coset(ax, ay, az), coset(0, by, bz)) + &
591 fx*ss(coset(ax - 1, ay, az), coset(0, by, bz)) + &
592 f2*kvec(1)*sc(coset(ax, ay, az), coset(0, by, bz))
598 sc(coset(ax, ay, az), coset(bx, by, bz)) = &
599 rbp(1)*sc(coset(ax, ay, az), &
601 fx*sc(coset(ax - 1, ay, az), coset(bx - 1, by, bz)) + &
602 f3*sc(coset(ax, ay, az), coset(bx - 2, by, bz)) - &
603 f2*kvec(1)*ss(coset(ax, ay, az), coset(bx - 1, by, bz))
604 ss(coset(ax, ay, az), coset(bx, by, bz)) = &
605 rbp(1)*ss(coset(ax, ay, az), &
607 fx*ss(coset(ax - 1, ay, az), coset(bx - 1, by, bz)) + &
608 f3*ss(coset(ax, ay, az), coset(bx - 2, by, bz)) + &
609 f2*kvec(1)*sc(coset(ax, ay, az), coset(bx - 1, by, bz))
728 DO ay = 0, la - ax
729 fay = REAL(ay, dp)
730 az = la - ax - ay
732 coa = coset(ax, ay, az)
733 coamx = coset(ax - 1, ay, az)
734 coamy = coset(ax, ay - 1, az)
735 coamz = coset(ax, ay, az - 1)
736 coapx = coset(ax + 1, ay, az)
737 coapy = coset(ax, ay + 1, az)
738 coapz = coset(ax, ay, az + 1)
810 INTEGER :: ax, ay, az, bx, by, bz, i, ipgf, j, & local
945 DO ay = 2, la
946 az = la - ay
947 s(coset(0, ay, az), 1, l) = rap(2)*s(coset(0, ay - 1, az), 1, l) + &
948 … f2*REAL(ay - 1, dp)*s(coset(0, ay - 2, az), 1, l)
949 IF (ly1 > 0) s(coset(0, ay, az), 1, l) = s(coset(0, ay, az), 1, l) + &
950 f2y*s(coset(0, ay - 1, az), 1, ly1)
955 DO ay = 0, la - 1
956 az = la - 1 - ay
957 s(coset(1, ay, az), 1, l) = rap(1)*s(coset(0, ay, az), 1, l)
958 IF (lx1 > 0) s(coset(1, ay, az), 1, l) = s(coset(1, ay, az), 1, l) + &
959 f2x*s(coset(0, ay, az), 1, lx1)
964 DO ay = 0, la - ax
965 az = la - ax - ay
966 s(coset(ax, ay, az), 1, l) = rap(1)*s(coset(ax - 1, ay, az), 1, l) + &
967 f3*s(coset(ax - 2, ay, az), 1, l)
968 IF (lx1 > 0) s(coset(ax, ay, az), 1, l) = s(coset(ax, ay, az), 1, l) + &
969 … f2x*s(coset(ax - 1, ay, az), 1, lx1)
999 DO ay = 0, la - ax
1000 az = la - ax - ay
1001 s(coset(ax, ay, az), 2, l) = s(coset(ax + 1, ay, az), 1, l) - &
1002 rab(1)*s(coset(ax, ay, az), 1, l)
1003 s(coset(ax, ay, az), 3, l) = s(coset(ax, ay + 1, az), 1, l) - &
1004 rab(2)*s(coset(ax, ay, az), 1, l)
1005 s(coset(ax, ay, az), 4, l) = s(coset(ax, ay, az + 1), 1, l) - &
1006 rab(3)*s(coset(ax, ay, az), 1, l)
1018 DO ay = 0, la_max - ax
1019 fy = f2*REAL(ay, dp)
1020 az = la_max - ax - ay
1023 s(coset(ax, ay, az), 2, l) = rbp(1)*s(coset(ax, ay, az), 1, l)
1025 s(coset(ax, ay, az), 2, l) = rbp(1)*s(coset(ax, ay, az), 1, l) + &
1026 fx*s(coset(ax - 1, ay, az), 1, l)
1028 IF (lx1 > 0) s(coset(ax, ay, az), 2, l) = s(coset(ax, ay, az), 2, l) + &
1029 … f2x*s(coset(ax, ay, az), 1, lx1)
1030 IF (ay == 0) THEN
1031 s(coset(ax, ay, az), 3, l) = rbp(2)*s(coset(ax, ay, az), 1, l)
1033 s(coset(ax, ay, az), 3, l) = rbp(2)*s(coset(ax, ay, az), 1, l) + &
1034 fy*s(coset(ax, ay - 1, az), 1, l)
1036 IF (ly1 > 0) s(coset(ax, ay, az), 3, l) = s(coset(ax, ay, az), 3, l) + &
1037 … f2y*s(coset(ax, ay, az), 1, ly1)
1039 s(coset(ax, ay, az), 4, l) = rbp(3)*s(coset(ax, ay, az), 1, l)
1041 s(coset(ax, ay, az), 4, l) = rbp(3)*s(coset(ax, ay, az), 1, l) + &
1042 fz*s(coset(ax, ay, az - 1), 1, l)
1044 IF (lz1 > 0) s(coset(ax, ay, az), 4, l) = s(coset(ax, ay, az), 4, l) + &
1045 … f2z*s(coset(ax, ay, az), 1, lz1)
1065 DO ay = 0, la - ax
1066 az = la - ax - ay
1070 s(coset(ax, ay, az), coset(0, 0, lb), l) = &
1071 s(coset(ax, ay, az + 1), coset(0, 0, lb - 1), l) - &
1072 rab(3)*s(coset(ax, ay, az), coset(0, 0, lb - 1), l)
1078 s(coset(ax, ay, az), coset(0, by, bz), l) = &
1079 s(coset(ax, ay + 1, az), coset(0, by - 1, bz), l) - &
1080 rab(2)*s(coset(ax, ay, az), coset(0, by - 1, bz), l)
1088 s(coset(ax, ay, az), coset(bx, by, bz), l) = &
1089 s(coset(ax + 1, ay, az), coset(bx - 1, by, bz), l) - &
1090 rab(1)*s(coset(ax, ay, az), coset(bx - 1, by, bz), l)
1105 DO ay = 0, la_max - ax
1106 fy = f2*REAL(ay, dp)
1107 az = la_max - ax - ay
1115 s(coset(ax, ay, az), coset(0, 0, lb), l) = &
1116 rbp(3)*s(coset(ax, ay, az), coset(0, 0, lb - 1), l) + &
1117 f3*s(coset(ax, ay, az), coset(0, 0, lb - 2), l)
1119 s(coset(ax, ay, az), coset(0, 0, lb), l) = &
1120 rbp(3)*s(coset(ax, ay, az), coset(0, 0, lb - 1), l) + &
1121 fz*s(coset(ax, ay, az - 1), coset(0, 0, lb - 1), l) + &
1122 f3*s(coset(ax, ay, az), coset(0, 0, lb - 2), l)
1124 IF (lz1 > 0) s(coset(ax, ay, az), coset(0, 0, lb), l) = &
1125 s(coset(ax, ay, az), coset(0, 0, lb), l) + &
1126 f2z*s(coset(ax, ay, az), coset(0, 0, lb - 1), lz1)
1130 IF (ay == 0) THEN
1132 s(coset(ax, ay, az), coset(0, 1, bz), l) = &
1133 rbp(2)*s(coset(ax, ay, az), coset(0, 0, bz), l)
1134 IF (ly1 > 0) s(coset(ax, ay, az), coset(0, 1, bz), l) = &
1135 s(coset(ax, ay, az), coset(0, 1, bz), l) + &
1136 f2y*s(coset(ax, ay, az), coset(0, 0, bz), ly1)
1140 s(coset(ax, ay, az), coset(0, by, bz), l) = &
1141 rbp(2)*s(coset(ax, ay, az), coset(0, by - 1, bz), l) + &
1142 f3*s(coset(ax, ay, az), coset(0, by - 2, bz), l)
1143 IF (ly1 > 0) s(coset(ax, ay, az), coset(0, by, bz), l) = &
1144 s(coset(ax, ay, az), coset(0, by, bz), l) + &
1145 f2y*s(coset(ax, ay, az), coset(0, by - 1, bz), ly1)
1149 s(coset(ax, ay, az), coset(0, 1, bz), l) = &
1150 rbp(2)*s(coset(ax, ay, az), coset(0, 0, bz), l) + &
1151 fy*s(coset(ax, ay - 1, az), coset(0, 0, bz), l)
1152 IF (ly1 > 0) s(coset(ax, ay, az), coset(0, 1, bz), l) = &
1153 s(coset(ax, ay, az), coset(0, 1, bz), l) + &
1154 f2y*s(coset(ax, ay, az), coset(0, 0, bz), ly1)
1158 s(coset(ax, ay, az), coset(0, by, bz), l) = &
1159 rbp(2)*s(coset(ax, ay, az), coset(0, by - 1, bz), l) + &
1160 fy*s(coset(ax, ay - 1, az), coset(0, by - 1, bz), l) + &
1161 f3*s(coset(ax, ay, az), coset(0, by - 2, bz), l)
1162 IF (ly1 > 0) s(coset(ax, ay, az), coset(0, by, bz), l) = &
1163 s(coset(ax, ay, az), coset(0, by, bz), l) + &
1164 f2y*s(coset(ax, ay, az), coset(0, by - 1, bz), ly1)
1173 s(coset(ax, ay, az), coset(1, by, bz), l) = &
1174 rbp(1)*s(coset(ax, ay, az), coset(0, by, bz), l)
1175 IF (lx1 > 0) s(coset(ax, ay, az), coset(1, by, bz), l) = &
1176 s(coset(ax, ay, az), coset(1, by, bz), l) + &
1177 f2x*s(coset(ax, ay, az), coset(0, by, bz), lx1)
1183 s(coset(ax, ay, az), coset(bx, by, bz), l) = &
1184 rbp(1)*s(coset(ax, ay, az), coset(bx - 1, by, bz), l) + &
1185 f3*s(coset(ax, ay, az), coset(bx - 2, by, bz), l)
1186 IF (lx1 > 0) s(coset(ax, ay, az), coset(bx, by, bz), l) = &
1187 s(coset(ax, ay, az), coset(bx, by, bz), l) + &
1188 f2x*s(coset(ax, ay, az), coset(bx - 1, by, bz), lx1)
1194 s(coset(ax, ay, az), coset(1, by, bz), l) = &
1195 rbp(1)*s(coset(ax, ay, az), coset(0, by, bz), l) + &
1196 fx*s(coset(ax - 1, ay, az), coset(0, by, bz), l)
1197 IF (lx1 > 0) s(coset(ax, ay, az), coset(1, by, bz), l) = &
1198 s(coset(ax, ay, az), coset(1, by, bz), l) + &
1199 f2x*s(coset(ax, ay, az), coset(0, by, bz), lx1)
1205 s(coset(ax, ay, az), coset(bx, by, bz), l) = &
1206 rbp(1)*s(coset(ax, ay, az), coset(bx - 1, by, bz), l) + &
1207 fx*s(coset(ax - 1, ay, az), coset(bx - 1, by, bz), l) + &
1208 f3*s(coset(ax, ay, az), coset(bx - 2, by, bz), l)
1209 IF (lx1 > 0) s(coset(ax, ay, az), coset(bx, by, bz), l) = &
1210 s(coset(ax, ay, az), coset(bx, by, bz), l) + &
1211 f2x*s(coset(ax, ay, az), coset(bx - 1, by, bz), lx1)