Lines Matching refs:j
484 do j = 1, nflow
486 Atau(:,i,j) =Atau(:,i,j) + A1(:,i,k)*PTau(:,k,j)
495 do j = 1, nflow
497 A1tauA0(:,i,j) =
498 & Atau(:,i,1)*A0(:,1,j) +
499 & Atau(:,i,2)*A0(:,2,j) +
500 & Atau(:,i,3)*A0(:,3,j) +
501 & Atau(:,i,4)*A0(:,4,j) +
502 & Atau(:,i,5)*A0(:,5,j)
508 do j = 1, nflow
510 stiff(:,i,j) = stiff(:,i,j) + (
511 & Atau(:,i,1)*A1(:,1,j)
512 & + Atau(:,i,2)*A1(:,2,j)
513 & + Atau(:,i,3)*A1(:,3,j)
514 & + Atau(:,i,4)*A1(:,4,j)
515 & + Atau(:,i,5)*A1(:,5,j)
522 do j = 1, nflow
524 stiff(:,i,j+5) = stiff(:,i,j+5) + (
525 & Atau(:,i,1)*A2(:,1,j)
526 & + Atau(:,i,2)*A2(:,2,j)
527 & + Atau(:,i,3)*A2(:,3,j)
528 & + Atau(:,i,4)*A2(:,4,j)
529 & + Atau(:,i,5)*A2(:,5,j)
536 do j = 1, nflow
538 stiff(:,i,j+10) = stiff(:,i,j+10) + (
539 & Atau(:,i,1)*A3(:,1,j)
540 & + Atau(:,i,2)*A3(:,2,j)
541 & + Atau(:,i,3)*A3(:,3,j)
542 & + Atau(:,i,4)*A3(:,4,j)
543 & + Atau(:,i,5)*A3(:,5,j)
564 do j = 1, nflow
566 Atau(:,i,j) = Atau(:,i,j) + A2(:,i,k)*PTau(:,k,j)
575 do j = 1, nflow
577 A2tauA0(:,i,j) =
578 & Atau(:,i,1)*A0(:,1,j) +
579 & Atau(:,i,2)*A0(:,2,j) +
580 & Atau(:,i,3)*A0(:,3,j) +
581 & Atau(:,i,4)*A0(:,4,j) +
582 & Atau(:,i,5)*A0(:,5,j)
588 do j = 1, nflow
590 stiff(:,i+5,j) = stiff(:,i+5,j) + (
591 & Atau(:,i,1)*A1(:,1,j)
592 & + Atau(:,i,2)*A1(:,2,j)
593 & + Atau(:,i,3)*A1(:,3,j)
594 & + Atau(:,i,4)*A1(:,4,j)
595 & + Atau(:,i,5)*A1(:,5,j)
602 do j = 1, nflow
604 stiff(:,i+5,j+5) = stiff(:,i+5,j+5) + (
605 & Atau(:,i,1)*A2(:,1,j)
606 & + Atau(:,i,2)*A2(:,2,j)
607 & + Atau(:,i,3)*A2(:,3,j)
608 & + Atau(:,i,4)*A2(:,4,j)
609 & + Atau(:,i,5)*A2(:,5,j)
616 do j = 1, nflow
618 stiff(:,i+5,j+10) = stiff(:,i+5,j+10) + (
619 & Atau(:,i,1)*A3(:,1,j)
620 & + Atau(:,i,2)*A3(:,2,j)
621 & + Atau(:,i,3)*A3(:,3,j)
622 & + Atau(:,i,4)*A3(:,4,j)
623 & + Atau(:,i,5)*A3(:,5,j)
644 do j = 1, nflow
646 Atau(:,i,j) = Atau(:,i,j) + A3(:,i,k)*PTau(:,k,j)
655 do j = 1, nflow
657 A3tauA0(:,i,j) =
658 & Atau(:,i,1)*A0(:,1,j) +
659 & Atau(:,i,2)*A0(:,2,j) +
660 & Atau(:,i,3)*A0(:,3,j) +
661 & Atau(:,i,4)*A0(:,4,j) +
662 & Atau(:,i,5)*A0(:,5,j)
668 do j = 1, nflow
670 stiff(:,i+10,j) = stiff(:,i+10,j) + (
671 & Atau(:,i,1)*A1(:,1,j)
672 & + Atau(:,i,2)*A1(:,2,j)
673 & + Atau(:,i,3)*A1(:,3,j)
674 & + Atau(:,i,4)*A1(:,4,j)
675 & + Atau(:,i,5)*A1(:,5,j)
682 do j = 1, nflow
684 stiff(:,i+10,j+5) = stiff(:,i+10,j+5) + (
685 & Atau(:,i,1)*A2(:,1,j)
686 & + Atau(:,i,2)*A2(:,2,j)
687 & + Atau(:,i,3)*A2(:,3,j)
688 & + Atau(:,i,4)*A2(:,4,j)
689 & + Atau(:,i,5)*A2(:,5,j)
696 do j = 1, nflow
698 stiff(:,i+10,j+10) = stiff(:,i+10,j+10) + (
699 & Atau(:,i,1)*A3(:,1,j)
700 & + Atau(:,i,2)*A3(:,2,j)
701 & + Atau(:,i,3)*A3(:,3,j)
702 & + Atau(:,i,4)*A3(:,4,j)
703 & + Atau(:,i,5)*A3(:,5,j)
730 do j = 1, nshl
731 j0 = nflow * (j - 1)
735 fact = shape(:,j) * WdetJ * almi/gami/alfi*dtgl