Lines Matching refs:EdgeBasis
2491 EdgeBasis, RotBasis, USolver ) RESULT(stat)
2509 …REAL(KIND=dp), OPTIONAL :: EdgeBasis(:,:) !< If present, the values of H(curl)-conforming bas… local
2532 IF(PRESENT(EdgeBasis)) THEN
2550 EdgeBasis=EdgeBasis,RotBasis=RotBasis,dBasisdx=dBasisdx,&
2555 CALL GetEdgeBasis(Element,EdgeBasis,RotBasis,Basis,dBasisdx)
5426 Basis, EdgeBasis, RotBasis, dBasisdx, SecondFamily, BasisDegree, &
5441 …REAL(KIND=dp) :: EdgeBasis(:,:) !< The basis functions b spanning the reference element … local
5479 EdgeBasis = 0.0d0
5866 EdgeBasis(i,1:3) = ReadyEdgeBasis(i,1:3)
5897 EdgeBasis(1,1) = -0.5d0
5899 EdgeBasis(1,1) = 0.5d0
5902 EdgeBasis(2,1) = 1.5d0 * u
5928 EdgeBasis(1,1) = -(3.0d0 + 3.0d0*Sqrt(3.0d0)*u - Sqrt(3.0d0)*v)/6.0d0
5929 EdgeBasis(1,2) = -(3.0d0 + Sqrt(3.0d0)*u - Sqrt(3.0d0)*v)/6.0d0
5932 EdgeBasis(2,1) = -(3.0d0 + Sqrt(3.0d0) - 3.0d0*(1.0d0 + Sqrt(3.0d0))*u - &
5934 EdgeBasis(2,2) = -(-3.0d0 - Sqrt(3.0d0) + u + Sqrt(3.0d0)*u + v + Sqrt(3.0d0)*v)/ &
5938 EdgeBasis(1,1) = (3.0d0 + Sqrt(3.0d0) - 3.0d0*(1.0d0 + Sqrt(3.0d0))*u - &
5940 EdgeBasis(1,2) = (-3.0d0 - Sqrt(3.0d0) + u + Sqrt(3.0d0)*u + v + Sqrt(3.0d0)*v)/ &
5944 EdgeBasis(2,1) = (3.0d0 + 3.0d0*Sqrt(3.0d0)*u - Sqrt(3.0d0)*v)/6.0d0
5945 EdgeBasis(2,2) = (3.0d0 + Sqrt(3.0d0)*u - Sqrt(3.0d0)*v)/6.0d0
5961 EdgeBasis(3,1) = ((3.0d0 + Sqrt(3.0d0))*v)/6.0d0
5962 … EdgeBasis(3,2) = -(-3.0d0 + Sqrt(3.0d0) + (-3.0d0 + Sqrt(3.0d0))*u + 2.0d0*Sqrt(3.0d0)*v)/6.0d0
5965 EdgeBasis(4,1) = ((-3.0d0 + Sqrt(3.0d0))*v)/6.0d0
5966 EdgeBasis(4,2) = -(2.0d0 + Sqrt(3.0d0) + (2.0d0 + Sqrt(3.0d0))*u - &
5970 EdgeBasis(3,1) = -((-3.0d0 + Sqrt(3.0d0))*v)/6.0d0
5971 EdgeBasis(3,2) = (2.0d0 + Sqrt(3.0d0) + (2.0d0 + Sqrt(3.0d0))*u - &
5975 EdgeBasis(4,1) = -((3.0d0 + Sqrt(3.0d0))*v)/6.0d0
5976 … EdgeBasis(4,2) = (-3.0d0 + Sqrt(3.0d0) + (-3.0d0 + Sqrt(3.0d0))*u + 2.0d0*Sqrt(3.0d0)*v)/6.0d0
5992 EdgeBasis(5,1) = ((-3.0d0 + Sqrt(3.0d0))*v)/6.0d0
5993 … EdgeBasis(5,2) = -(-3.0d0 - Sqrt(3.0d0) + (3.0d0 + Sqrt(3.0d0))*u + 2.0d0*Sqrt(3.0d0)*v)/6.0d0
5996 EdgeBasis(6,1) = ((3.0d0 + 2.0d0*Sqrt(3.0d0))*v)/(3.0d0*(1.0d0 + Sqrt(3.0d0)))
5997 EdgeBasis(6,2) = ((-1.0d0 + u + v + Sqrt(3.0d0)*v)/(3.0d0 + Sqrt(3.0d0)))
6000 EdgeBasis(5,1) = -((3.0d0 + 2.0d0*Sqrt(3.0d0))*v)/(3.0d0*(1.0d0 + Sqrt(3.0d0)))
6001 EdgeBasis(5,2) = -((-1.0d0 + u + v + Sqrt(3.0d0)*v)/(3.0d0 + Sqrt(3.0d0)))
6004 EdgeBasis(6,1) = -((-3.0d0 + Sqrt(3.0d0))*v)/6.0d0
6005 … EdgeBasis(6,2) = (-3.0d0 - Sqrt(3.0d0) + (3.0d0 + Sqrt(3.0d0))*u + 2.0d0*Sqrt(3.0d0)*v)/6.0d0
6023 EdgeBasis(1,1) = (3.0d0 - Sqrt(3.0d0)*v)/6.0d0
6024 EdgeBasis(1,2) = u/(2.0d0*Sqrt(3.0d0))
6027 EdgeBasis(1,:) = -EdgeBasis(1,:)
6031 EdgeBasis(2,1) = -(u*(-3.0d0 + Sqrt(3.0d0)*v))/2.0d0
6032 EdgeBasis(2,2) = (Sqrt(3.0d0)*u**2)/2.0d0
6041 EdgeBasis(4,1) = ((Sqrt(3.0d0) + Sqrt(3.0d0)*u - 3.0d0*v)*v)/4.0d0
6042 EdgeBasis(4,2) = (Sqrt(3.0d0)*(1.0d0 + u)*(-1.0d0 - u + Sqrt(3.0d0)*v))/4.0d0
6053 EdgeBasis(k,1) = -v/(2.0d0*Sqrt(3.0d0))
6054 EdgeBasis(k,2) = (1 + u)/(2.0d0*Sqrt(3.0d0))
6057 EdgeBasis(k,:) = -EdgeBasis(k,:)
6066 EdgeBasis(6,1) = (v*(-Sqrt(3.0d0) + Sqrt(3.0d0)*u + 3.0d0*v))/4.0d0
6067 EdgeBasis(6,2) = -(Sqrt(3.0d0)*(-1.0d0 + u)*(-1.0d0 + u + Sqrt(3.0d0)*v))/4.0d0
6078 EdgeBasis(k,1) = -v/(2.0d0*Sqrt(3.0d0))
6079 EdgeBasis(k,2) = (-1 + u)/(2.0d0*Sqrt(3.0d0))
6082 EdgeBasis(k,:) = -EdgeBasis(k,:)
6113 EdgeBasis(7,:) = D1 * WorkBasis(I1,:)
6115 EdgeBasis(8,:) = D2 * WorkBasis(I2,:)
6139 EdgeBasis(1,1) = 0.1D1 / 0.4D1 - v / 0.4D1
6142 EdgeBasis(1,:) = -EdgeBasis(1,:)
6145 EdgeBasis(2,1) = 0.3D1 * u * (0.1D1 / 0.4D1 - v / 0.4D1)
6154 EdgeBasis(3,2) = 0.1D1 / 0.4D1 + u / 0.4D1
6157 EdgeBasis(3,:) = -EdgeBasis(3,:)
6160 EdgeBasis(4,2) = 0.3D1 * v * (0.1D1 / 0.4D1 + u / 0.4D1)
6169 EdgeBasis(5,1) = -0.1D1 / 0.4D1 - v / 0.4D1
6172 EdgeBasis(5,:) = -EdgeBasis(5,:)
6175 EdgeBasis(6,1) = -0.3D1 * u * (-0.1D1 / 0.4D1 - v / 0.4D1)
6184 EdgeBasis(7,2) = -0.1D1 / 0.4D1 + u / 0.4D1
6187 EdgeBasis(7,:) = -EdgeBasis(7,:)
6190 EdgeBasis(8,2) = -0.3D1 * v * (-0.1D1 / 0.4D1 + u / 0.4D1)
6224 EdgeBasis(9,:) = D1 * WorkBasis(2*(I1-1)+1,:)
6226 EdgeBasis(10,:) = WorkBasis(2*(I1-1)+2,:)
6228 EdgeBasis(11,:) = D2 * WorkBasis(2*(I2-1)+1,:)
6230 EdgeBasis(12,:) = WorkBasis(2*(I2-1)+2,:)
6245 EdgeBasis(1,1) = ((-1.0d0 + v)*v)/4.0d0
6246 EdgeBasis(1,2) = 0.0d0
6249 EdgeBasis(1,:) = -EdgeBasis(1,:)
6259 EdgeBasis(2,1) = 0.0d0
6260 EdgeBasis(2,2) = (u*(1.0d0 + u))/4.0d0
6263 EdgeBasis(2,:) = -EdgeBasis(2,:)
6273 EdgeBasis(3,1) = -(v*(1.0d0 + v))/4.0d0
6274 EdgeBasis(3,2) = 0.0d0
6277 EdgeBasis(3,:) = -EdgeBasis(3,:)
6287 EdgeBasis(4,1) = 0.0d0
6288 EdgeBasis(4,2) = -((-1 + u)*u)/4.0d0
6291 EdgeBasis(4,:) = -EdgeBasis(4,:)
6324 EdgeBasis(5,:) = D1 * WorkBasis(I1,:)
6326 EdgeBasis(6,:) = D2 * WorkBasis(I2,:)
6349 … EdgeBasis(1,1) = -(6.0d0 + 6.0d0*Sqrt(3.0d0)*u - 2.0d0*Sqrt(3.0d0)*v - Sqrt(6.0d0)*w)/12.0d0
6350 … EdgeBasis(1,2) = -(6.0d0 + 2.0d0*Sqrt(3.0d0)*u - 2.0d0*Sqrt(3.0d0)*v - Sqrt(6.0d0)*w)/12.0d0
6351 … EdgeBasis(1,3) = -(3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0)*u - Sqrt(6.0d0)*v - Sqrt(3.0d0)*w)/12.0d0
6356 EdgeBasis(2,1) = (-6.0d0 - 2.0d0*Sqrt(3.0d0) + 6.0d0*(1.0d0 + Sqrt(3.0d0))*u + &
6358 EdgeBasis(2,2) = -(-6.0d0 - 2.0d0*Sqrt(3.0d0) + 2.0d0*(1.0d0 + Sqrt(3.0d0))*u + &
6360 … EdgeBasis(2,3) = -(-3.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) + (Sqrt(2.0d0) + Sqrt(6.0d0))*u + &
6366 EdgeBasis(1,1) = -(-6.0d0 - 2.0d0*Sqrt(3.0d0) + 6.0d0*(1.0d0 + Sqrt(3.0d0))*u + &
6368 EdgeBasis(1,2) = (-6.0d0 - 2.0d0*Sqrt(3.0d0) + 2.0d0*(1.0d0 + Sqrt(3.0d0))*u + &
6370 … EdgeBasis(1,3) = (-3.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) + (Sqrt(2.0d0) + Sqrt(6.0d0))*u + &
6376 … EdgeBasis(2,1) = (6.0d0 + 6.0d0*Sqrt(3.0d0)*u - 2.0d0*Sqrt(3.0d0)*v - Sqrt(6.0d0)*w)/12.0d0
6377 … EdgeBasis(2,2) = (6.0d0 + 2.0d0*Sqrt(3.0d0)*u - 2.0d0*Sqrt(3.0d0)*v - Sqrt(6.0d0)*w)/12.0d0
6378 … EdgeBasis(2,3) = (3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0)*u - Sqrt(6.0d0)*v - Sqrt(3.0d0)*w)/12.0d0
6396 EdgeBasis(3,1) = (3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w)/24.0d0
6397 EdgeBasis(3,2) = -(4.0d0*(-3.0d0 + Sqrt(3.0d0))*u + 8.0d0*Sqrt(3.0d0)*v + &
6399 … EdgeBasis(3,3) = -(3.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) - Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6405 EdgeBasis(4,1) = (-3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w)/24.0d0
6406 EdgeBasis(4,2) = (-4.0d0*(2.0d0 + Sqrt(3.0d0))*u + 4.0d0*(1.0d0 + Sqrt(3.0d0))*v + &
6408 … EdgeBasis(4,3) = -(-2.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) - Sqrt(2.0d0)*(2.0d0 + Sqrt(3.0d0))*u + &
6414 EdgeBasis(3,1) = -(-3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w)/24.0d0
6415 EdgeBasis(3,2) = -(-4.0d0*(2.0d0 + Sqrt(3.0d0))*u + 4.0d0*(1.0d0 + Sqrt(3.0d0))*v + &
6417 … EdgeBasis(3,3) = (-2.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) - Sqrt(2.0d0)*(2.0d0 + Sqrt(3.0d0))*u + &
6423 EdgeBasis(4,1) = -((3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w))/24.0d0
6424 EdgeBasis(4,2) = (4.0d0*(-3.0d0 + Sqrt(3.0d0))*u + 8.0d0*Sqrt(3.0d0)*v + &
6426 … EdgeBasis(4,3) = (3.0d0*Sqrt(2.0d0) - Sqrt(6.0d0) - Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6445 EdgeBasis(5,1) = ((-3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w))/24.0d0
6446 EdgeBasis(5,2) = -(4.0d0*(3.0d0 + Sqrt(3.0d0))*u + 8.0d0*Sqrt(3.0d0)*v + &
6448 … EdgeBasis(5,3) = -(3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) - Sqrt(2.0d0)*(3.0d0 + Sqrt(3.0d0))*u + &
6454 EdgeBasis(6,1) = ((3.0d0 + 2.0d0*Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w))/ &
6456 … EdgeBasis(6,2) = -(4.0d0 - 4.0d0*u - 4.0d0*(1.0d0 + Sqrt(3.0d0))*v + Sqrt(2.0d0)*w)/ &
6458 … EdgeBasis(6,3) = -(-Sqrt(2.0d0) + Sqrt(2.0d0)*u - Sqrt(2.0d0)*(2.0d0 + Sqrt(3.0d0))*v + &
6464 EdgeBasis(5,1) = -((3.0d0 + 2.0d0*Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w))/ &
6466 EdgeBasis(5,2) = (4.0d0 - 4.0d0*u - 4.0d0*(1.0d0 + Sqrt(3.0d0))*v + Sqrt(2.0d0)*w)/ &
6468 … EdgeBasis(5,3) = (-Sqrt(2.0d0) + Sqrt(2.0d0)*u - Sqrt(2.0d0)*(2.0d0 + Sqrt(3.0d0))*v + &
6474 EdgeBasis(6,1) = -((-3.0d0 + Sqrt(3.0d0))*(4.0d0*v - Sqrt(2.0d0)*w))/24.0d0
6475 EdgeBasis(6,2) = (4.0d0*(3.0d0 + Sqrt(3.0d0))*u + 8.0d0*Sqrt(3.0d0)*v + &
6477 … EdgeBasis(6,3) = (3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) - Sqrt(2.0d0)*(3.0d0 + Sqrt(3.0d0))*u + &
6496 EdgeBasis(7,1) = -((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6497 EdgeBasis(7,2) = -((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6498 … EdgeBasis(7,3) = -(-3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) - Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6504 EdgeBasis(8,1) = -((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6505 EdgeBasis(8,2) = -((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6506 … EdgeBasis(8,3) = -((-3.0d0 + Sqrt(3.0d0))*w - (Sqrt(2.0d0)*(3.0d0 + 2.0d0*Sqrt(3.0d0))* &
6513 EdgeBasis(7,1) = ((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6514 EdgeBasis(7,2) = ((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6515 … EdgeBasis(7,3) = ((-3.0d0 + Sqrt(3.0d0))*w - (Sqrt(2.0d0)*(3.0d0 + 2.0d0*Sqrt(3.0d0))* &
6522 EdgeBasis(8,1) = ((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6523 EdgeBasis(8,2) = ((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6524 … EdgeBasis(8,3) = (-3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) - Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6543 EdgeBasis(9,1) = ((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6544 EdgeBasis(9,2) = -((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6545 … EdgeBasis(9,3) = -(-3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) + Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6551 EdgeBasis(10,1) = ((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6552 EdgeBasis(10,2) = -((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6553 … EdgeBasis(10,3) = -((-3.0d0 + Sqrt(3.0d0))*w - (Sqrt(2.0d0)*(3.0d0 + 2.0d0*Sqrt(3.0d0))*&
6560 EdgeBasis(9,1) = -((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6561 EdgeBasis(9,2) = ((-3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6562 … EdgeBasis(9,3) = ((-3.0d0 + Sqrt(3.0d0))*w - (Sqrt(2.0d0)*(3.0d0 + 2.0d0*Sqrt(3.0d0))*&
6569 EdgeBasis(10,1) = -((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(2.0d0))
6570 EdgeBasis(10,2) = ((3.0d0 + Sqrt(3.0d0))*w)/(4.0d0*Sqrt(6.0d0))
6571 … EdgeBasis(10,3) = (-3.0d0*Sqrt(2.0d0) + Sqrt(6.0d0) + Sqrt(2.0d0)*(-3.0d0 + Sqrt(3.0d0))*u + &
6590 EdgeBasis(11,1) = 0.0d0
6591 EdgeBasis(11,2) = ((1.0d0 + Sqrt(3.0d0))*w)/(2.0d0*Sqrt(2.0d0))
6592 EdgeBasis(11,3) = -(6.0d0*Sqrt(2.0d0)*v - 4.0d0*Sqrt(6.0d0)*v - 3.0d0*w + &
6598 EdgeBasis(12,1) = 0.0d0
6599 EdgeBasis(12,2) = ((-3.0d0 + Sqrt(3.0d0))*w)/(2.0d0*Sqrt(6.0d0))
6600 EdgeBasis(12,3) = -((Sqrt(2.0d0) + Sqrt(6.0d0))*v - Sqrt(3.0d0)*w)/4.0d0
6605 EdgeBasis(11,1) = 0.0d0
6606 EdgeBasis(11,2) = -((-3.0d0 + Sqrt(3.0d0))*w)/(2.0d0*Sqrt(6.0d0))
6607 EdgeBasis(11,3) = ((Sqrt(2.0d0) + Sqrt(6.0d0))*v - Sqrt(3.0d0)*w)/4.0d0
6612 EdgeBasis(12,1) = 0.0d0
6613 EdgeBasis(12,2) = -((1.0d0 + Sqrt(3.0d0))*w)/(2.0d0*Sqrt(2.0d0))
6614 EdgeBasis(12,3) = (6.0d0*Sqrt(2.0d0)*v - 4.0d0*Sqrt(6.0d0)*v - 3.0d0*w + &
6635 EdgeBasis(1,1) = (6.0d0 - 2.0d0*Sqrt(3.0d0)*v - Sqrt(6.0d0)*w)/24.0d0
6636 EdgeBasis(1,2) = u/(4.0d0*Sqrt(3.0d0))
6637 EdgeBasis(1,3) = u/(4.0d0*Sqrt(6.0d0))
6642 EdgeBasis(1,:) = -EdgeBasis(1,:)
6646 EdgeBasis(2,1) = -(u*(-6.0d0 + 2.0d0*Sqrt(3.0d0)*v + Sqrt(6.0d0)*w))/4.0d0
6647 EdgeBasis(2,2) = (Sqrt(3.0d0)*u**2)/2.0d0
6648 EdgeBasis(2,3) = (Sqrt(1.5d0)*u**2)/2.0d0
6659 … EdgeBasis(4,1) = ((Sqrt(3.0d0) + Sqrt(3.0d0)*u - 3.0d0*v)*(4.0d0*v - Sqrt(2.0d0)*w))/16.0d0
6660 EdgeBasis(4,2) = -((1.0d0 + u - Sqrt(3.0d0)*v)*&
6662 EdgeBasis(4,3) = -((Sqrt(3.0d0) + Sqrt(3.0d0)*u - 3.0d0*v)*&
6677 EdgeBasis(k,1) = (-4.0d0*v + Sqrt(2.0d0)*w)/(16.0d0*Sqrt(3.0d0))
6678 EdgeBasis(k,2) = (4.0d0*Sqrt(3.0d0) + 4.0d0*Sqrt(3.0d0)*u - 3.0d0*Sqrt(2.0d0)*w)/48.0d0
6679 EdgeBasis(k,3) = -(Sqrt(3.0d0) + Sqrt(3.0d0)*u - 3.0d0*v)/(24.0d0*Sqrt(2.0d0))
6684 EdgeBasis(k,:) = -EdgeBasis(k,:)
6693 EdgeBasis(6,1) = ((-Sqrt(3.0d0) + Sqrt(3.0d0)*u + 3.0d0*v)*&
6695 EdgeBasis(6,2) = -((-1.0d0 + u + Sqrt(3.0d0)*v)*&
6697 EdgeBasis(6,3) = ((-Sqrt(3.0d0) + Sqrt(3.0d0)*u + 3.0d0*v)*&
6712 EdgeBasis(k,1) = (-4.0d0*v + Sqrt(2.0d0)*w)/(16.0d0*Sqrt(3.0d0))
6713 … EdgeBasis(k,2) = (-4.0d0*Sqrt(3.0d0) + 4.0d0*Sqrt(3.0d0)*u + 3.0d0*Sqrt(2.0d0)*w)/48.0d0
6714 EdgeBasis(k,3) = (Sqrt(6.0d0) - Sqrt(6.0d0)*u - 3.0d0*Sqrt(2.0d0)*v)/48.0d0
6719 EdgeBasis(k,:) = -EdgeBasis(k,:)
6728 … EdgeBasis(8,1) = (3.0d0*w*(-Sqrt(6.0d0) + Sqrt(6.0d0)*u + Sqrt(2.0d0)*v + 4.0d0*w))/16.0d0
6729 EdgeBasis(8,2) = (w*(-3.0d0*Sqrt(2.0d0) + 3.0d0*Sqrt(2.0d0)*u + Sqrt(6.0d0)*v + &
6731 EdgeBasis(8,3) = -((-Sqrt(3.0d0) + Sqrt(3.0d0)*u + v)*&
6747 EdgeBasis(k,1) = (Sqrt(1.5d0)*w)/8.0d0
6748 EdgeBasis(k,2) = w/(8.0d0*Sqrt(2.0d0))
6749 EdgeBasis(k,3) = (Sqrt(6.0d0) - Sqrt(6.0d0)*u - Sqrt(2.0d0)*v)/16.0d0
6754 EdgeBasis(k,:) = -EdgeBasis(k,:)
6763 … EdgeBasis(10,1) = (3.0d0*(Sqrt(6.0d0) + Sqrt(6.0d0)*u - Sqrt(2.0d0)*v - 4.0d0*w)*w)/16.0d0
6764 EdgeBasis(10,2) = (w*(-3.0d0*Sqrt(2.0d0) - 3.0d0*Sqrt(2.0d0)*u + &
6766 EdgeBasis(10,3) = ((Sqrt(6.0d0) + Sqrt(6.0d0)*u - Sqrt(2.0d0)*v)*&
6782 EdgeBasis(k,1) = -(Sqrt(1.5d0)*w)/8.0d0
6783 EdgeBasis(k,2) = w/(8.0d0*Sqrt(2.0d0))
6784 EdgeBasis(k,3) = (Sqrt(6.0d0) + Sqrt(6.0d0)*u - Sqrt(2.0d0)*v)/16.0d0
6789 EdgeBasis(k,:) = -EdgeBasis(k,:)
6798 EdgeBasis(12,1) = 0.0d0
6799 EdgeBasis(12,2) = (Sqrt(3.0d0)*(Sqrt(2.0d0)*v - 2.0d0*w)*w)/4.0d0
6800 EdgeBasis(12,3) = (Sqrt(1.5d0)*v*(-v + Sqrt(2.0d0)*w))/2.0d0
6814 EdgeBasis(k,1) = 0.0d0
6815 EdgeBasis(k,2) = -w/(4.0d0*Sqrt(2.0d0))
6816 EdgeBasis(k,3) = v/(4.0d0*Sqrt(2.0d0))
6821 EdgeBasis(k,:) = -EdgeBasis(k,:)
6831 EdgeBasis(2*(k-1)+1,:) = 2.0d0 * EdgeBasis(2*(k-1)+1,:)
6836 EdgeBasis(k,:) = 2.0d0 * EdgeBasis(k,:)
6890 EdgeBasis(13,:) = D1 * WorkBasis(I1,:)
6892 EdgeBasis(14,:) = D2 * WorkBasis(I2,:)
6935 EdgeBasis(15,:) = D1 * WorkBasis(I1,:)
6937 EdgeBasis(16,:) = D2 * WorkBasis(I2,:)
6980 EdgeBasis(17,:) = D1 * WorkBasis(I1,:)
6982 EdgeBasis(18,:) = D2 * WorkBasis(I2,:)
7025 EdgeBasis(19,:) = D1 * WorkBasis(I1,:)
7027 EdgeBasis(20,:) = D2 * WorkBasis(I2,:)
7065 EdgeBasis(1,1) = 0.1D1 / 0.4D1 - v / 0.4D1 - w * sqrt(0.2D1) / 0.8D1
7066 EdgeBasis(1,2) = 0.0d0
7067 EdgeBasis(1,3) = sqrt(0.2D1) * u * (w * sqrt(0.2D1) + 2.0D0 * v - 0.2D1) / &
7074 EdgeBasis(1,:) = -EdgeBasis(1,:)
7079 EdgeBasis(2,1:3) = 3.0d0 * u * EdgeBasis(1,1:3)
7095 EdgeBasis(k,1) = 0.0d0
7096 EdgeBasis(k,2) = 0.1D1 / 0.4D1 + u / 0.4D1 - w * sqrt(0.2D1) / 0.8D1
7097 EdgeBasis(k,3) = sqrt(0.2D1) * v * (w * sqrt(0.2D1) - 2.0D0 * u - 0.2D1) / &
7104 EdgeBasis(k,:) = -EdgeBasis(k,:)
7109 EdgeBasis(k+1,1:3) = 3.0d0 * v * EdgeBasis(k,1:3)
7125 EdgeBasis(k,1) = 0.1D1 / 0.4D1 + v / 0.4D1 - w * sqrt(0.2D1) / 0.8D1
7126 EdgeBasis(k,2) = 0.0d0
7127 EdgeBasis(k,3) = sqrt(0.2D1) * u * (w * sqrt(0.2D1) - 2.0D0 * v - 0.2D1) / &
7135 EdgeBasis(k,:) = -EdgeBasis(k,:)
7140 EdgeBasis(k+1,1:3) = 3.0d0 * u * EdgeBasis(k,1:3)
7157 EdgeBasis(k,1) = 0.0d0
7158 EdgeBasis(k,2) = 0.1D1 / 0.4D1 - u / 0.4D1 - w * sqrt(0.2D1) / 0.8D1
7159 EdgeBasis(k,3) = sqrt(0.2D1) * v * (w * sqrt(0.2D1) + 2.0D0 * u - 0.2D1) / &
7167 EdgeBasis(k,:) = -EdgeBasis(k,:)
7172 EdgeBasis(k+1,1:3) = 3.0d0 * v * EdgeBasis(k,1:3)
7189 EdgeBasis(k,1) = w * sqrt(0.2D1) * (w * sqrt(0.2D1) + 2.0D0 * v - 0.2D1) / &
7191 EdgeBasis(k,2) = w * sqrt(0.2D1) * (w * sqrt(0.2D1) + 2.0D0 * u - 0.2D1) / &
7193 EdgeBasis(k,3) = -sqrt(0.2D1)/ 0.4D1 * (0.2D1 * sqrt(0.2D1) * u * v * w - &
7204 EdgeBasis(k,:) = -EdgeBasis(k,:)
7209 … EdgeBasis(k+1,1:3) = 3.0d0 * EdgeSign(k) * EdgeBasis(k,1:3) * ( LBasis(5)-LBasis(1)+LBasis(3) )
7236 EdgeBasis(k,1) = -w * sqrt(0.2D1) * (w * sqrt(0.2D1) + 2.0D0 * v - 0.2D1) / &
7238 EdgeBasis(k,2) = w * sqrt(0.2D1) * (w * sqrt(0.2D1) - 2.0D0 * u - 0.2D1) / &
7240 EdgeBasis(k,3) = sqrt(0.2D1)/ 0.4D1 * (0.2D1 * sqrt(0.2D1) * u * v * w - 0.2D1 * &
7250 EdgeBasis(k,:) = -EdgeBasis(k,:)
7255 … EdgeBasis(k+1,1:3) = 3.0d0 * EdgeSign(k) * EdgeBasis(k,1:3) * ( LBasis(5)-LBasis(2)+LBasis(4) )
7282 EdgeBasis(k,1) = -w * sqrt(0.2D1)/ 0.8D1 * (w * sqrt(0.2D1) - 2.0D0 * v - 0.2D1) / &
7284 EdgeBasis(k,2) = -w * sqrt(0.2D1) / 0.8D1 * (w * sqrt(0.2D1) - 2.0D0 * u - 0.2D1) / &
7286 EdgeBasis(k,3) = -sqrt(0.2D1)/ 0.4D1 * (0.2D1 * sqrt(0.2D1) * u * v * w + 0.2D1 * &
7297 EdgeBasis(k,:) = -EdgeBasis(k,:)
7302 … EdgeBasis(k+1,1:3) = 3.0d0 * EdgeSign(k) * EdgeBasis(k,1:3) * ( LBasis(5)-LBasis(3)+LBasis(1) )
7329 EdgeBasis(k,1) = w * sqrt(0.2D1) / 0.8D1 * (w * sqrt(0.2D1) - 2.0D0 * v - 0.2D1) / &
7331 EdgeBasis(k,2) = -w * sqrt(0.2D1) / 0.8D1 * (w * sqrt(0.2D1) + 2.0D0 * u - 0.2D1) / &
7333 EdgeBasis(k,3) = sqrt(0.2D1) / 0.4D1 * (0.2D1 * sqrt(0.2D1) * u * v * w + &
7343 EdgeBasis(k,:) = -EdgeBasis(k,:)
7348 … EdgeBasis(k+1,1:3) = 3.0d0 * EdgeSign(k) * EdgeBasis(k,1:3) * ( LBasis(5)-LBasis(4)+LBasis(2) )
7369 WorkBasis(1,1:3) = 2.0d0 * ( EdgeSign(1) * EdgeBasis(1,1:3) * Beta(4) + &
7370 EdgeSign(5) * EdgeBasis(5,1:3) * Beta(2) ) / (1.0d0 - LBasis(5))
7383 WorkBasis(3,1:3) = 2.0d0 * ( EdgeSign(3) * EdgeBasis(3,1:3) * Beta(1) + &
7384 EdgeSign(7) * EdgeBasis(7,1:3) * Beta(3) ) / (1.0d0 - LBasis(5))
7410 EdgeBasis(17,:) = D1 * WorkBasis(2*(I1-1)+1,:)
7412 EdgeBasis(18,:) = WorkBasis(2*(I1-1)+2,:)
7414 EdgeBasis(19,:) = D2 * WorkBasis(2*(I2-1)+1,:)
7416 EdgeBasis(20,:) = WorkBasis(2*(I2-1)+2,:)
7435 WorkBasis(1,1:3) = LBasis(5) * EdgeSign(1) * EdgeBasis(1,1:3)
7442 WorkBasis(2,1:3) = Beta(3) * EdgeSign(9) * EdgeBasis(9,1:3)
7455 WorkBasis(3,1:3) = Beta(1) * EdgeSign(11) * EdgeBasis(11,1:3)
7468 EdgeBasis(21,:) = D1 * WorkBasis(I1,:)
7470 EdgeBasis(22,:) = D2 * WorkBasis(I2,:)
7488 WorkBasis(1,1:3) = LBasis(5) * EdgeSign(3) * EdgeBasis(3,1:3)
7494 WorkBasis(2,1:3) = Beta(4) * EdgeSign(11) * EdgeBasis(11,1:3)
7506 WorkBasis(3,1:3) = Beta(2) * EdgeSign(13) * EdgeBasis(13,1:3)
7518 EdgeBasis(23,:) = D1 * WorkBasis(I1,:)
7520 EdgeBasis(24,:) = D2 * WorkBasis(I2,:)
7538 WorkBasis(1,1:3) = -LBasis(5) * EdgeSign(5) * EdgeBasis(5,1:3)
7544 WorkBasis(2,1:3) = Beta(1) * EdgeSign(13) * EdgeBasis(13,1:3)
7556 WorkBasis(3,1:3) = Beta(3) * EdgeSign(15) * EdgeBasis(15,1:3)
7568 EdgeBasis(25,:) = D1 * WorkBasis(I1,:)
7570 EdgeBasis(26,:) = D2 * WorkBasis(I2,:)
7588 WorkBasis(1,1:3) = -LBasis(5) * EdgeSign(7) * EdgeBasis(7,1:3)
7594 WorkBasis(2,1:3) = Beta(2) * EdgeSign(15) * EdgeBasis(15,1:3)
7606 WorkBasis(3,1:3) = Beta(4) * EdgeSign(9) * EdgeBasis(9,1:3)
7618 EdgeBasis(27,:) = D1 * WorkBasis(I1,:)
7620 EdgeBasis(28,:) = D2 * WorkBasis(I2,:)
7626 EdgeBasis(29,1:3) = LBasis(5) * Beta(4) * EdgeSign(1) * EdgeBasis(1,1:3)
7633 EdgeBasis(30,1:3) = LBasis(5) * Beta(3) * EdgeSign(7) * EdgeBasis(7,1:3)
7640 EdgeBasis(31,1:3) = Beta(3) * Beta(4) * EdgeSign(9) * EdgeBasis(9,1:3)
7667 EdgeBasis(1,1) = (v*(-1 + (2*v)/(2 - Sqrt(2.0d0)*w)))/4.0d0
7668 EdgeBasis(1,2) = 0.0d0
7669 … EdgeBasis(1,3) = (u*v*(-Sqrt(2.0d0) + Sqrt(2.0d0)*v + w))/(2.0d0*(-2 + Sqrt(2.0d0)*w)**2)
7674 EdgeBasis(1,:) = -EdgeBasis(1,:)
7684 EdgeBasis(2,1) = 0.0d0
7685 EdgeBasis(2,2) = (u*(1 + (2*u)/(2 - Sqrt(2.0d0)*w)))/4.0d0
7686 … EdgeBasis(2,3) = (u*v*(Sqrt(2.0d0) + Sqrt(2.0d0)*u - w))/(2.0d0*(-2 + Sqrt(2.0d0)*w)**2)
7691 EdgeBasis(2,:) = -EdgeBasis(2,:)
7701 EdgeBasis(3,1) = (v*(1 + (2*v)/(2 - Sqrt(2.0d0)*w)))/4.0d0
7702 EdgeBasis(3,2) = 0.0d0
7703 … EdgeBasis(3,3) = (u*v*(Sqrt(2.0d0) + Sqrt(2.0d0)*v - w))/(2.0d0*(-2 + Sqrt(2.0d0)*w)**2)
7708 EdgeBasis(3,:) = -EdgeBasis(3,:)
7718 EdgeBasis(4,1) = 0.0d0
7719 EdgeBasis(4,2) = (u*(-1 + (2*u)/(2 - Sqrt(2.0d0)*w)))/4.0d0
7720 … EdgeBasis(4,3) = (u*v*(-Sqrt(2.0d0) + Sqrt(2.0d0)*u + w))/(2.0d0*(-2 + Sqrt(2.0d0)*w)**2)
7725 EdgeBasis(4,:) = -EdgeBasis(4,:)
7735 EdgeBasis(5,1) = (w*(-Sqrt(2.0d0) + Sqrt(2.0d0)*v + w))/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7736 EdgeBasis(5,2) = (w*(-Sqrt(2.0d0) + Sqrt(2.0d0)*u + w))/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7737 … EdgeBasis(5,3) = (u*(-2*Sqrt(2.0d0) + 2*v*(Sqrt(2.0d0) - 2*w) + 4*w - Sqrt(2.0d0)*w**2) - &
7745 EdgeBasis(5,:) = -EdgeBasis(5,:)
7755 EdgeBasis(6,1) = (w*(-Sqrt(2.0d0) + Sqrt(2.0d0)*v + w))/(8.0d0 - 4*Sqrt(2.0d0)*w)
7756 EdgeBasis(6,2) = (w*(-Sqrt(2.0d0) - Sqrt(2.0d0)*u + w))/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7757 EdgeBasis(6,3) = (-((-1 + v)*(2*Sqrt(2.0d0) - 4*w + Sqrt(2.0d0)*w**2)) + &
7766 EdgeBasis(6,:) = -EdgeBasis(6,:)
7776 EdgeBasis(7,1) = ((Sqrt(2.0d0) + Sqrt(2.0d0)*v - w)*w)/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7777 EdgeBasis(7,2) = ((Sqrt(2.0d0) + Sqrt(2.0d0)*u - w)*w)/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7778 EdgeBasis(7,3) = ((1 + v)*(2*Sqrt(2.0d0) - 4*w + Sqrt(2.0d0)*w**2) + &
7787 EdgeBasis(7,:) = -EdgeBasis(7,:)
7797 EdgeBasis(8,1) = (w*(-Sqrt(2.0d0) - Sqrt(2.0d0)*v + w))/(-8.0d0 + 4*Sqrt(2.0d0)*w)
7798 EdgeBasis(8,2) = (w*(-Sqrt(2.0d0) + Sqrt(2.0d0)*u + w))/(8.0d0 - 4*Sqrt(2.0d0)*w)
7799 EdgeBasis(8,3) = ((1 + v)*(2*Sqrt(2.0d0) - 4*w + Sqrt(2.0d0)*w**2) - &
7808 EdgeBasis(8,:) = -EdgeBasis(8,:)
7847 EdgeBasis(9,:) = D1 * WorkBasis(I1,:)
7849 EdgeBasis(10,:) = D2 * WorkBasis(I2,:)
7892 EdgeBasis(1,1:2) = WorkBasis(1,1:2) * h1
7896 EdgeBasis(2,1:2) = WorkBasis(2,1:2) * h1
7901 EdgeBasis(1,1:2) = -EdgeBasis(1,1:2)
7911 EdgeBasis(7,1:2) = WorkBasis(1,1:2) * h2
7915 EdgeBasis(8,1:2) = WorkBasis(2,1:2) * h2
7920 EdgeBasis(7,1:2) = -EdgeBasis(7,1:2)
7942 EdgeBasis(3,1:2) = WorkBasis(1,1:2) * h1
7946 EdgeBasis(4,1:2) = WorkBasis(2,1:2) * h1
7951 EdgeBasis(3,1:2) = -EdgeBasis(3,1:2)
7961 EdgeBasis(9,1:2) = WorkBasis(1,1:2) * h2
7965 EdgeBasis(10,1:2) = WorkBasis(2,1:2) * h2
7970 EdgeBasis(9,1:2) = -EdgeBasis(9,1:2)
7992 EdgeBasis(5,1:2) = WorkBasis(1,1:2) * h1
7996 EdgeBasis(6,1:2) = WorkBasis(2,1:2) * h1
8001 EdgeBasis(5,1:2) = -EdgeBasis(5,1:2)
8011 EdgeBasis(11,1:2) = WorkBasis(1,1:2) * h2
8015 EdgeBasis(12,1:2) = WorkBasis(2,1:2) * h2
8020 EdgeBasis(11,1:2) = -EdgeBasis(11,1:2)
8036 EdgeBasis(12+(q-1)*2+1,3) = 0.5d0 * TriangleNodalPBasis(q, u, v)
8039 EdgeBasis(12+(q-1)*2+2,3) = 3.0d0 * EdgeBasis(12+(q-1)*2+1,3) * w
8044 EdgeBasis(12+(q-1)*2+1,3) = -EdgeBasis(12+(q-1)*2+1,3)
8074 EdgeBasis(19,1:2) = D1 * WorkBasis(I1,1:2) * h1
8079 EdgeBasis(20,1:2) = D2 * WorkBasis(I2,1:2) * h1
8099 EdgeBasis(21,1:2) = D1 * WorkBasis(I1,1:2) * h2
8104 EdgeBasis(22,1:2) = D2 * WorkBasis(I2,1:2) * h2
8146 EdgeBasis(23,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8148 EdgeBasis(24,:) = WorkBasis(2*(I1-1)+2,:)
8150 EdgeBasis(25,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8152 EdgeBasis(26,:) = WorkBasis(2*(I2-1)+2,:)
8192 EdgeBasis(27,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8194 EdgeBasis(28,:) = WorkBasis(2*(I1-1)+2,:)
8196 EdgeBasis(29,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8198 EdgeBasis(30,:) = WorkBasis(2*(I2-1)+2,:)
8238 EdgeBasis(31,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8240 EdgeBasis(32,:) = WorkBasis(2*(I1-1)+2,:)
8242 EdgeBasis(33,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8244 EdgeBasis(34,:) = WorkBasis(2*(I2-1)+2,:)
8250 EdgeBasis(35,1) = (v*(1.0d0 + u - v/Sqrt(3.0d0)))/(4.0d0*Sqrt(3.0d0)) * h3
8251 … EdgeBasis(35,2) = ((-1.0d0 + u)*(-3.0d0 - 3.0d0*u + Sqrt(3.0d0)*v))/(12.0d0*Sqrt(3.0d0)) * h3
8252 CurlBasis(35,1) = -EdgeBasis(35,2)/h3 * dh3
8253 CurlBasis(35,2) = EdgeBasis(35,1)/h3 * dh3
8256 EdgeBasis(36,1) = (v*(-3.0d0 + 3.0d0*u + Sqrt(3.0d0)*v))/(12.0d0*Sqrt(3.0d0)) * h3
8257 … EdgeBasis(36,2) = -((1.0d0 + u)*(-3.0d0 + 3.0d0*u + Sqrt(3.0d0)*v))/(12.0d0*Sqrt(3.0d0)) * h3
8258 CurlBasis(36,1) = -EdgeBasis(36,2)/h3 * dh3
8259 CurlBasis(36,2) = EdgeBasis(36,1)/h3 * dh3
8275 EdgeBasis(1,1) = -((-3.0d0 + Sqrt(3.0d0)*v)*(-1.0d0 + w)*w)/12.0d0
8276 EdgeBasis(1,2) = (u*(-1.0d0 + w)*w)/(4.0d0*Sqrt(3.0d0))
8277 EdgeBasis(1,3) = 0.0d0
8282 EdgeBasis(1,:) = -EdgeBasis(1,:)
8292 EdgeBasis(2,1) = -(v*(-1.0d0 + w)*w)/(4.0d0*Sqrt(3.0d0))
8293 EdgeBasis(2,2) = ((1.0d0 + u)*(-1.0d0 + w)*w)/(4.0d0*Sqrt(3.0d0))
8294 EdgeBasis(2,3) = 0.0d0
8299 EdgeBasis(2,:) = -EdgeBasis(2,:)
8309 EdgeBasis(3,1) = -(v*(-1.0d0 + w)*w)/(4.0d0*Sqrt(3.0d0))
8310 EdgeBasis(3,2) = ((-1.0d0 + u)*(-1.0d0 + w)*w)/(4.0d0*Sqrt(3.0d0))
8311 EdgeBasis(3,3) = 0.0d0
8316 EdgeBasis(3,:) = -EdgeBasis(3,:)
8326 EdgeBasis(4,1) = -((-3.0d0 + Sqrt(3.0d0)*v)*w*(1.0d0 + w))/12.0d0
8327 EdgeBasis(4,2) = (u*w*(1.0d0 + w))/(4.0d0*Sqrt(3.0d0))
8328 EdgeBasis(4,3) = 0.0d0
8333 EdgeBasis(4,:) = -EdgeBasis(4,:)
8343 EdgeBasis(5,1) = -(v*w*(1.0d0 + w))/(4.0d0*Sqrt(3.0d0))
8344 EdgeBasis(5,2) = ((1.0d0 + u)*w*(1.0d0 + w))/(4.0d0*Sqrt(3.0d0))
8345 EdgeBasis(5,3) = 0.0d0
8350 EdgeBasis(5,:) = -EdgeBasis(5,:)
8360 EdgeBasis(6,1) = -(v*w*(1.0d0 + w))/(4.0d0*Sqrt(3.0d0))
8361 EdgeBasis(6,2) = ((-1.0d0 + u)*w*(1.0d0 + w))/(4.0d0*Sqrt(3.0d0))
8362 EdgeBasis(6,3) = 0.0d0
8367 EdgeBasis(6,:) = -EdgeBasis(6,:)
8377 EdgeBasis(7,1) = 0.0d0
8378 EdgeBasis(7,2) = 0.0d0
8379 EdgeBasis(7,3) = (3*u**2 + v*(-Sqrt(3.0d0) + v) + u*(-3.0d0 + 2*Sqrt(3.0d0)*v))/12.0d0
8384 EdgeBasis(7,:) = -EdgeBasis(7,:)
8394 EdgeBasis(8,1) = 0.0d0
8395 EdgeBasis(8,2) = 0.0d0
8396 EdgeBasis(8,3) = (3*u**2 + v*(-Sqrt(3.0d0) + v) + u*(3.0d0 - 2*Sqrt(3.0d0)*v))/12.0d0
8401 EdgeBasis(8,:) = -EdgeBasis(8,:)
8411 EdgeBasis(9,1) = 0.0d0
8412 EdgeBasis(9,2) = 0.0d0
8413 EdgeBasis(9,3) = (v*(-Sqrt(3.0d0) + 2*v))/6.0d0
8418 EdgeBasis(9,:) = -EdgeBasis(9,:)
8454 EdgeBasis(10,:) = D1 * WorkBasis(I1,:)
8456 EdgeBasis(11,:) = D2 * WorkBasis(I2,:)
8484 EdgeBasis(12,:) = D1 * WorkBasis(I1,:)
8486 EdgeBasis(13,:) = D2 * WorkBasis(I2,:)
8514 EdgeBasis(14,:) = D1 * WorkBasis(I1,:)
8516 EdgeBasis(15,:) = D2 * WorkBasis(I2,:)
8545 EdgeBasis(2*(k-1)+1,1) = 0.5d0 * LineNodalPBasis(1,w) * LineNodalPBasis(q,v)
8548 EdgeBasis(2*(k-1)+2,1) = 1.5d0 * LineNodalPBasis(1,w) * u * LineNodalPBasis(q,v)
8552 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8567 EdgeBasis(2*(k-1)+1,1) = 0.5d0 * LineNodalPBasis(2,w) * LineNodalPBasis(q,v)
8570 EdgeBasis(2*(k-1)+2,1) = 1.5d0 * LineNodalPBasis(2,w) * u * LineNodalPBasis(q,v)
8574 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8589 EdgeBasis(2*(k-1)+1,2) = 0.5d0 * LineNodalPBasis(1,w) * LineNodalPBasis(3-q,u)
8592 EdgeBasis(2*(k-1)+2,2) = 1.5d0 * LineNodalPBasis(1,w) * v * LineNodalPBasis(3-q,u)
8596 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8611 EdgeBasis(2*(k-1)+1,2) = 0.5d0 * LineNodalPBasis(2,w) * LineNodalPBasis(3-q,u)
8614 EdgeBasis(2*(k-1)+2,2) = 1.5d0 * LineNodalPBasis(2,w) * v * LineNodalPBasis(3-q,u)
8618 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8633 EdgeBasis(2*(k-1)+1,3) = 0.5d0 * LineNodalPBasis(1,u) * LineNodalPBasis(q,v)
8636 EdgeBasis(2*(k-1)+2,3) = 1.5d0 * LineNodalPBasis(1,u) * w * LineNodalPBasis(q,v)
8640 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8655 EdgeBasis(2*(k-1)+1,3) = 0.5d0 * LineNodalPBasis(2,u) * LineNodalPBasis(q,v)
8658 EdgeBasis(2*(k-1)+2,3) = 1.5d0 * LineNodalPBasis(2,u) * w * LineNodalPBasis(q,v)
8662 EdgeBasis(2*(k-1)+1,:) = -EdgeBasis(2*(k-1)+1,:)
8710 EdgeBasis(k+4*(q-1)+1,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8712 EdgeBasis(k+4*(q-1)+2,:) = WorkBasis(2*(I1-1)+2,:)
8714 EdgeBasis(k+4*(q-1)+3,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8716 EdgeBasis(k+4*(q-1)+4,:) = WorkBasis(2*(I2-1)+2,:)
8760 EdgeBasis(k+1,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8762 EdgeBasis(k+2,:) = WorkBasis(2*(I1-1)+2,:)
8764 EdgeBasis(k+3,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8766 EdgeBasis(k+4,:) = WorkBasis(2*(I2-1)+2,:)
8810 EdgeBasis(k+1,:) = D1 * WorkBasis(2*(I1-1)+1,:)
8812 EdgeBasis(k+2,:) = WorkBasis(2*(I1-1)+2,:)
8814 EdgeBasis(k+3,:) = D2 * WorkBasis(2*(I2-1)+1,:)
8816 EdgeBasis(k+4,:) = WorkBasis(2*(I2-1)+2,:)
8822 EdgeBasis(49,1) = 8.0d0 * LineNodalPBasis(1,w) * LineNodalPBasis(2,w) * &
8827 EdgeBasis(50,1) = 24.0d0 * LineNodalPBasis(1,w) * LineNodalPBasis(2,w) * u * &
8833 EdgeBasis(51,2) = 8.0d0 * LineNodalPBasis(1,w) * LineNodalPBasis(2,w) * &
8838 EdgeBasis(52,2) = 24.0d0 * LineNodalPBasis(1,w) * LineNodalPBasis(2,w) * v * &
8843 EdgeBasis(53,3) = 8.0d0 * LineNodalPBasis(1,v) * LineNodalPBasis(2,v) * &
8848 EdgeBasis(54,3) = 24.0d0 * LineNodalPBasis(1,v) * LineNodalPBasis(2,v) * w * &
8866 EdgeBasis(1,1) = ((-1.0d0 + v)*v*(-1.0d0 + w)*w)/8.0d0
8867 EdgeBasis(1,2) = 0.0d0
8868 EdgeBasis(1,3) = 0.0d0
8873 EdgeBasis(1,:) = -EdgeBasis(1,:)
8883 EdgeBasis(2,1) = 0.0d0
8884 EdgeBasis(2,2) = (u*(1.0d0 + u)*(-1.0d0 + w)*w)/8.0d0
8885 EdgeBasis(2,3) = 0.0d0
8890 EdgeBasis(2,:) = -EdgeBasis(2,:)
8900 EdgeBasis(3,1) = (v*(1.0d0 + v)*(-1.0d0 + w)*w)/8.0d0
8901 EdgeBasis(3,2) = 0.0d0
8902 EdgeBasis(3,3) = 0.0d0
8907 EdgeBasis(3,:) = -EdgeBasis(3,:)
8917 EdgeBasis(4,1) = 0.0d0
8918 EdgeBasis(4,2) = ((-1.0d0 + u)*u*(-1.0d0 + w)*w)/8.0d0
8919 EdgeBasis(4,3) = 0.0d0
8924 EdgeBasis(4,:) = -EdgeBasis(4,:)
8934 EdgeBasis(5,1) = ((-1.0d0 + v)*v*w*(1.0d0 + w))/8.0d0
8935 EdgeBasis(5,2) = 0.0d0
8936 EdgeBasis(5,3) = 0.0d0
8941 EdgeBasis(5,:) = -EdgeBasis(5,:)
8951 EdgeBasis(6,1) = 0.0d0
8952 EdgeBasis(6,2) = (u*(1.0d0 + u)*w*(1.0d0 + w))/8.0d0
8953 EdgeBasis(6,3) = 0.0d0
8958 EdgeBasis(6,:) = -EdgeBasis(6,:)
8968 EdgeBasis(7,1) = (v*(1.0d0 + v)*w*(1.0d0 + w))/8.0d0
8969 EdgeBasis(7,2) = 0.0d0
8970 EdgeBasis(7,3) = 0.0d0
8975 EdgeBasis(7,:) = -EdgeBasis(7,:)
8985 EdgeBasis(8,1) = 0.0d0
8986 EdgeBasis(8,2) = ((-1.0d0 + u)*u*w*(1.0d0 + w))/8.0d0
8987 EdgeBasis(8,3) = 0.0d0
8992 EdgeBasis(8,:) = -EdgeBasis(8,:)
9002 EdgeBasis(9,1) = 0.0d0
9003 EdgeBasis(9,2) = 0.0d0
9004 EdgeBasis(9,3) = ((-1.0d0 + u)*u*(-1.0d0 + v)*v)/8.0d0
9009 EdgeBasis(9,:) = -EdgeBasis(9,:)
9019 EdgeBasis(10,1) = 0.0d0
9020 EdgeBasis(10,2) = 0.0d0
9021 EdgeBasis(10,3) = (u*(1.0d0 + u)*(-1.0d0 + v)*v)/8.0d0
9026 EdgeBasis(10,:) = -EdgeBasis(10,:)
9036 EdgeBasis(11,1) = 0.0d0
9037 EdgeBasis(11,2) = 0.0d0
9038 EdgeBasis(11,3) = (u*(1.0d0 + u)*v*(1.0d0 + v))/8.0d0
9043 EdgeBasis(11,:) = -EdgeBasis(11,:)
9053 EdgeBasis(12,1) = 0.0d0
9054 EdgeBasis(12,2) = 0.0d0
9055 EdgeBasis(12,3) = ((-1.0d0 + u)*u*v*(1.0d0 + v))/8.0d0
9060 EdgeBasis(12,:) = -EdgeBasis(12,:)
9099 EdgeBasis(13,:) = D1 * WorkBasis(I1,:)
9101 EdgeBasis(14,:) = D2 * WorkBasis(I2,:)
9129 EdgeBasis(15,:) = D1 * WorkBasis(I1,:)
9131 EdgeBasis(16,:) = D2 * WorkBasis(I2,:)
9159 EdgeBasis(17,:) = D1 * WorkBasis(I1,:)
9161 EdgeBasis(18,:) = D2 * WorkBasis(I2,:)
9189 EdgeBasis(19,:) = D1 * WorkBasis(I1,:)
9191 EdgeBasis(20,:) = D2 * WorkBasis(I2,:)
9219 EdgeBasis(21,:) = D1 * WorkBasis(I1,:)
9221 EdgeBasis(22,:) = D2 * WorkBasis(I2,:)
9249 EdgeBasis(23,:) = D1 * WorkBasis(I1,:)
9251 EdgeBasis(24,:) = D2 * WorkBasis(I2,:)
9257 EdgeBasis(25,1) = ((-1.0d0 + v**2)*(-1.0d0 + w**2))/2.0d0
9258 EdgeBasis(25,2) = 0.0d0
9259 EdgeBasis(25,3) = 0.0d0
9264 EdgeBasis(26,1) = 0.0d0
9265 EdgeBasis(26,2) = ((-1.0d0 + u**2)*(-1.0d0 + w**2))/2.0d0
9266 EdgeBasis(26,3) = 0.0d0
9271 EdgeBasis(27,1) = 0.0d0
9272 EdgeBasis(27,2) = 0.0d0
9273 EdgeBasis(27,3) = ((-1.0d0 + u**2)*(-1.0d0 + v**2))/2.0d0
9292 B(k) = SUM( LG(k,1:dim) * EdgeBasis(j,1:dim) )
9294 EdgeBasis(j,1:dim) = B(1:dim)
9345 WorkBasis(1,1:2) = EdgeBasis(j,1:2)
9346 EdgeBasis(j,1) = WorkBasis(1,2)
9347 EdgeBasis(j,2) = -WorkBasis(1,1)
9352 B(k) = SUM( LF(k,1:dim) * EdgeBasis(j,1:dim) ) / DetJ
9354 EdgeBasis(j,1:cdim) = B(1:cdim)
9361 B(k) = SUM( LG(k,1:dim) * EdgeBasis(j,1:dim) )
9363 EdgeBasis(j,1:cdim) = B(1:cdim)
9814 EdgeBasis(8,3), CurlBasis(8,3) local
9864 EdgeBasis(1,1) = (1.0d0 - v - w)/4.0d0
9865 EdgeBasis(1,2) = 0.0d0
9866 EdgeBasis(1,3) = (u*(-1.0d0 + v + w))/(4.0d0*(-1.0d0 + w))
9871 EdgeBasis(2,1) = 0.0d0
9872 EdgeBasis(2,2) = (1.0d0 + u - w)/4.0d0
9873 EdgeBasis(2,3) = (v*(1.0d0 + u - w))/(4.0d0 - 4.0d0*w)
9878 EdgeBasis(3,1) = (1.0d0 + v - w)/4.0d0
9879 EdgeBasis(3,2) = 0.0d0
9880 EdgeBasis(3,3) = (u*(1.0d0 + v - w))/(4.0d0 - 4.0d0*w)
9885 EdgeBasis(4,1) = 0.0d0
9886 EdgeBasis(4,2) = (1.0d0 - u - w)/4.0d0
9887 EdgeBasis(4,3) = (v*(-1.0d0 + u + w))/(4.0d0*(-1.0d0 + w))
9892 EdgeBasis(5,1) = (w*(-1.0d0 + v + w))/(4.0d0*(-1.0d0 + w))
9893 EdgeBasis(5,2) = (w*(-1.0d0 + u + w))/(4.0d0*(-1.0d0 + w))
9894 EdgeBasis(5,3) = (-((-1.0d0 + v)*(-1.0d0 + w)**2) + u*(v - (-1.0d0 + w)**2 - 2.0d0*v*w))/&
9900 EdgeBasis(6,1) = -(w*(-1.0d0 + v + w))/(4.0d0*(-1.0d0 + w))
9901 EdgeBasis(6,2) = (w*(-1.0d0 - u + w))/(4.0d0*(-1.0d0 + w))
9902 … EdgeBasis(6,3) = (-((-1.0d0 + v)*(-1.0d0 + w)**2) + u*((-1.0d0 + w)**2 + v*(-1.0d0 + 2.0d0*w)))/&
9908 EdgeBasis(7,1) = ((1.0d0 + v - w)*w)/(4.0d0*(-1.0d0 + w))
9909 EdgeBasis(7,2) = ((1.0d0 + u - w)*w)/(4.0d0*(-1.0d0 + w))
9910 EdgeBasis(7,3) = ((1.0d0 + v)*(-1.0d0 + w)**2 + u*(v + (-1.0d0 + w)**2 - 2.0d0*v*w))/&
9916 EdgeBasis(8,1) = (w*(-1.0d0 - v + w))/(4.0d0*(-1.0d0 + w))
9917 EdgeBasis(8,2) = -(w*(-1.0d0 + u + w))/(4.0d0*(-1.0d0 + w))
9918 EdgeBasis(8,3) = ((1.0d0 + v)*(-1.0d0 + w)**2 - u*(v + (-1.0d0 + w)**2 - 2.0d0*v*w))/&
9952 WBasis(i,k) = SUM( G(1:3,k) * EdgeBasis(i,1:3) )